Gutenberg Gait Database, a ground reaction force database of level overground walking in healthy individuals

The Gutenberg Gait Database comprises data of 350 healthy individuals recorded in our laboratory over the past seven years. The database contains ground reaction force (GRF) and center of pressure (COP) data of two consecutive steps measured - by two force plates embedded in the ground - during level overground walking at self-selected walking speed. The database includes participants of varying ages, from 11 to 64 years. For each participant, up to eight gait analysis sessions were recorded, with each session comprising at least eight gait trials. The database provides unprocessed (raw) and processed (ready-to-use) data, including three-dimensional GRF and two-dimensional COP signals during the stance phase. These data records offer new possibilities for future studies on human gait, e.g., the application as a reference set for the analysis of pathological gait patterns, or for automatic classification using machine learning. In the future, the database will be expanded continuously to obtain an even larger and well-balanced database with respect to age, sex, and other gait-specific factors. Measurement(s) ground reaction force • centre of pressure • walking behavior • gait measurement • Normal Gait Technology Type(s) force plate • Sensor Device Factor Type(s) age • sex • walking speed Sample Characteristic - Organism Homo sapiens Sample Characteristic - Environment laboratory environment Measurement(s) ground reaction force • centre of pressure • walking behavior • gait measurement • Normal Gait Technology Type(s) force plate • Sensor Device Factor Type(s) age • sex • walking speed Sample Characteristic - Organism Homo sapiens Sample Characteristic - Environment laboratory environment Machine-accessible metadata file describing the reported data: https://doi.org/10.6084/m9.figshare.15112488

