Skip to main content


Swipe om te navigeren naar een ander artikel

01-12-2018 | Methodology | Uitgave 1/2018 Open Access

Journal of Foot and Ankle Research 1/2018

A minimal markerset for three-dimensional foot function assessment: measuring navicular drop and drift under dynamic conditions

Journal of Foot and Ankle Research > Uitgave 1/2018
Patric Eichelberger, Angela Blasimann, Nicole Lutz, Fabian Krause, Heiner Baur
Belangrijke opmerkingen

Electronic supplementary material

The online version of this article (https://​doi.​org/​10.​1186/​s13047-018-0257-2) contains supplementary material, which is available to authorized users.
Skin surface marker calcaneus
Dynamic navicular drift
Dynamic navicular drop
Foot strike
Intraclass correlation coefficient
Limits of agreement
First assessment on day 1
Second assessment on day 1
Assessment on day 2
Skin surface marker 5th metatarsal head
Skin surface marker 1st metatarsal head
Skin surface marker navicular tuberosity
Navicular height
Navicular height at foot strike
Minimum navicular height during stance
Navicular height during static trial
Navicular width
Navicular width at foot strike
Maximum navicular width during stance
Navicular width during static trial
Root mean square error
Standard deviation
Standard deviation of the test-retest differences
Standard error of measurement
Stance phase


Consisting of numerous bones, joints and ligaments, the foot is one of the mechanically more complex structures in the human body [1]. As it is the link to the supporting surface, it is subjected to high impact forces and moments, and proper function is crucial for locomotion since the foot influences the whole kinematic chain [2]. The concept, classifying feet into planus, recuts and cavus based on morphological criteria [1], is still common in foot function evaluation. It is based on the common understanding of the ’excessively pronated’ foot being dysfunctional and the degree of pronation being classified based on structural measures like malleolar valgus index, arch height index or arch height flexibility [1]. Although foot problems, lower leg injury and lower extremity kinematics are often attributed to foot type, causality is not conclusive [37]. This is maybe because there is no clear consensus as to whether prediction of dynamic foot function based on structural criteria is reasonable [1, 5, 8], and that clinical measurement of foot pronation is inherently difficult owing to the complex interaction of the joints of the foot when this movement occurs [6, 9]. Furthermore, current methodological limitations in relation to measurement of dynamic foot pronation [5] may be another reason. Conventional measures, like maximum and minimum rearfoot eversion or eversion excursion, are thought to represent foot pronation. These parameters acquired at the rearfoot are likely to be limited for predicting a motion that mainly occurs at the midfoot. It was proposed to merge data on structure with information on foot function under dynamic loading to more closely relate to the functional behaviour of the foot during locomotion [8]. Although numerous foot models to study multi-segment foot kinematics are currently available [10] they do not come up with the requirements of clinical assessments such as being simple and efficient [11] owing to complexity [12] and the knowledge about their reliability [10] is limited.
The vertical navicular drop measures the difference in navicular height between a subtalar neutral and a weight-bearing position [13] and serves as a quick clinical tool for the assessment of a surrogate parameter of midfoot mobility. Similarly, the medial navicular drift was proposed to provide further insight into the mechanics of the talonavicular joint [6]. Both, vertical and mediolateral movements, should be assessed to understand the movement of the midfoot [11] and may be regarded as meaningful clinical indicators to represent overall foot function [6]. Measuring navicular mobility under dynamic conditions by the dynamic navicular drop and drift (dNDrop, dNDrift) seems to be a promising and intuitive approach to characterise foot function. Evidence shows that the navicular drop [13] is a poor predictor of the dynamic navicular drop [9, 1416] and that it is therefore necessary to measure the navicular drop dynamically in order to be representative for foot function. Studies that addressed the dynamic navicular drop used either two-dimensional (2D) video technique [14, 1620], single- or biplane X-ray [9, 21], three-dimensional (3D) motion capture [15] or a stretch sensor device [22, 23]. Although important to describe navicular mobility comprehensively, the dynamic navicular drift has only been examined by Cornwall and McPoil [24]. While some of the 2D studies concerned repeatability [18, 19], the only study that used 3D methods [15] did not consider any reliability issues and did not make use of the potential to measure dynamic navicular drift simultaneously to dynamic navicular drop. The authors recently presented a concept [25] to provide a practical tool to measure clinically accepted foot parameters under dynamic conditions. This study aimed to evaluate the intra- and interday reliability of measuring dynamic navicular drop and drift and the associated time points using the novel concept [25] during the stance phase of walking in healthy adults.


Participants and experimental procedure

The investigation was implemented as a descriptive, cross-sectional laboratory study analysing the intra- and interday reliability of dynamic navicular drop and drift measurements during walking gait. The study was approved by the ethics committee of the canton of Bern (KEK number Z007/12) and all participants signed informed consent. The examined sample consisted of twenty-one healthy participants, 12 females (foot length 239 (SD 11) mm, age 27.9 (SD 5.0) years, height 166 (SD 6) cm, body mass 62.3 (SD 6.9) kg) and nine males (foot length 266 (SD 12) mm, age 32.8 (SD 10.3) years, height 181 (SD 7) cm, body mass 75 (SD 7.4) kg). Exclusion criteria were back pain or other spinal disorders, any kind of musculoskeletal affections, neurological injuries or diseases, thrombosis or fractures in the past 12 months, implants or surgery on the lower extremity in the past 12 months, bone tumours, angiopathies, alcohol abuse, dementia, acute infections (such as common cold) or high-intensity exercise the day before the measurements.
Each subject underwent a test-retest study work flow, where the navicular height (NH) and width (NW) during the stance phase (SP) were assessed with a set of four markers (Fig. 1) on the left and right foot. A testing session consisted of ten barefoot walking trials and all participants completed two subsequent (no pause between) sessions on a first day (M1a, M1b) to assess intraday reliabiltiy and a third session (M2) one week apart from their first day to assess interday reliability. An optical motion capture system (cameras: 2× Bonita10, 8× Bonita3; volume: 5.5 m × 1.2 m × 2 m; 200 Hz; Nexus 1.8.5 software, Vicon Motion Systems Ltd., Oxford, UK) was used for kinematic measurements. Initial contact and toe off events were determined based on ground reaction force measurements with two force plates (AMTI OR 6, Watertown, USA, 1000 Hz, threshold 25 N).

Markerset and model outputs

Reflective skin surface markers (diameter 16 mm) were placed by a single experienced physical therapist on the left and the right foot to track the following anatomical landmarks (Fig. 1):
first Metatarso-Phalangeal joint, metatarsal head, dorso-Medial aspect.
fifth Metatarso-Phalangeal joint, metatarsal head, dorso-Lateral aspect.
CAlcaneal tuberosity, posterior aspect of the calcaneus (Achilles’ tendon insertion).
NAvicular tuberosity, most medial aspect of the navicular bone.
MPM, MPL and CA markers served to define the XY-plane (Fig. 1) which was calibrated to be parallel to the foot’s plantar surface based on a static trial and CA marked the origin of a Cartesian coordinate system. The foot’s longitudinal axis (Y) was the bisecting line [26] of the MPM, MPL and CA triangle (Fig. 1). According to the right-hand rule convention, the foot’s vertical axis (Z) was perpendicular to the XY-plane pointing towards cranial and the foot’s lateral axis (X) pointed laterally for the right foot and medially for the left foot. The navicular height described the distance of the navicular marker from the XY-plane and the navicular width described the medial distance of the navicular maker from the YZ-plane.

Post-processing and statistics

Post-processing and statistical analysis were performed using Matlab (Version R2016a, The MathWorks Inc., Natick, USA). Navicular height and width were calculated and low-pass filtered (4th order zero-lag Butterworth, 4 Hz). Stance phase intervals were extracted, resampled to 200 samples and subsequently averaged among ten trials. Navicular height and width from the static trials (NHS and NWS) were taken as reference positions to express the navicular height and width from gait trials. The curve features, dynamic navicular drop and drift, were extracted from the navicular height and width time-series, respectively, and calculated as differences between minimum (for navicular height, Eq. 1) or maximum (for navicular width, Eq. 2) excursions during the stance phase and values at foot strike (NHFS, NWFS):
$$\begin{array}{*{20}l} dNDrop & = NH_{\text{FS}} - NH_{\text{Min}} \end{array} $$
$$\begin{array}{*{20}l} dNDrift & = NW_{\text{FS}} - NW_{\text{Max}} \end{array} $$
Respective time points of dynamic navicular drop and drift were determined in %SP and denoted tdNDrop and tdNDrift. Reliability was primary analyzed by the Bland-Altman method [27] and completed by the intraclass correlation coefficient ICC21 and the standard error of measurement (SEM, Eq. 3) [28] with SDd being the standard deviation of the differences from the Bland-Altman analysis.
$$ SEM = \frac{SDd}{\sqrt{2}} $$
Left and right feet measurements were not treated separately but treated as independent samples making up a total of 42 feet that were assessed at the time points M1a, M1b and M2. Reliability of NHFSNHMinNWFSNWMaxtdNDrift and tdNDrop was also considered to explore how each parameter contributed to intra- and interday variations in dynamic navicular drop and drift, respectively. All variables and the Bland-Altman differences were tested for normal distribution by the Kolmogorow-Smirnow test. The significance of the bias from the Bland-Altman analysis was evaluated by paired t-tests. The level of statistical significance was set as p<0.05.

Instrumental error assessment

To estimate the instrumental errors, the marker set was applied to a ski boot, which served as a rigid frame. A person wearing the ski boot completed ten walking trials. Model outputs were calculated as described in the “Post-processing” section based on marker distances measured with a caliper. Under the assumption that the markers built a rigid cluster, all marker distances would have been constant throughout the stance phase and deviations from the rigid cluster model outputs were taken for instrumental error estimation by calculating the root mean square errors among the stance phase.


At all three testing sessions the proposed markerset delivered distinct movement patterns of navicular height and width throughout the stance phase among the 21 healthy individuals examined in the present study (Fig. 2a and c). Similar to the mean navicular height and width among all individuals, the standard deviations of 10 walking trials from each subject were averaged among all participants to describe the trial-to-trial errors for each testing session (Fig. 2b and d). Step-to-step variations were relatively constant throughout the stance phase and were between 0.65 and 1.17 mm for navicular height and between 0.38 and 0.88 mm for navicular width. Both, navicular height and width, showed pronounced minima and maxima, respectively, between 74.87 and 77.39% SP (Figs. 2a and c, 3b). A slow caudal drop and medial drift up to around 76% SP followed by a more rapid cranial rise and lateral shift until toe off could have been observed. Except the time point of navicular drop and drift at session M1b, all variables extracted from navicular height and width time-series followed a normal distribution. Complete numerical values of the descriptive statistics are given in the Tables 1, 2, 3 and 4. Related to the static pose, the mean navicular height at foot strike was between 0.9 and 1.2 mm (SD 3.0 to 3.5 mm) and the mean minimum navicular height during the stance phase between -4.8 and -5.1 mm (SD 2.5 to 2.7 mm) (Fig. 3). The resulting dynamic navicular drop was between -5.9 and -6 mm (SD 2.3 mm) (Fig. 3). The mean navicular width at foot strike was between -2.9 and -3.1 mm (SD 1.8 to 1.9 mm) and the mean maximum navicular width during the stance phase was 2.5 mm (SD 1.7 to 2.3 mm) (Fig. 3). The resulting dynamic navicular drift was between 5.4 and 5.6 mm (SD 2.3 to 2.7 mm) (Fig. 3).
Table 1
Dynamic navicular drop and associated time points during stance
dNDrop M1a
dNDrop M1b
dNDrop M2
tdNDrop M1a
tdNDrop M1b
tdNDrop M2
− 5.94
− 5.96
− 6.00
− 5.63
− 5.54
− 5.84
− 11.59
− 11.82
− 12.48
− 0.48
− 0.38
− 1.51
< 0.001
< 0.001
< 0.001
KSp denotes the p-value from the Kolmogorow-Smirnow normality distribution test
Table 2
Dynamic navicular drift and associated time points during stance
dNDrift M1a
dNDrift M1b
dNDrift M2
tdNDrift M1a
tdNDrift M1b
tdNDrift M2
KSp denotes the p-value from the Kolmogorow-Smirnow normality distribution test
Table 3
Navicular height and width at foot strike
− 3.05
− 3.06
− 2.88
− 2.79
− 3.22
− 2.76
− 9.76
− 9.05
− 6.49
− 8.83
− 8.91
− 6.40
KSp denotes the p-value from the Kolmogorow-Smirnow normality distribution test
Table 4
Minimum and maximum navicular height and width during stance
NH MinM1a
NH MinM1b
NH MinM2
NW MaxM1a
NW MaxM1b
NW MaxM2
− 5.08
− 4.76
− 4.85
− 4.85
− 4.48
− 4.46
− 10.24
− 9.97
− 10.46
− 1.33
− 0.55
− 0.53
KSp denotes the p-value from the Kolmogorow-Smirnow normality distribution test


Since the article focuses on measuring navicular drop and drift, Bland-Altman diagrams are presented for these outcomes only in the main part of the manuscript. For the Bland-Altman diagrams of the other variables the reader is referred to the Additional file 1. None of the considered variables had a significant intra- or interday bias.
The intra- and interday repeatability of the dynamic navicular drop were found to be 1.1 and 2.3 mm, those of the dynamic navicular drift 1.5 and 5.2 mm (Fig. 4), respectively. Associated intraclass correlation coefficients were found to be 0.97 and 0.83 for the dynamic navicular drop and 0.96 and 0.46 for the navicular drift (Fig. 4).
Intra- and interday repeatability of navicular height at foot strike were 2.7 mm (ICC21 0.91) and 4.9 mm (ICC21 0.68) and those of the minimum navicular height were 2.8 mm (ICC21 0.85) and 4.7 mm (ICC21 0.55). The intra- and interday repeatability of the navicular width at foot strike were 2.3 mm (ICC21 0.81) and 3.7 mm (ICC21 0.49) and those of the maximum navicular width were 2.7 mm (ICC21 0.72) and 4.7 mm (ICC21 0.35).
The time point of the dynamic navicular drop had an intra- and interday repeatability of 22% SP (ICC21 0.79) and 21% SP (ICC21 0.67), respectively. The time point of the dynamic navicular drift had an intra- and interday repeatability of 26% SP (ICC21 0.36) and 20% SP (ICC21 0.43), respectively.

Instrumental errors

The mean deviations from the reference, i.e. the instrumental errors, had some kind of systematic patters over the stance phase, especially the navicular width model output (Fig. 5). The instrumental errors were in minimum and maximum -0.38 mm and 0.48 mm, respectively, for navicular height and -1.40 mm and 0.67 mm, respectively, for navicular width (Fig. 5). Quantified as root mean square errors over the stance phase, the instrumental error was RMSENH = 0.25 mm for the navicular height and RMSENW = 0.86 mm for the navicular width (Fig. 5).


We developed a minimal markerset to dynamically assess foot function by means of navicular mobility and examined the intra- and interday reliability of the dynamic navicular drop and drift as measures for the maximum vertical and medial displacements of the navicular during the stance phase of walking. The investigation revealed navicular height stance phase patterns and dynamic navicular drop values that are consistent with current knowledge. The observed navicular height patterns resembled those previously presented by other authors [15, 22, 23], who also found the navicular height being decreasing after foot strike with a characteristic minimum around 75% SP, followed by a subsequent increase during the remaining stance phase. The mean dynamic navicular drop from this study (6.0 mm) was below those which Dicharry [15] found in hypomobile, neutral and hypermobile feet (7.9 to 8.3 mm) with a similar 3D measurement approach but were slightly higher or similar to those Nielsen et al. [17] measured with a 2D approach in supinated, neutral and pronated feet. Compared to values that were presented as reference data based on 79 healthy participants measured with 2D video (mean 5.4±1.7 mm, range [1.3,9.2] mm), the values from this study were on average slightly higher but exhibited an extended range (0.4 to 12.5 mm). Compared with the only reference for the dynamic navicular drift [24] (4.7±2.0 mm) known to the authors, the values from this study were slightly higher. We found the following stance phase pattern for the navicular: caudal drop and medial drift up to 75% SP followed by subsequent cranial rise and lateral drift. This complies with the common understanding of normal foot function with pronation and arch flattening during loading response and midstance followed by supination and arch rise during push-off to stiffen the foot for effective propulsion power generation [2]. Discrepancies in dynamic navicular drop may be explained by uncontrolled foot lengths which are thought to influence the dynamic navicular drop [18], inclusion of participants without classifying feet and the limited sample size of only 21 participants. Hence the presented data should not be regarded as representative normative data but the outlined measurements can be considered with well provided external validity. Further, as for every skin-surface marker based movement analysis protocol, whether 2D or 3D, soft tissue artifacts (STAs) should always be borne in mind in terms of internal validity. The navicular height has been shown to be overestimated of about 1-2 mm during the first half of stance but up to 12 mm during terminal stance [29]. Hence, the proposed protocol, as well as similar protocols already known from literature, do most likely not exactly measure the position of the navicular tuberosity. However, the displacement patterns seem to be well represented by skin-surface marker based approaches, which may deliver more clinically important information about dynamic foot function compared with static assessments.
For both, the dynamic navicular drop and drift, the interday reliability was worse than the intraday reliability: repeatability and SEM were larger and ICC21 values smaller. The ICC values indicate that the dynamic navicular drop can be reliably measured within and between days but that the dynamic navicular drift is only reliably measurable within day. However the conclusions based on ICCs have to be drawn carefully, because ICCs are dependent from the sample heterogeneity [28], i.e. the actual values of dynamic navicular drop and drift, respectively. The values reported as repeatability are absolute measures of measurement error and quantify the minimal detectable changes for within and between day assessments. Hence, changes in dynamic navicular drop must be in minimum 1.1 mm (19%) and 2.3 mm (38%) to be detectable in within day and between day assessments, respectively. For the dynamic navicular drift, changes must be in minimum 1.5 mm (28%) for within day assessments to be detectable. Reliability of the time points of dynamic navicular drop and drift was negatively influenced by outliers resulting from participants that did not show a clear minimum navicular height or maximum navicular width around 76% SP (see Additional file 1). Disregarding these outliers, the time points of dynamic navicular drop and drift would have been measurable very consistently. All curve parameters from which the dynamic navicular drop and drift were calculated showed poorer reliability than the dynamic navicular drop and drift themselves (repeatability of about factor two increased). Their repeatability was larger for interday measurements but respectively equal for foot strike and minimum/maximum measures. Hence, whether values at foot strike nor minimum or maximum values were predominating the reliability of the dynamic navicular drop and drift and the reliability of the range of motion measures was superior than single features of navicular height and width during the stance phase.
The study was the first that investigated the reliability of 3D approaches to measure dynamic navicular drop and drift and found that both are measurable with similar and good reliability during intraday assessments. For interday assessments, only the dynamic navicular drop can be recommended.
Reliability is influenced by the examiner (marker placement errors), the laboratory setup [30] (instrumental error) and the participants (step-to-step and intersession errors) which together make up the total error in gait analysis based on multiple trials [31]. Instrumental errors, which were previously optimized [30], and step-to-step errors contributed to the reported trial-to-trial errors (Fig. 2). Intraday reliability was influenced by trial-to-trial and intersession errors, interday reliability additionally by potential marker placement errors. The latter are thought to be of minor importance for the current investigation because: (1) a calibration routine guaranteed the reference planes for measuring navicular height and width being properly aligned and (2) the same experienced examiner placed the markers. Hence, the larger interday reliabilities would have been primary induced by increased intersession errors during interday compared with intraday assessments. This may be explained by learning effects or natural day-to-day variations in movement patterns. For example, if participants were initially not familiar with walking barefoot they could have become more aware about how to walk barefoot from participating in the study that potentially affected movement patterns during the repeated assessments. Further, it has to be acknowledged that the experiments were carried out under laboratory conditions, where for example targeting the force plates might have initially bothered participants from walking naturally. Hence, that participants learned to deal with the laboratory conditions with increasing practice might have been another learning effect that affected movement patterns.


In contrast to recent multi-segment models consisting of many markers, relying on rigid body assumptions and delivering rather complex 3D rotations to characterize foot function, the presented minimal approach uses projections of the trajectory of the navicular tuberosity to the foot’s transverse and sagital planes. This allows to transfer the known concept of navicular drop and drift assessment from quasi-static clinical test conditions to functional tasks, which is recommended to more closely relate assessments to the functional behavior of the foot [8]. We hope that simultaneously measuring medio-lateral in addition to vertical navicular movement will add value to foot function assessment in the future. We propose a concept which allows adding intrinsic foot movement measurements to standard 3D gait analysis protocols. Using a minimal markerset consisting of only four markers provides a “lightweight” approach in terms of preparation time, laboratory infrastructure requirements and post-processing efforts. Validating the proposed measurements of navicular mobility against a multi-segment foot model would be valuable in the future to further substantiate the use of dynamic naviular drop and drift as parameters for foot function. Further, it should be evaluated if the proposed model outputs may discriminate patient groups, reveal intervention effects or identify risk factors for injury and dysfunction.


The authors would like to thank Duncan Cummins for professional proof-reading.


This study was supported by a grant from the Swiss National Science Foundation (SNSF), Project 140928.

Availability of data and materials

The dataset is available online in text files with comma separated values (.csv) bundled in a ZIP-archive provided as Additional file 2.

Ethics approval and consent to participate

The study was approved by the ethics committee of the canton of Bern (KEK number Z007/12) and all participants signed informed consent.

Consent for publication

All participants were informed that the findings of the work may be published based on their data in fully anonymized form.

Competing interests

The funding source had no involvement in study design, data collection, analysis and interpretation of data, in writing the manuscript and in the decision to submit the manuscript for publication. All authors state no conflict of interest.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://​creativecommons.​org/​licenses/​by/​4.​0/​), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://​creativecommons.​org/​publicdomain/​zero/​1.​0/​) applies to the data made available in this article, unless otherwise stated.

Onze productaanbevelingen

BSL Podotherapeut Totaal

Binnen de bundel kunt u gebruik maken van boeken, tijdschriften, e-learnings, web-tv's en uitlegvideo's. BSL Podotherapeut Totaal is overal toegankelijk; via uw PC, tablet of smartphone.

Extra materiaal
Additional file 1 Bland-Altman diagrams for navicular height and width at foot strike and minimum and maximum of navicular height and width, respectively, during stance. Bland-Altman diagrams for the time points of minimum navicular height and maximum navicular width are also presented. (PDF 711 kb)
Additional file 2 Dataset. The dataset used for this article. A ZIP-archive with one text file (comma separated values;.csv) per participant, containing the stance-normalized navicular height and width time series from all analyzed steps. (ZIP 8335 kb)
Additional file 3 Individual model outputs. Navicular height and width during stance from all individuals. (PDF 2324 kb)
Over dit artikel

Andere artikelen Uitgave 1/2018

Journal of Foot and Ankle Research 1/2018 Naar de uitgave