# Arc-Length Re-Parametrization and Signal Registration to Determine a Characteristic Average and Statistical Response Corridors of Biomechanical Data

- Department of Mechanical Engineering, University of Waterloo, Waterloo, ON, Canada

A characteristic average and biofidelity response corridors are commonly used to represent the average behaviour and variability of biomechanical signal data for analysis and comparison to surrogates such as anthropometric test devices and computational models. However, existing methods for computing the characteristic average and corresponding response corridors of experimental data are often customized to specific types or shapes of signal and therefore limited in general applicability. In addition, simple methods such as point-wise averaging can distort or misrepresent important features if signals are not well aligned and highly correlated. In this study, an improved method of computing the characteristic average and response corridors of a set of experimental signals is presented based on arc-length re-parameterization and signal registration. The proposed arc-length corridor method was applied to three literature datasets demonstrating a range of characteristics common to biomechanical data, such as monotonic increasing force-displacement responses with variability, oscillatory acceleration-time signals, and hysteretic load-unload data. The proposed method addresses two challenges in assessing experimental data: arc-length re-parameterization enables the assessment of complex-shaped signals, including hysteretic load-unload data, while signal registration aligned signal features such as peaks and valleys to prevent distortion when determining the characteristic average response. The arc-length corridor method was shown to compute the characteristic average and response corridors for a wide range of biomechanical data, while providing a consistent statistical framework to characterize variability in the data. The arc-length corridor method is provided to the community in the freely available and open-source software package, ARCGen.

## Introduction

Experimental data collected from biological tissues, volunteer studies, or post-mortem human subject testing is used to assess or validate surrogates, such as computational models and anthropomorphic test devices, but is often highly variable (Schmitt et al., 2014). This variability complicates the assessment and validation of biomechanical surrogates (Yoganandan et al., 2015), where an average response and statistical representation of variability are needed. Quantitative curve matching assessments such as r-squared, root mean square error, or more sophisticated techniques like dynamic time warping require a single representative or average curve to compare experimental data to surrogate models (Kim et al., 2013). In contrast, biofidelity response corridors have long been used to represent the variability of experimental data (e.g., CORA (Gehre et al., 2009)). However, the fundamental challenge of computing an average curve and corridors remains ill-defined and is often tailored to specific applications such that existing methods are not generally applicable.

Documented methods of computing a characteristic average and biomechanical response corridors have been broadly divided into time-based and cross-variable (representing one measured variable with respect to a second measured variable, such as force-displacement) methods (Kim et al., 2013). While all signals are, to one extent or another, time-based, the time component of a signal may not be applicable or useful, such as in the case of quasi-static mechanical characterization testing or in cases where time information may not be available, such as reproducing force-displacement responses from publications.

Due to the prevalence of time-series data in experimental work, time-based methods are among the most mature and widely used to determine characteristic curves from experimental data (Barker and Cronin, 2021). Time-based methods have generally relied on point-wise calculation of average values and statistics such as a point-wise standard deviation (Kim et al., 2013), requiring all experimental signals to be collected at a uniform rate and to have the same number of points. Unequally space data or signals of different sampling rates can be incorporated but must be resampled to match the sampling rate and the number of data points as the other signals being analyzed.

However, a challenge for time-based methods is that point-wise calculations can misrepresent or distort characteristics of the signals if critical features such as peaks and valleys are not aligned in time. Such misalignment could result from physical variation between individual experimental tests, or experimental difficulties such as synchronizing the start of time-series data between different tests. A range of methods have been presented to align signals based on techniques such as alignment of a characteristic time (Cavanaugh et al., 1986; Maltese et al., 2002) and alignment to minimize variation (Morgan et al., 1981; Eppinger et al., 1984) or maximized correlation between signals (Nusholtz et al., 2009; Donnelly et al., 2017), generally by shifting signals in time relative to one another. These phase-shift techniques have seen success in minimizing distortion during averaging but can become less effective if signals contain more than one critical feature (Ramsay, 2018).

Signal registration (also called curve registration) is a technique borrowed from the field of functional data analysis (Ramsay and Li, 1998; Tang and Muller, 2008; Gasser and Kneip, 2012; Ramsay, 2018) for aligning multiple features of time-series data simultaneously to compute a characteristic average. The principal behind signal registration is the introduction of warping functions,