classification tasks (e.g., healthy controls vs. Parkinson's disease). The question of whether it is possible to train ML models that meet clinical requirements in terms of robustness, transparency, and generalizability has rarely been investigated. This has so far hindered broader clinical application and acceptance of ML models. The availability of sufficient and high-quality data is an important prerequisite for the training of reliable ML models. However, the availability of 3DGA data is often a limitation in practice. Among other authors 24,25 , we also made 3DGA data available to the public [26][27][28][29][30] in previous studies 17, [31][32][33][34] . However, different data processing procedures and data structures were used in these studies, making collaborative use of the data difficult. In recent years, a rather small number of annotated large-scale datasets have been made publicly accessible 35 . Publicly available datasets can be used to train more robust models. In practice, even with a large-scale dataset, such as GaitRec 35 , data from individuals without gait pathology (healthy controls) represent a bottleneck. One reason for the scarcity of such data is that most gait analysis laboratories are located at clinics and usually record and examine only patients with pathological gait patterns.
In order to address this shortcoming, we provide -with the Gutenberg Gait Database -the gait data from healthy controls collected in our laboratory over the past seven years. The data is provided in a uniform format to allow for a continuously growing and publicly accessible database. The overall goal is to bridge existing gaps in publicly available gait datasets. Thereby, we aim at creating a basis for reliable ML models that can be used as decision-support system in clinical practice and research. Based on this goal, we prepared the processed data in such a way that it can be merged and used in conjunction with the GaitRec dataset 35 . In addition, the size and quality of the database allow it to serve as an extension of the study population in gait-related research areas, e.g., shoe and insole research 36 , security systems based on biometric recognition 37 , gait-based fatigue 38 and emotion 39 detection in psychological and sport-related contexts. In this setting our database can be used in various ways, e.g., as reference data or as source for automatic outlier detection. From a more epistemological point of view, the continuously growing database will also allow increasing flexibility in dealing with much more diverse questions related to human gait. Questions concerning population-motivated research 40 , problems of specific groups 41 , or the complexity of individual case-oriented time series 31,32 will be put on a broader data-based foundation over time.
The Gutenberg Gait Database provides exclusively force plate data, namely ground reaction force (GRF) and center of pressure (COP) signals. The current best practice in clinical gait analysis describes a patient's gait using a combination of force plate data with kinematic and electromyographic data. However, kinematic and electromyographic data are prone to several difficulties, such as inconsistencies due to differences in anthropometric characteristics of participants, experience of investigators, measurement protocols, and laboratory settings [42][43][44] . This makes it more difficult to create a homogeneous, large-scale, and high-quality dataset compared to using less interference-prone data, such as GRF signals 45,46 . Therefore, the use of force plate data offers advantages for the development of ML models for gait analysis, although the provided information appears to be reduced in comparison to kinematic data. However, previous studies 23,47 investigating ML methods for automated classification of gait impairments based on force plate data showed promising results suggesting their suitability for clinical applications.

Methods
Datasets. The Gutenberg Gait Database combines datasets from five already published studies on human gait 17,31-34 and data from five unpublished studies. A total sample of 350 participants (142 female, 205 male, and 3 unknown) aged between 11 and 64 years is included. Prior to the recording, all participants reported that they did not have any gait pathology and were not suffering from any injuries or diseases that affected gait. Table 1 summarizes demographic details for each individual dataset and the total database. Figure 1 shows the overall and sex-specific distributions of age, body mass, body height, and walking speed for the database.
All studies (published and unpublished) were carried out according to the Declaration of Helsinki at the Johannes Gutenberg-University in Mainz (Germany). All participants were informed about the experimental protocol and provided their written informed consent to participate in the study. The approval from the ethical committee of the medical association Rhineland-Palatinate in Mainz (Germany) was received.  www.nature.com/scientificdata www.nature.com/scientificdata/ Data recording & Experimental protocol. Bi-lateral analog force plate signals were recorded by asking participants to walk at their preferred (self-selected) walking speed on a level and approximately 10 m long walkway. Two force plate configurations were used: (i) an inline configuration using two centrally embedded force plates (Kistler, Type 9287CA, Switzerland) and (ii) a staggered configuration using two force plates (Kistler, Type 9286AA, Switzerland) integrated in a wooden walkway.
For both force plate configurations, the analog force plate signals were amplified (Kistler, Type 5233 A, Switzerland) and converted to digital signals using a sampling frequency of 1,000 Hz. A data acquisition system (Kistler, Type 5695, Switzerland) with a 16-bit analog-digital converter (Measurement Computing Corporation, Type USB-2533, USA) was used with a signal input range of ±10 V. Depending on the underlying experimental protocol, the walking speed was either estimated using (i) two light barriers with two photoelectric sensors (Imhof Timing, Germany) at a sampling frequency of 1,000 Hz or (ii) the three-dimensional pelvis marker trajectories captured by nine infrared cameras (Qualisys AB, Type Oqus 310, Sweden) at a sampling frequency of 250 Hz.
Participants were asked to perform gait trials to familiarize with the experimental setup and to determine an individual starting position for the gait analysis session. The number of familiarization trials differed between the experimental protocols. The exact number is specified for each study in Table 2. This procedure has already been shown to minimize the impact of targeting the force plates on the observed gait variables 48,49 . In addition, the participants were instructed to look at a symbol (neutral smiley) on the opposing wall of the laboratory to direct their attention away from the force plates and ensure a natural walk with an upright body position.
During one gait analysis session, participants walked until a predefined number of valid gait trials were available. These gait trials were defined as valid by the assessor if the participant walked "naturally" (e.g., with respect to force plate targeting) and both force plates were hit cleanly. The predefined number of gait trials per session varied between the experimental protocols and ranged from 8 to 20 gait trials. The exact number for each experimental protocol is specified in Table 2. Depending on the experimental study design, one to eight gait analysis sessions were recorded per participant. www.nature.com/scientificdata www.nature.com/scientificdata/ Data processing. The three-dimensional GRFs (vertical, anterior-posterior, and medio-lateral) and the two-dimensional COPs (anterior-posterior and medio-lateral) were calculated based on the analog force plate signals. The database provides unprocessed (raw) and processed (ready-to-use) GRF and COP signals during the stance phase. The data processing procedure was coordinated with Horsak et al. 35 so that the processing of the data in the Gutenberg Gait Database is identical to the GaitRec dataset. Thereby, we were able to prevent the obstacles that often exist in practice when using different datasets jointly. The main benefit for the community is the combined use of both data sources. We have, thus, eliminated a major disadvantage of the GaitRec dataset, namely ensuring that the number of healthy control participants is no longer a bottleneck.
For both settings, i.e., unprocessed and processed data, following pre-processing steps were performed. The offset of each analog force plate signal was corrected using the mean value of the first ten frames. The analog force plate signals were down-sampled to 250 Hz. The orientation of the medio-lateral and anterior-posterior GRF and COP signals were unified. Thus, medial and anterior forces were transformed to positive and lateral and posterior to negative values.
For the unprocessed (raw) data, we determined the signals in the following way. The stance phase was determined using a vertical GRF threshold of 25 N. The cropped GRF signals of the stance phase were used to calculate the COP signals.  Table 3. Description of the data stored in the "GRF_*.csv" files. "*" for the associated file name is a placeholder for "right" and "left" (adapted from Horsak et al. 35 ). n is either the number of frames during one step across the force plate for the unprocessed data ("RAW") or a time-normalized vector of 101 points for the processed ("PRO") data. Note that the first four columns of each file hold the DATASET_ID, SUBJECT_ID, SESSION_ID, and TRIAL_ID. www.nature.com/scientificdata www.nature.com/scientificdata/ For the processed (ready-to-use) data, we filtered the GRF signals using a second-order Butterworth bidirectional low-pass filter at a cut-off frequency of 20 Hz. The stance phase was determined based on the filtered GRF signals using a vertical GRF threshold of 25 N. For the processed COP signals, we filtered the unprocessed (raw) COP signals as well with a second-order Butterworth bidirectional low-pass filter at a cut-off frequency of 20 Hz. Furthermore, we cropped the filtered COP signals with a vertical GRF threshold of 80 N to avoid artifacts in COP calculation at small GRF signal values. In addition, the medio-lateral COP signals were mean-centered and anterior-posterior COP signals zero-centered. Each GRF and COP signal was time-normalized to 101 data points, corresponding to 100% stance phase. The GRF signals were normalized to the body weight, measured before each gait analysis session. The whole data processing was performed within the Matlab 2019a (The MathWorks, Inc., Natick, Massachusetts, USA) framework.

Data Records
All published data are fully anonymized and are available online from figshare 50 . As already pointed out, we decided to follow the data processing procedure and data structure as well as the naming of the files according to the GaitRec dataset 35 . The data records consist of twenty files containing the GRF data for each gait trial (see Table 3) and one file containing the measured walking speed for each gait trial. In addition, we provide one file containing metadata for each gait analysis session, including additional participants' information, e.g., class label, sex, age, body mass. All files are available as comma-separated value files (.csv). The twenty GRF data files are organized according to the following naming convention: "GRF-type-processing-side.csv". The type denotes, whether the file holds the vertical ("F_V"), anterior-posterior ("F_AP"), medio-lateral ("F_ML") or the anterior-posterior or medio-lateral COP ("COP_AP", "COP_ML") time-series. Processing denotes, if the files hold the unprocessed (raw) data ("RAW") or the processed (ready-to-use) data ("PRO"). The side denotes, if the data are from the "left" or "right" body side. The common prefix for all files is "GRF-". An example filename is: "GRF_F_V_RAW_left.csv".
Each of the "GRF-type-processing-side.csv" files is structured as a matrix with T rows × K columns (T = 8,819; K = 105 for "PRO" and K = 216 for "RAW"). Each row holds the data of one gait trial. The first column identifies www.nature.com/scientificdata www.nature.com/scientificdata/ each dataset ("DATASET_ID"), the second column each participant ("SUBJECT_ID"), the third column each gait analysis session ("SESSION_ID"), and the fourth column each single gait trial within a session ("TRIAL_ ID"). The remaining columns contain the values of the GRF signals for each gait trial. Note that due to the non-normalized nature of the data and the resulting different time-series lengths in the "RAW" files, non-available numbers have been replaced by "NaN" to maintain a constant matrix-dimension.
The file holding the measured walking speed for each gait trial is named "GRF_walking_speed.csv". The file is structured as a matrix with T rows × L columns (T = 8,819; L = 5). Each row holds the data of one gait trial. The first column identifies each dataset ("DATASET_ID"), the second column each participant ("SUBJECT_ ID"), the third column each gait analysis session ("SESSION_ID"), and the fourth column each single gait trial within a session ("TRIAL_ID"). The fifth column contains the measured walking speed for each gait trial ("WALKING_SPEED"). The walking speed was not measured in datasets 6, 8, and 10. Non-available numbers have been replaced by "NaN" to maintain a constant matrix-dimension.
The metadata file, which contains additional participant and session-related information is named "GRF_ metadata.csv" (see Table 4). The file is structured as a matrix with S rows × M columns (S = 661; M = 21). Here, the first three columns hold the DATASET_ID, SUBJECT_ID, and SESSION_ID, the other columns hold information such as sex, body mass, and age (see Table 4 for more details). Non-available numbers have been replaced by "NaN" to maintain a constant matrix-dimension.

technical Validation
The force plates and the measurement equipment were calibrated by the manufacturer (Kistler, Switzerland) and regularly checked and serviced during laboratory practice. No specific procedure (e.g., such as the CalTester method) was used.
In addition, on each day when measurements were conducted, the proper functioning of the force plates and measuring equipment was ensured by the following procedure: (i) A 30 s recording without load on the force plates was taken and ensured that the signal noise was below ±1 N. (ii) The assessor performed a weight measurement to verify the proper amplification of the analog channels. (iii) The assessor walked along the 10 m analysis walkway with one foot contact on each force plate and verified that the GRF signals showed the characteristic curves.
For an impression of data integrity, the processed data for each dataset is shown in Fig. 2 Table 4. Description of the information stored in the metadata file (adapted from Horsak et al. 35 ). *The metadata items highlighted by an asterisk were included primarily to ensure a consistent data structure between the Gutenberg Gait Database and the GaitRec dataset 35 . www.nature.com/scientificdata www.nature.com/scientificdata/

Usage Notes
The data are stored in *.csv files and can be easily imported into any software framework for further data analysis. We provide two scripts that allow a straightforward data import for Matlab (The MathWorks, Inc., Natick, Massachusetts, United States, 2019a) and Python (Python Software Foundation, 3.7). Additionally, two scripts (for Matlab and Python) are available for merging the GaitRec dataset 35 and the Gutenberg Gait Database. For the GaitRec dataset the DATASET_ID is set to 0. Since the metadata files and the data files have the same structure, a simple consolidation can be achieved. The GaitRec dataset has a bottleneck in terms of healthy control participants. Merging the two datasets can compensate for this limitation and allow the data to be much more useful for future research. Merging the two data sources would increase the number of healthy controls from 211 to 561, which approximately corresponds to the cardinality of the gait disorder classes: hip (N = 450), knee (N = 625), ankle (N = 627), calcaneus (N = 382).  Fig. 3 Visualization of zero-centered anterior-posterior (left panel) and mean-centered medio-lateral (right panel) components of the COP measurements per dataset. Mean and standard deviation signals (calculated per dataset) are highlighted as solid and dashed colored lines. We carefully inspected the gait trials where the signals differed considerably and made sure that these differences were not the result of measurement or calculation errors. Using the kinematic data, we were able to verify that the deviating signals were from gait trials of forefoot or midfoot walking participants.