Most digital models of the equine distal limb that are available in the community are static and/or subject specific; hence, they have limited applications in veterinary research. In this paper, we present an articulatable model of the entire equine distal limb based on statistical shape modeling. The model describes the inter-subject variability in bone geometry while maintaining proper jointspace distances to support model articulation toward different poses. Shape variation modes are explained in terms of common biometrics in order to ease model interpretation from a veterinary point of view. The model is publicly available through a graphical user interface (https://github.com/jvhoutte/equisim) in order to facilitate future digitalization in veterinary research, such as computer-aided designs, three-dimensional printing of bone implants, bone fracture risk assessment through finite element methods, and data registration and segmentation problems for clinical practices.

@en

Current classifications of midshaft clavicle fractures are based on radiography. The aim of the study was to evaluate the fracture pattern of clavicle fractures using 3-dimensional computed tomography (3D CT). A retrospective analysis was performed on CT scans of 65 acute clavicle fractures. Using quantitative 3D CT reconstruction techniques, the fracture of the clavicle was virtually reduced. Based on these reconstructions, a group-based fracture heat map and small fragment heat map, and the location of the most common fracture line were determined. Also, the direction and amount of displacement were evaluated. Three fracture patterns could be distinguished. The primary fracture line in type 1 is going from posteromedial to anterolateral and located between 50% and 68% of the clavicle’s length. In type 2, a transverse fracture line is located around 55%, and in type 3, a superolateral to inferomedial line is located between 47% and 56%. Wedged fracture fragments can be seen in types 1 and 2 and are mainly situated inferiorly. The displacement is similar in all types, but the main direction of displacement is specific for the different types (posterior, anterior, inferior). We can conclude that several fracture patterns can be seen in clavicle fractures. Most fractures are located laterally at the midshaft of the clavicle. Wedged segments are mainly located inferiorly, and at the posterior part of the clavicle, no comminution is ever seen. The direction of displacement depends on the fracture pattern.@en

In this paper, we investigate the heating function of the nasal cavity qualitatively, using a high-quality, large-scale statistical shape model. This model consists of a symmetrical and an asymmetrical part and provides a new and unique way of examining changes in nasal heating function resulting from natural variations in nasal shape (as obtained from 100 clinical CT scans). Data collected from patients suffering from different nasal or sinus-related complaints are included. Parameterized models allow us to investigate the effect of continuous deviations in shape from the mean nasal cavity. This approach also enables us to avoid many of the compounded effects on flow and heat exchange, which one would encounter when comparing different patient-specific models. The effects of global size, size-related features, and turbinate size are investigated using the symmetrical shape model. The asymmetrical model is used to investigate different types of septal deviation using Mladina's classification. The qualitative results are discussed and compared with findings from the existing literature.@en

The high prevalence of foot pain, and its relation to foot shape, indicates the need for an expert system to identify foot shape abnormalities. Yet, to date, no such expert system exists that examines the full 3D foot shape and produces an intuitive explanation of why a foot is abnormal. In this work, we present the first such expert system that satisfies those goals. The system is based on the concept of model-based outlier detection: a statistical shape model (SSM) is generated from 186 3D optical foot scans of healthy feet. This model acts as a knowledge base which is then parameterized by one's demographic characteristics (e.g., age, weight, height, shoe size) through a multivariate regression. This regression introduces model flexibility as it allows the model to be fine tuned to a specific individual. This fine tuned model is then used as a baseline to which the individual's 3D foot scan can be compared using standard statistical tests (e.g. t-tests). These statistical tests are performed at each vertex along the foot surface to identify the degree and location of shape outliers. Our expert system was validated on foot scans from patients with hallux valgus and abnormal foot arches. As expected, our results varied per patient, confirming that feet with the same clinical classification still show high shape variability. Additionally, the foot shape abnormalities identified by our system not only agreed with the expected location and severity of the tested foot deformities, but our analysis of the full 3D foot shape was able to completely characterize the extent of those abnormalities for the first time. These results show that the combination of statistical shape modelling, multivariate regression, and statistical testing is powerful enough to perform outlier detection for 3D foot shapes. The resulting insights provided by this system could prove useful in both shoe design and clinical diagnosis.@en

In times of online shopping, it is a challenge to select the right size of the desired clothing without fitting it before ordering. Therefore, this study describes three techniques to predict 3D upper body dimensions. The first method used basic personal parameters (gender, age, weight and length), the second technique used also the shoulder width and the last method used a 3D Styku scan to add extra input parameters. The accuracy of the three prediction methods was compared against hand measurements for 17 upper body dimensions of 37 subjects. The Intraclass Correlation Coefficient increases with 11.2% for the Styku method compared to the other methods. For chest, hip and waist measurements, the basic method is reliable to predict 3D body dimensions and indicate the right size from an existing collection. For more accurate upper body dimensions as needed for producing custom made clothing, a 3D Styku scan can supply the desired input.

@en

Data reduction techniques are commonly applied to dynamic plantar pressure measurements, often prior to the measurement's analysis. In performing these data reductions, information is discarded from the measurement before it can be evaluated, leading to unkonwn consequences. In this study, we aim to provide the first assessment of what impact data reduction techniques have on plantar pressure measurements. Specifically, we quantify the extent to which information of any kind is discarded when performing common data reductions. Plantar pressure measurements were collected from 33 healthy controls, 8 Hallux Valgus patients, and 10 Metatarsalgia patients. Eleven common data reductions were then applied to the measurements, and the resulting datasets were compared to the original measurement in three ways. First, information theory was used to estimate the information content present in the original and reduced datasets. Second, principal component analysis was used to estimate the number of intrinsic dimensions present. Finally, a permutational multivariate ANOVA was performed to evaluate the significance of group differences between the healthy controls, Hallux Valgus, and Metatarsalgia groups. The evaluated data reductions showed a minimum of 99.1% loss in information content and losses of dimensionality between 20.8% and 83.3%. Significant group differences were also lost after each of the 11 data reductions (α=0.05), but these results may differ for other patient groups (especially those with highly-deformed footprints) or other region of interest definitions. Nevertheless, the existence of these results suggest that the diagnostic content of dynamic plantar pressure measurements is yet to be fully exploited.

@en

For product developers that design near-body products, virtual mannequins that represent realistic body shapes, are valuable tools. With statistical shape modelling, the variability of such body shapes can be described. Shape variation captured by statistical shape models (SSMs) is often polluted by posture variations, leading to less compact models. In this paper, we propose a framework that has low computational complexity to build a posture invariant SSM, by capturing and correcting the posture of an instance. The posture-normalised SSM is shown to be substantially more compact than the non-posture-normalised SSM. Practitioner summary: Statistical shape modelling is a technique to map out the variability of (body) shapes. This variability is often polluted by variations in posture. In this paper, we propose a framework to build a posture invariant statistical shape model. Abbreviations: SSM: statistical shape model; 1D: one-dimensional; 3D: three-dimensional; DHM: digital human model; LBS: linear blend skinning; PCA: princial component analysis; PC: principal component; TTR: thumb tip reach.

@en

Multi-post-labeling-delay pseudo-continuous arterial spin labeling (multi-PLD PCASL) allows for absolute quantification of the cerebral blood flow (CBF) as well as the arterial transit time (ATT). Estimating these perfusion parameters from multi-PLD PCASL data is a non-linear inverse problem, which is commonly tackled by fitting the single-compartment model (SCM) for PCASL, with CBF and ATT as free parameters. The longitudinal relaxation time of tissue T1t is an important parameter in this model, as it governs the decay of the perfusion signal entirely upon entry in the imaging voxel. Conventionally, T1t is fixed to a population average. This approach can cause CBF quantification errors, as T1t can vary significantly inter- and intra-subject. This study compares the impact on CBF quantification, in terms of accuracy and precision, of either fixing T1t, the conventional approach, or estimating it alongside CBF and ATT. It is shown that the conventional approach can cause a significant bias in CBF. Indeed, simulation experiments reveal that if T1t is fixed to a value that is 10% off its true value, this may already result in a bias of 15% in CBF. On the other hand, as is shown by both simulation and real data experiments, estimating T1t along with CBF and ATT results in a loss of CBF precision of the same order, even if the experiment design is optimized for the latter estimation problem. Simulation experiments suggest that an optimal balance between accuracy and precision of CBF estimation from multi-PLD PCASL data can be expected when using the two-parameter estimator with a fixed T1t value between population averages of T1t and the longitudinal relaxation time of blood T1b.@en

During a transcatheter aortic valve implantation, an axisymmetric implant is placed in an irregularly shaped aortic root. Implanting an incorrect size can cause complications such as leakage of blood alongside or through the implant. The aim of this study was to construct a method that determines the optimal size of the implant based on the three-dimensional shape of the aortic root. Based on the pre-interventional computed tomography scan of 89 patients, a statistical shape model of their aortic root was constructed. The weights associated with the principal components and the volume of calcification in the aortic valve were used as parameters in a classification algorithm. The classification algorithm was trained using the patients with no or mild leakage after their intervention. Subsequently, the algorithms were applied to the patients with moderate to severe leakage. Cross validation showed that a random forest classifier assigned the same size in 65 ± 7% of the training cases, while 57 ± 8% of the patients with moderate to severe leakage were assigned a different size. This initial study showed that this semi-automatic method has the potential to correctly assign an implant size. Further research is required to assess whether the different size implants would improve the outcome of those patients. @en

This chapter presents digital human models (DHMs) as supporting tools in user-centered design for smart clothing. Smart clothing represents the new class of design era 2.0 with interactive technologies intended to be attractive, comfortable, and fit for the purpose of the identified user (Scataglini et al., 2019). While DHMs offer the possibility to analyze or predict human posture, motion, and other functions in a high-fidelity, physics-based, three-dimensional, real-time environment, smart clothing allows for monitoring vital (heart rate, breathing rate, temperature) and emotional (heart rate variability, electrodermal activity) parameters, movements, and postures. Effectiveness of functional wear is verified on the integration of DHMs into the design of a smart clothing system. Smart clothing represents a “second skin” that has a close “intimate” relation with the human body. This relation is physiological, psychological, biomechanical, and ergonomical. Physiological, biomechanical, and ergonomic aspects represent a retro-feedback of the user’s function in the iterative co-design workflow for designing smart garments, permitting the redesign and the technological refinement of it.

@en

In this paper, we present a new framework to integrate movement acquired by a motion capture system to a statistical body shape model using Blender. This provides a visualization of a digital human model based upon anthropometry and biomechanics of the subject. A moving statistical body shape model helps to visualize physical tasks with inter-individual variability in body shapes as well as anthropometric dimensions. This parametric modeling approach is useful for reliable prediction and simulation of the body shape movement of a specific population with a few given predictors such as stature, body mass index and age.@en

Parametric 3D human body models are valuable tools for ergonomic product design and statistical shape modelling (SSM) is a powerful technique to build realistic body models from a database of 3D scans. Like the underlying 3D scans, body models built from SSMs are typically represented with triangle meshes. Unfortunately, triangle meshes are not well supported by CAD software where spline geometry dominates. Therefore, we propose a methodology to convert databases of pre-corresponded triangle meshes into multi-patch B-spline SSMs. An evaluation on four 3D scan databases shows that our method is able to generate accurate and water-tight models while preserving inter-subject correspondences by construction. In addition, we demonstrate that such SSMs can be used to generate design manikins which can be readily used in SolidWorks for designing well conforming product parts.@en

Background: Foot morphology has received increasing attention from both biomechanics researches and footwear manufacturers. Usually, the morphology of the foot is quantified by 2D footprints. However, footprint quantification ignores the foot's vertical dimension and hence, does not allow accurate quantification of complex 3D foot shape. Methods: The shape variation of healthy 3D feet in a population of 31 adult women and 31 adult men who live in Belgium was studied using geometric morphometric methods. The effect of different factors such as sex, age, shoe size, frequency of sport activity, Body Mass Index (BMI), foot asymmetry, and foot loading on foot shape was investigated. Correlation between these factors and foot shape was examined using multivariate linear regression. Results: The complex nature of a foot's 3D shape leads to high variability in healthy populations. After normalizing for scale, the major axes of variation in foot morphology are (in order of decreasing variance): arch height, combined ball width and inter-toe distance, global foot width, hallux bone orientation (valgus-varus), foot type (e.g. Egyptian, Greek), and midfoot width. These first six modes of variation capture 92.59% of the total shape variation. Higher BMI results in increased ankle width, Achilles tendon width, heel width and a thicker forefoot along the dorsoplantar axis. Age was found to be associated with heel width, Achilles tendon width, toe height and hallux orientation. A bigger shoe size was found to be associated with a narrow Achilles tendon, a hallux varus, a narrow heel, heel expansion along the posterior direction, and a lower arch compared to smaller shoe size. Sex was found to be associated with differences in ankle width, Achilles tendon width, and heel width. Frequency of sport activity was associated with Achilles tendon width and toe height. Conclusion: A detailed analysis of the 3D foot shape, allowed by geometric morphometrics, provides insights in foot variations in three dimensions that can not be obtained from 2D footprints. These insights could be applied in various scientific disciplines, including orthotics and shoe design.

@en

Background

Pedobarography produces large sets of plantar pressure samples that are routinely subsampled (e.g. using regions of interest) or aggregated (e.g. center of pressure trajectories, peak pressure images) in order to simplify statistical analysis and provide intuitive clinical measures.

Research Question

We hypothesize that these data reductions discard gait information that can be used to differentiate between groups or conditions.


Methods

To test the hypothesis of null information loss, we created an implementation of statistical parametric mapping (SPM) for dynamic plantar pressure datasets (i.e. plantar pressure videos). Our SPM software framework brings all plantar pressure videos into anatomical and temporal correspondence, then performs statistical tests at each sampling location in space and time. Novelly, we introduce non-linear temporal registration into the framework in order to normalize for timing differences within the stance phase. We refer to our software framework as STAPP: spatiotemporal analysis of plantar pressure measurements. Using STAPP, we tested our hypothesis on plantar pressure videos from 33 healthy subjects walking at different speeds.

Results

As walking speed increased, STAPP was able to identify significant decreases in plantar pressure at mid-stance from the heel through the lateral forefoot. The extent of these plantar pressure decreases have not previously been observed using existing plantar pressure analysis techniques.

Significance

We therefore conclude that the subsampling of plantar pressure videos - a task which led to the discarding of gait information in our study - can be avoided using STAPP.
@en

In quantitative magnetic resonance T1 mapping, the Variable Flip Angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution T1 weighted images in a clinically feasible time. Fast, linear methods that estimate T1 maps from these weighted images have been proposed, such as DESPOT1 and iterative reweighted linear least squares (IRWLLS). More accurate, nonlinear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this work, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR T1 mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise T1 map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared to conventional gradient-based NLLS estimators, while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than efficient implementations of the VARPRO method. Furthermore, NOVIFAST is shown to be robust against initialization.@en

We present a new approach to estimate geometry parameters of glass fibers in glass fiber-reinforced polymers from simulated X-ray micro-computed tomography scans. Traditionally, these parameters are estimated using a multi-step procedure including image reconstruction, pre-processing, segmentation and analysis of features of interest. Each step in this chain introduces errors that propagate through the pipeline and impair the accuracy of the estimated parameters. In the approach presented in this paper, we reconstruct volumes from a low number of projection angles using an iterative reconstruction technique and then estimate position, direction and length of the contained fibers incorporating a priori knowledge about their shape, modeled as a geometric representation, which is then optimized. Using simulation experiments, we show that our method can estimate those representations even in presence of noisy data and only very few projection angles available.@en

The human nose is a complex organ that shows large morphological variations and has many important functions. However, the relation between shape and function is not yet fully understood. In this work, we present a high quality statistical shape model of the human nose based on clinical CT data of 46 patients. A technique based on cylindrical parametrization was used to create a correspondence between the nasal shapes of the population. Applying principal component analysis on these corresponded nasal cavities resulted in an average nasal geometry and geometrical variations, known as principal components, present in the population with a high precision. The analysis led to 46 principal components, which account for 95% of the total geometrical variation captured. These variations are first discussed qualitatively, and the effect on the average nasal shape of the first five principal components is visualized. Hereafter, by using this statistical shape model, two application examples that lead to quantitative data are shown: nasal shape in function of age and gender, and a morphometric analysis of different anatomical regions. Shape models, as the one presented here, can help to get a better understanding of nasal shape and variation, and their relationship with demographic data.@en

An important factor influencing the quality of magnetic resonance (MR) images is the reconstruction method that is employed, and specifically, the type of prior knowledge that is exploited during reconstruction. In this work, we introduce a new type of prior knowledge, partial discreteness (PD), where a small number of regions in the image are assumed to be homogeneous and can be well represented by a constant magnitude. In particular, we mathematically formalize the partial discreteness property based on a Gaussian Mixture Model (GMM) and derive a partial discreteness image representation that characterizes the salient features of partially discrete images: a constant intensity in homogeneous areas and texture in heterogeneous areas. The partial discreteness representation is then used to construct a novel prior dedicated to the reconstruction of partially discrete MR images. The strength of the proposed prior is demonstrated on various simulated and real k-space data-based experiments with partially discrete images. Results demonstrate that the PD algorithm performs competitively with state-of-the-art reconstruction methods, being flexible and easy to implement.

@en

Purpose: Quantitative T₁ mapping is a magnetic resonance imaging technique that estimates the spin-lattice relaxation time of tissues. Even though T₁ mapping has a broad range of potential applications, it is not routinely used in clinical practice as accurate and precise high resolution T₁ mapping requires infeasibly long acquisition times. Method: To improve the trade-off between the acquisition time, signal-to-noise ratio and spatial resolution, we acquire a set of low resolution T₁-weighted images and directly estimate a high resolution T₁ map by means of super-resolution reconstruction. Results: Simulation and in vivo experiments show an increased spatial resolution of the T₁ map, while preserving a high signal-to-noise ratio and short scan time. Moreover, the proposed method outperforms conventional estimation in terms of root-mean-square error. Conclusion: Super resolution T₁ estimation enables resolution enhancement in T₁ mapping with the use of standard (inversion recovery) T₁ acquisition sequences. Magn Reson Med, 2016.

@en

In quantitative MR T1 mapping, the spin-lattice relaxation time T1 of tissues is estimated from a series of T1-weighted images. As the T1 estimation is a voxel-wise estimation procedure, correct spatial alignment of the T1-weighted images is crucial. Conventionally, the T1-weighted images are first registered based on a general-purpose registration metric, after which the T1 map is estimated. However, as demonstrated in this paper, such a two-step approach leads to a bias in the final T1 map. In our work, instead of considering motion correction as a preprocessing step, we recover the motion-free T1 map using a unified estimation approach. In particular, we propose a unified framework where the motion parameters and the T1 map are simultaneously estimated with a Maximum Likelihood (ML) estimator. With our framework, the relaxation model, the motion model as well as the data statistics are jointly incorporated to provide substantially more accurate motion and T1 parameter estimates. Experiments with realistic Monte Carlo simulations show that the proposed unified ML framework outperforms the conventional two-step approach as well as state-of-the-art modelbased approaches, in terms of both motion and T1 map accuracy and mean-square error. Furthermore, the proposed method was additionally validated in a controlled experiment with real T1-weighted data and with two in vivo human brain T1-weighted data sets, showing its applicability in real-life scenarios.

@en