Computing corridors for cross-variable signals presents a number of challenges. While researchers have successfully combined the average and standard deviation of two time-series signals to produce cross-variable response corridors (Cavanaugh et al., 1986), this method is only effective when signals are highly correlated and when the original sampling time data is available. Further, point-wise averages of signals based on one of the cross-variable signal axes can result in loss of the characteristic shape of the underlying data (Kim et al., 2013; Mattucci and Cronin, 2015) or entirely fail if signals do not extend to the same point on each axis (Lessley et al., 2004; Mattucci and Cronin, 2015), such as differences in displacement to failure.

A historical method of generating corridors for cross-variable data is the so-called ‘eyeball average’, wherein a characteristic average and corridors are drawn based on how the researcher interprets the signals (Lobdell et al., 1972). Unfortunately, this methodology is highly subjective and does not characterize the average curve or statistical variability of the data. The work of Lessley et al. (2004) presents one of the first rigorous methods for computing a characteristic average and corridors by normalizing each cross-variable signal against itself, then computing the average and standard deviation for each point in the normalized curves Shaw et al. (2006) later implemented an improved statistical representation of signals based on a two-dimensional normal distribution. However, both the methods were limited to monotonic data, requiring non-monotonic signals to be divided into strictly monotonic segments. Consistent segmentation can be challenging for highly variable signals or signals without a well-defined characteristic shape Perez-Rapela et al. (2018) eliminated the need for segmentation by re-parameterizing signals based on their arc-length, rather than ordinate and abscissa values. Perez-Rapela further utilized a simple method for aligning signal features by normalizing arc-length calculations to a feature shared between signals, such as a characteristic peak that appears in every signal.

While biomechanical data takes the forms of time-based or cross-variable signals, the overall shape and characteristics of the signals are equally important metrics as these factors place limitations on the methods used to analyze the signals. To that end, biomechanical data can be broadly divided into three categories based on signal shape. The first category includes strictly monotonic signals, such as the mechanical response of biological tissues such as ligament force-displacement signals (Mattucci et al., 2012; Mattucci and Cronin, 2015). Second, signals that exhibit oscillatory or otherwise non-monotonic behaviour in one axis. Acceleration-time signals, such as head acceleration under rapid deceleration being a classic example (Ewing and Thomas, 1972). Third, signals that exhibit hysteresis or non-monotonic behavior in both axes. The load-unload force-displacement response from impact testing, such pendulum impacts on the thorax of post-mortem human subjects (Kroell et al., 1971), fall under this category.

In the present study, an improved general methodology for computing a characteristic average and statistical response corridors from biomechanical data is presented to address many of the limitations encountered in previously documented methods. The proposed method incorporates arc-length re-parameterization and signal registration to enable general applicability, and was applied to three distinct signal types: monotonic increasing data (force-displacement) with variability, oscillatory acceleration-time data, and hysteretic load-unload force-displacement data.

## Methodology

The proposed method of computing a characteristic average and response corridor outputs from a set of input signals presented here can be divided into three stages: arc-length re-parameterization, signal registration, and statistical analysis (Figure 1). The analysis considers a set of ^{th} point of the ^{th} signal is defined as

**FIGURE 1**. An overview of the three stages of the arc-length corridor method used to compute the characteristic average and response corridors of a set of inputs signals.

### Stage 1: Arc-Length Re-Parameterization of Input Signals

In the first stage, individual signals were re-parametrized with respect to arc-length, a monotonically increasing metric intrinsically tied to the shape of each signal (Figure 2). Each signal was first scaled based on the total range of the data (Eqs 1, 2) such that the ordinate and abscissa values were of the same order in magnitude.

where

**FIGURE 2**. A graphical depiction of the arc-length re-parameterization process, where input signals **(A)** defined **(B)**.

In the scaling equations,

Following scaling, the arc-length of each ^{th} data point (^{th} signal was computed with zero arc length start corresponding to the first point of the signal. Linear behavior was assumed between discrete points of the signal to compute arc-length for simplicity and was found to be reasonable for the data sets considered in this study (Eq. 7). Higher-order interpolations can be applied in cases where the data is sparse. The arc-length corresponding to each point in the signal was then normalized using the total arc-length of each signal, where

### Stage 2: Signal Registration to Align Critical Features

Signal registration was applied to align features of all

**FIGURE 3**. A graphical depiction of the signal registration process to align critical features **(A)** Signals in the arc-length space may not have key features aligned. **(B)** Peaks occurring a different normalized arc-lengths **(C)** Warping functions are introduced for each signal **(D)** to align features. Signal registration operates on arc-length, not the underlying

Signal registration used of strictly monotonic, signal-specific warping functions,

The number of interior control points, m, required for a signal was problem-specific; however, some guidance was developed in the current study. It was found that signal registration was generally not required for monotonic data (Figure 4A). However, for non-monotonic data, the number of interior control points should be equal to the number of prominent inflection points present in the characteristic shape of the signal. This number should not include inflection points resulting from small-scale, high-frequency oscillations caused by noise or experimental variation (Figures 4B, C). In addition, some signals may demonstrate load-unload response requiring control points (Figure 4D) corresponding to inflection points in the arc-length space.

**FIGURE 4**. Exemplar signals to aid in selecting the number of control points (**(A)** Monotonic signal, no warping needed. **(B)** Non-monotonic, a single control point is used in this signal due to the single inflection point **(C)** Oscillatory signal, three control points corresponding to three inflection points. **(D)** Load-unload signal, one warping point is required for a single inflection point in

The location of the m control points in arc length, for all n signals were determined simultaneously by maximizing the cross-correlation score,

As C approaches unity, the cross-correlation between all signals improves. Correlation scores were computed for

In order to prevent shape distortion, warping was limited by a penalty function,

For the penalty factor

### Stage 3: Statistical Analysis of Re-Parametrized and Registered Signals

After arc-length re-parameterization and signal registration to ensure all signals have the same number of points and have critical features aligned, point-wise statistical analysis was conducted using normalized arc-length. This statistical analysis assumed all

The characteristic average of the signals was constructed by computing the mean of all

**FIGURE 5**. A graphical depiction of the statistical analysis of re-parameterized signals and construction of response corridors **(A)** The mean and standard deviation of signals are computed point-wise in the arc-length space for **(B)** Mean

From this distribution, a two-dimensional confidence region can be defined as

where

Response corridors were defined as the envelope of all elliptical confidence regions at each resampled arc-length value. As determining the envelope of a set of ellipses does not necessarily have a closed-form analytical solution, a marching squares algorithm (Lorensen and Cline, 1987) was used to extract this envelope numerically.

### Evaluation of Three Experimental Datasets Using the Arc-Length Corridor Method

Three data sets were assessed using the proposed arc-length corridor method. The first data set, comprising force-displacement data for cervical spine ligaments (Mattucci et al., 2012), demonstrated monotonic behavior with relatively large variability and a lack of shared endpoints. Five ligaments were tested, including the anterior longitudinal ligament (Figure 6), posterior longitudinal ligament, ligamentum flavum, capsular ligament, and interspinous ligament (Supplementary Material S1;Section 1). Force-displacement signals were reproduced from the original experimental data provided by the authors. In keeping with the method used by Mattucci et al., signals were truncated at the end of the traumatic region, excluding post-ultimate load response. Although Mattucci et al. also excluded the traumatic region of some experimental signals because the signals were considered uncharacteristic in shape, no portion of any signal was excluded in the current study. Since the signals were monotonic in nature, no control points were defined (

**FIGURE 6**. The quasi-static force-displacement response of 12 anterior longitudinal ligament specimens (Mattucci et al., 2012).

The second data set, head acceleration data from human volunteer 15 g frontal sled tests (Ewing and Thomas, 1972), was non-monotonic and oscillatory (Figure 7). This study reported head center of gravity kinematics in the global coordinate system, where the

**FIGURE 7**. Head acceleration **(A)** and displacement **(B)** in the *z*-direction for a 15 g frontal acceleration. Signals were collected from seven human volunteers.

The third dataset was force-displacement data from pendulum impacts on the thorax of post-mortem human subjects (Kroell et al., 1971), demonstrating non-monotonic and hysteretic response (Figure 8). Signals used in this work were digitized from the original (Kroell et al., 1971) publication with no additional processing aside from converting measurements to metric units.

**FIGURE 8**. Kroell et al. thoracic impact response signals collected from eight post-mortem human subjects.

## Results

The characteristic average and response corridors were computed for the three literature datasets to demonstrate the efficacy of the arc-length corridor method.

### Dataset 1: Monotonic Cervical Ligament Response Data

The arc-length corridor method was applied to five cervical ligament quasi-static experimental datasets (Mattucci and Cronin, 2015). Only the results of the Anterior Longitudinal ligament are documented here (Figure 9), with the remainder presented in the Supplemental Material S1;Supplementary Figures S1–5.

**FIGURE 9**. Anterior longitudinal ligament force-displacement responses superimposed with Mattucci average (blue) and arc-length average (black) and

In all cases, the arc-length corridor method produced a characteristic average that captured the three physiologic regions of ligament behavior: a toe region, a predominantly linear region, and a concave-down traumatic region. The characteristic average curve from the arc-length method agreed well with that calculated by Mattucci and Cronin (2015). The corridors computed by the arc-length method were defined as ±1 one standard deviation, and demonstrated broadening width with increased displacement.

### Dataset 2: Oscillatory Head Kinematics Data

The arc-length corridor method was applied to compute the characteristic average and response corridors for head acceleration (Figure 10) and displacement (Figure 11) in the *x*-direction and rotation about the ^{−2} was used during signal registration. Four warping control points were used to register the acceleration signals in the *z*-direction while two control points were used for displacement signals. A full accounting of warping control points is provided in Supplementary Table S1 for all kinematics. The calculated corridors were

**FIGURE 10**. Head z-acceleration for seven subjects, with the characteristic average curve (black) compared to the point-wise average (blue). Corridors are

**FIGURE 11**. Head z-displacement for seven subjects, with the characteristic average curve (black) compared to the point-wise average (blue). Corridors are

### Dataset 3: Hysteretic Thoracic Impact Response Data

Application of the arc-length corridor method to the Kroell dataset (Kroell et al., 1971) used two control points (

**FIGURE 12**. Kroell force-displacement signals overlaid with arc-length average (black) and corridors (gold) as well as the Lobdell corridors (blue). Arc-length corridors are

## Discussion

A consistent and reproducible method is needed to determine an average characteristic curve and quantify variability in biomechanical data that may be used as input for computational human models or to assess the output and biofidelity of a model. Such curves are used as input for tissue material properties (Mattucci and Cronin, 2015) in detailed human body models, and can be used to assess the biofidelity of human computational models using methods such as cross-correlation (Barker and Cronin, 2021). Contemporary data averaging and corridor methods are often tailored to a specific data set and require some judgement in application to other signals, such that no general, consistent methodology is available to address the wide variety of data observed in biological tests.

In a review of published methods of computing response corridors, Kim et al. (2013) determined that while most time- and cross-variable-based methods produce broadly similar characteristic averages, the computed response corridors differ to a large extent, likely owing to differences in formulation between each method. Kim et al. judged that a relatively simple time-based, point-wise averaging method produced the best corridors. However, it is important that Kim et al. used a reasonably simplistic synthetic dataset that did not include features unavoidable in experimental data such as noise, non-monotonic behaviour, or signal misalignment. The arc-length corridor method proposed in the current study employs arc-length re-parameterization to enable calculations for non-monotonic data, including load-unload type data without manual signal segmentation. Although differences in magnitude between the axes can skew arc-length calculations and misrepresent features, scaling of signals to compute arc-length minimized distortion. Such differences in magnitude are common in experimental data like force-displacement or stress-strain curves.

The integration of signal registration through the use of warping functions allowed for the automatic alignment of critical signal features (such as peaks and valleys) based on a user-defined number of control points and warping penalty factor. Signal registration addresses a limitation in other arc-length parameterization and point-wise averaging methods, which assume critical features occur at approximately the same value of the dependent variable in all signals. This is a valid assumption for signals that are monotonic in both axes or highly correlated signals that exhibit very little variability. However, this assumption is violated when input signals exhibit variability or may not share all the same features. Further, the utilization of signal registration allowed the characteristic average and response corridors to be computed without smearing or distorting the underlying shape of the data.

The proposed arc-length corridor method was used to analyze three datasets with distinctly varying signal characteristics, a generalization not possible with current methodologies. Monotonic force-displacement human ligament data (Mattucci and Cronin, 2015) demonstrating variability and lack of shared signal endpoints could not be assessed using point-wise averaging. Mattucci and Cronin (2015) proposed a tailored method based on the assumed shape of ligament behaviour. The data was analyzed by segmenting each signal into three characteristic shapes or regions, then normalizing and averaging curves within each region. One consequence of this method was an enforced stationary point at the end of the traumatic region. While this is a valid assumption from a physiologic perspective, such a trait was not readily apparent in the experimental signals.

In contrast, the arc-length corridor method did not contain a priori assumptions of signal shape. The method generated similar average curves to the Mattucci method, which is to be expected since the aim of the Mattucci method was to generate representative average curves tailored to the data set. However, the average curves did differ between methods for the interspinous (Supplementary Figure S5) and Ligamentum Flavum (Supplementary Figure S3). The arc-length corridor method resulted in lower average curves for the traumatic region of these two ligaments. The higher average produced by the Mattucci method may be the result of the assumed curve shape and magnitude difference for a small number of curves in each dataset. Importantly, the arc-length corridor method enabled the generation of statistical corridors, which was not possible using the method developed by Mattucci.

The head center of gravity data in the human volunteer sled data (Ewing and Thomas, 1972) were oscillatory in nature and demonstrated low variability. Applying arc-length re-parametrization without signal registration for the head z-acceleration, it was noted that the average and corridors did not capture the true mean height or location of the intermediate peaks centred at a time of approximately 115 m (Figure 13). These peaks were smeared out somewhat due to variations in arc-length between signals resulting from noise and experimental variability (Figure 13A). After performing signal registration with

**FIGURE 13**. Effect of signal registration on the characteristic average and response corridors for head acceleration in the *z*-direction. Signal registration **(B)** utilized **(A)**.

It was noted that the experimental signals exhibited divergent behaviour later in time (Figure 11) made signal registration difficult, where the differences in curve length and diminished cross-correlation between signals. Conversely, the arc-length corridor method produced characteristic averages and corridors that reflected experimental behaviour well in cases where signals were better correlated.

The low variability of the Ewing and Thomas dataset enabled point-wise methods (e.g., Barker and Cronin (2021)) to represent the average response and corridors. The arc-length corridor method produced broadly similar results compared with the time-based point-wise average and corridors computed for each set of signals. However, it is important to consider that this dataset is ideal for using a time-based point-wise average, as all signals were sampled at exactly the same times and published with key features well aligned in time. Not all time-based acceleration data reported in the literature is as well-behaved and highly correlated as the Ewing and Thomas data.

The Lobdell thorax impact response signals are a highly influential dataset in biomechanics, particularly in occupant safety, and serve as the basis of some of the earliest cross-variable response corridors (Lobdell et al., 1972). Of the three datasets examined, the application of the arc-length corridor method to the Kroell et al. thoracic impact dataset presented the largest improvement over other reported techniques. The arc-length corridor method provided the ability to compute a characteristic average and response corridors without segmentation of the signals (Lessley et al., 2004). Furthermore, the statistical underpinning of the corridors from the current methodology provided a significant improvement over the original

The arc-length corridor method presents several major advantages over existing techniques documented in the literature as a generalized methodology that can be applied across a wide range of biomechanical signals. While this methodology shares features from the arc-length method documented by Perez-Rapela et al. (2018), it is set apart by a more robust method of re-parameterizing input signals, the incorporation of signal registration to continuously align signal features, and the statistical framework and automated techniques used to extract response corridors. However, the arc-length corridor method does have some limitations. The method relies on the arc-length of a signal, such that high frequency oscillatory signals, such as those arising from noise, could skew results by increasing the computed arc-length of a signal. While signal registration does mitigate this issue to a degree, the arc-length method could produce skewed averages and corridors with extremely oscillatory or noisy input signals. Skewed results can be remedied in one of two ways: filtering the signals (recommended in most cases) or increasing the number of warping control points. In a similar vein, corridors could be skewed should one input signal be significantly longer than the others, such as if one signal in a set of time series responses was collected for twice as long as the others. The simplest solution to this issue is to crop the longest signal; however, cropping should only be applied when physically justified.

While signal registration was shown to be valuable in aligning input signals, especially as seen in the Ewing and Thomas data discussed earlier, signal registration could potentially introduce a degree of subjectivity, as the number of warping control points and penalty factor is somewhat specific to a given dataset. In this study, the number of control points was recommended to be equal to the number of prominent inflection points in the shape of the signal, and this rule of thumb worked well for the data sets considered in the current study. Furthermore, a penalty factor of

Figure 14 demonstrates how varying the penalty factor *.*

**FIGURE 14**. The effect of penalty factor *z*-direction. If the penalty factor is high **(A)**, signal registration is not performed. If the penalty factor is low **(C)**, signal registration can skew points towards certain features, producing areas of poor discretization and signal resolution. Appropriate penalty factor selection **(B)** produces a curve that captures signal features without significant changes in point density.

Finally, while the arc-length corridor technique can accommodate varying sampling rates, the original signal must have a sufficiently high sampling rate to ensure that no distortion is introduced during the re-parameterization and signal registration. Although linear interpolation between discrete signal points was applicable in the current study, if signals were not sampled at a sufficiently high rate, there is a possibility that the resampled signals may not reflect the underlying physics of the original signal. Re-parameterization with higher-order interpolation methods may help prevent distortion of the original signal in these cases. However, such techniques applied to sparsely sampled data present the risk of introducing artifacts if signals are not sufficiently sampled to resolve critical features.

While not a limitation, the arc-length corridor method is agnostic to the physical context of the data. To that end, anthropomorphic scaling, which is commonly applied to data collected from human volunteer or post-mortem human subject trials to account for variation between subject weight, height, and other dimensions, was not performed. Anthropomorphic scaling, such as the equal stress equal velocity method, impulse-momentum approach (Petitjean et al., 2015) or even more advanced techniques (Sun et al., 2016), would generally reduce variability between signals and result in tighter corridors. If required, anthropomorphic scaling can be considered prior to performing corridor generation using the methodology presented in this work. However, care must be taken to ensure that scaling is appropriately applied to avoid the risk of combining data from incongruous subjects. Additionally, while anthropomorphic scaling tends to reduce the variability of the dataset, it is important to note that the dataset no longer represents the subject population, but rather the variability with respect to a specific set of human dimensions.”

## Conclusion

This paper has presented a method of computing the characteristic average and bio-fidelity response corridors based on arc-length re-parameterization. Compared to the existing corridor creation techniques, the method presented in this paper has several advantages. First, re-parameterization of experimental signals allows for the analysis of data that may be hysteretic or may terminate at different locations. Second, variation and misalignment of signal features can be accommodated by using signal registration, preventing distortion or smearing of signal features, a common issue in point-wise averaging methods. Finally, this technique permits the ability to compute and interpret response corridors within a statistical framework when assessing the variability of a dataset.

Overall, the generalized nature of the proposed method allowed for consistent, statistically-based, and repeatable analysis of signals across many data types with a single analysis framework.

## Data Availability Statement

The datasets presented in this study, along with the code used to implement the arc-length corridor methodology ca be found at github.com/DCHartlen/ARCGen.

## Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

## Conflict of Interest

The 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, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Acknowledgments

The authors would like to thank the generous support of the Natural Science and Engineering Research Council of Canada and the Global Human Body Models Consortium.

## Supplementary Material

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

## References

Barker, J. B., and Cronin, D. S. (2021). Multilevel Validation of a Male Neck Finite Element Model with Active Musculature. *J. Biomechanical Eng.* 143, 1085501. doi:10.1115/1.4047866/1085501

Cavanaugh, J. M., Nyquist, G. W., Goldberg, S. J., and King, A. I. (1986). *Lower Abdominal Tolerance and Response*. SAE Technical Papers. SAE International. doi:10.4271/861878

Donnelly, B. R., Rhule, H. H., Moorhouse, K. M., Kang, Y. S., and Stammen, J. A. (2017). “An Improved Deflection Energy Method to Normalise PMHS Thoracic Response Data,” in Conference proceedings International Research Council on the Biomechanics of Injury, IRCOBI, Antwerp, Belgium. Zurich, Switzerland: International Research Council on the Biomechanics of Injury, 577–591.

Eppinger, R. H., Marcus, J. H., and Morgan, R. M. (1984). *Development of Dummy and Injury index for NHTSA's Thoracic Side Impact Protection Research Program*. SAE Technical Papers. SAE International. doi:10.4271/840885

Ewing, C. L., and Thomas, D. J. (1972). *Human Head and Neck Response to Impact Acceleration*. Pensacola, Fl: Naval Aerospace Medical RESEARCH Lab Pensacola Fl.

Gasser, T., and Kneip, A. (2012). Searching for Structure in Curve Samples. *J. Amer. Stat. Assoc.* 90(432), 1179-1188. doi:10.1080/01621459.1995.10476624

Gehre, C., Gades, H., and Wernicke, P. (2009). “Objective Rating of Signals Using Test and Simulation Responses,” in Proceedings: International Technical Conference on the Enhanced Safety of Vehicles, Stuttgart, Germany.

Kim, T., Shin, J., Ye, X., Crandall, J., Knospe, C., and Funk, J. (2013). Evaluation of Methods for the Development of Representative Responses and Corridors from Biomechanical Data Using Mechanical Models. *Int. J. Crashworthiness* 18, 633–646. doi:10.1080/13588265.2013.830946

Kroell, C. K., Schneider, D. C., and Nahum, A. M. (1971). *Impact Tolerance and Response of the Human Thorax*. SAE Technical Papers. SAE International. doi:10.4271/710851

Krumbholz, K., Hardy, A. J., and de Boer, J. (2020). Automated Extraction of Auditory Brainstem Response Latencies and Amplitudes by Means of Non-linear Curve Registration. *Comput. Methods Programs Biomed.* 196, 105595. doi:10.1016/J.CMPB.2020.105595

Lessley, D., Crandall, J., Shaw, G., Kent, R., and Funk, J. (2004). *A Normalization Technique for Developing Corridors from Individual Subject Responses*. SAE Technical Papers. SAE International. doi:10.4271/2004-01-0288

Lobdell, T. E., Kroell, C. K., Schneider, D. C., Hering, W. E., and Hahum, A. M. (1972). “Impact Response of the Human Thorax,” in *Human Impact Response: Measurement and Simulation*. Editors W. F. King, and H. J. Mertz (New York: Springer Science+Business Media, 201

Lorensen, W. E., and Cline, H. E. (1987). Marching Cubes: A High Resolution 3D Surface Construction Algorithm. *SIGGRAPH Comput. Graph.* 21, 163–169. doi:10.1145/37402.37422

Maltese, M. R., Eppinger, R. H., Rhule, H. H., Donnelly, B. R., Pintar, F. A., and Yoganandan, N. (2002). *Response Corridors of Human Surrogates in Lateral Impacts*. SAE Technical Papers. SAE International. doi:10.4271/2002-22-0017

Mattucci, S. F. E., and Cronin, D. S. (2015). A Method to Characterize Average Cervical Spine Ligament Response Based on Raw Data Sets for Implementation into Injury Biomechanics Models. *J. Mech. Behav. Biomed. Mater.* 41, 251–260. doi:10.1016/j.jmbbm.2014.09.023

Mattucci, S. F. E., Moulton, J. A., Chandrashekar, N., and Cronin, D. S. (2012). Strain Rate Dependent Properties of Younger Human Cervical Spine Ligaments. *J. Mech. Behav. Biomed. Mater.* 10, 216–226. doi:10.1016/j.jmbbm.2012.02.004

Morgan, R. M., Marcus, J. H., and Eppinger, R. H. (1981). *Correlation of Side Impact Dummy/Cadaver Tests*. SAE Technical Papers. SAE International. doi:10.4271/811008

National Highway Traffic Safety Administration (2017). Biomechanics Test Database. Available at: https://www.nhtsa.gov/research-data/databases-and-software (Accessed December 9, 2021).

Nusholtz, G. S., Hsu, T. P., Yibing, S., Kochekseraii, S. B., and Luna, M. A. G. (2009). “Creating Representative Curves from Multiple Time Histories of Vehicle, ATD and Biomechanics Tests,” in Proceedings of the 21st International Technical Conference on the Enhanced Safety of Vehicles, Stuttgart, Germany.

Perez-Rapela, D., Markusic, C., Whitcomb, B., Pipkorn, B., Forman, J. L., Donlon, J. P., et al. (2018). “Comparison of WorldSID to PMHS Kinematics in Far-Side Impact,” in Conference proceedings International Research Council on the Biomechanics of Injury, IRCOBI, Athens, Greece. International Research Council on the Biomechanics of Injury, 630. IRCOBI is centered in Zurich, Switzerland.

Petitjean, A., Trosseille, X., Yoganandan, N., and Pintar, F. A. (2015). “Normalization and Scaling for Human Response Corridors and Development of Injury Risk Curves,” in *Accidental Injury* (New York, NY: Springer New York), 769–792. doi:10.1007/978-1-4939-1732-7_26

Ramsay, J. (2018). “Curve Registration,” in *The Oxford Handbook of Functional Data Analysis*. Editors F. Ferraty, and Y. Romain (Oxford: Oxford University Press). doi:10.1093/oxfordhb/9780199568444.013.9

Ramsay, J. O., and Li, X. (1998). Curve Registration. *J. R. Stat. Soc. Ser. B (Statistical Methodology)* 60, 351–363. doi:10.1111/1467-9868.00129

Schmitt, K.-U., Niederer, P. F., Cronin, D. S., Muser, M. H., and Walz, F. (2014). Introduction. *Trauma Biomech. Introduction Inj. Biomech.* 9783642539206, 1–14. doi:10.1007/978-3-642-53920-6

Shaw, J. M., Herriott, R. G., McFadden, J. D., Donnelly, B. R., and Bolte, J. H. (2006). Oblique and Lateral Impact Response of the PMHS Thorax. *Stapp Car Crash J.* 50, 147–167. doi:10.4271/2006-22-0007

Sun, W., Jin, J. H., Reed, M. P., Gayzik, F. S., Danelson, K. A., Bass, C. R., et al. (2016). A Method for Developing Biomechanical Response Corridors Based on Principal Component Analysis. *J. Biomech.* 49, 3208–3215. doi:10.1016/J.JBIOMECH.2016.07.034

Tang, R., and Muller, H.-G. (2008). Pairwise Curve Synchronization for Functional Data. *Biometrika* 95, 875–889. doi:10.1093/biomet/asn047

Yoganandan, N., Yoganandan, N., and Nahum, A. M. (2015). *Accidental Injury: Accidental Injury: Biomechanics and Prevention, New York, NY*, 1–851. doi:10.1007/978-1-4939-1732-7

## Nomenclature

^{th} point of the horizontal and vertical components of signal

Keywords: arc-length re-parameterization, biomechanical data, statistical response corridors, characteristic average curve, biofidelity, force-displacement, acceleration-time

Citation: Hartlen DC and Cronin DS (2022) Arc-Length Re-Parametrization and Signal Registration to Determine a Characteristic Average and Statistical Response Corridors of Biomechanical Data. *Front. Bioeng. Biotechnol.* 10:843148. doi: 10.3389/fbioe.2022.843148

Received: 24 December 2021; Accepted: 23 February 2022;

Published: 24 March 2022.

Edited by:

Carlos D. Maciel, University of São Paulo, BrazilReviewed by:

Antonio Apicella, Università della Campania Luigi Vanvitelli, ItalyFrancisco J. Lopez-Valdes, Comillas Pontifical University, Spain

Giovana Yuko Nakashima, Federal Institute of São Paulo, Brazil

Copyright © 2022 Hartlen and Cronin. 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: Duane S. Cronin, duane.cronin@uwaterloo.ca