Foot modeling and smart plantar pressure reconstruction from three sensors.

In order to monitor pressure under feet, this study presents a biomechanical model of the human foot. The main elements of the foot that induce the plantar pressure distribution are described. Then the link between the forces applied at the ankle and the distribution of the plantar pressure is established. Assumptions are made by defining the concepts of a 3D internal foot shape, which can be extracted from the plantar pressure measurements, and a uniform elastic medium, which describes the soft tissues behaviour. In a second part, we show that just 3 discrete pressure sensors per foot are enough to generate real time plantar pressure cartographies in the standing position or during walking. Finally, the generated cartographies are compared with pressure cartographies issued from the F-SCAN system. The results show 0.01 daN (2% of full scale) average error, in the standing position.

INTRODUCTION

Foot modeling and plantar pressure measurements have widely focused interest of researchers. Plantar pressure measurements are used in the diagnostics of foot related problems such as plantar foot pain or diabetic foot [1-3]. In the quiet standing position or during walking, the plantar pressures highlight the compression of the soft tissues (skin, fat, ligaments, and muscles) under the feet. The plantar pressure distribution can be measured by using insoles made of a pressure sensor matrix. F-SCAN® mobile System [from Teckscan, Inc. Boston, MA. 02127-1309, USA] with 954 sensors or Pedar® with 256 sensors [from Novel, Gmbh novel GmbH, Ismaninger Str. 51, D-81675 Munich] are the most popular instrumented insole systems used by podiatrists. Information processed from the plantar pressure distributions are used to design a special footwear or insole. This footwear (or foot insole) inserted between the foot and the ground may help in the treatment and prevention of foot pathology [4-6]. The design of the appropriate footwear by the podiatrist is generally based on pressure measurements. This kind of design is rather empirical and needs simulation tools to predict the plantar pressure cartography as a function of the thickness, stiffness and shape of the footwear.

Foot modeling may target different applications: prediction of plantar pressure, behaviour during balance, prevention of ulcer formation, gait or running analyses, etc. Many researchers focus on the biomechanical foot model. One can mention 3D finite element models of the foot and ankle, which have been developed to study ulcer formation [7, 8], or to predict the forces between the foot and different types of support, such as moulded insoles, insoles of varying stiffness [9-12]. Models of coupled foot-boot [13] have also been proposed.

This modeling technique requires prior knowledge of the foot shape in 3D, provided by Magnetic Resonance Imaging (MRI) measurements, Computed Tomography (CT) scanners or video reconstruction [14]. Unfortunately, most podiatrists do not have these expensive and sophisticated devices at their disposal. With the aim of offering the podiatrist a practical tool able to simulate plantar pressure distribution, we have derived a new model of the foot, taking into account the internal foot shape and position, which is suitable for dynamic balance description. The model must also take into account the properties of the shoe insole. Of course this model will have to work effectively with only the data that podiatrists can measure on their own.

The objective is to provide a model generating the plantar pressure cartography for different foot positions (seat) and different elastic media under the foot. At first, the model of a bare foot on planar ground is presented. The extraction of the foot parameters from the plantar pressure cartographies is detailed in the standing position. The results show the extrapolation of the plantar pressure cartography when using three pressure sensors. In the second part, the case of a foot in a shoe is presented and the results during walking are shown.

METHODS

Foot Model and Hypothesis

In the standing position, equilibrium is controlled by the muscle action, which moves the distribution of the plantar pressure by changing the rotation of the foot around the ankle (anterior-posterior) and the repartition of the total weight of the body on both feet (lateral).

The force applied to the ankle compresses the foot, which can be considered foremost as a skeleton covered by an elastic medium made of soft tissues. The compression naturally increases plantar pressure.

One can consider that the foot skeleton can be divided into three parts, from back to front: the heel, the metatarsus and the phalanges. The soft tissues of the plant constitute a complex medium. The thickness and elasticity of the soft tissues depend on their location under the foot, but the heel and the metatarsal heads support the biggest pressure.

A simplified model has been developed in order to describe the relationship between posture and plantar pressure through mechanical foot characteristics.

In standing position, at rest, the sum of the forces perceived on the surface of the two feet corresponds to the force resulting from the total body mass under the acceleration of gravity. When walking, the force resulting from the vertical acceleration of the centre of mass has to be added.

Bare Foot on Planar Ground

At first, one considers the case of a bare foot on a planar ground.

In order to describe the relationship between the force applied to the ankle and the plantar pressure distribution, we can define three elements for the mathematical model: FS, the foot seat, IS the internal foot shape, and EM the elastic medium. The parameters of each element are foot dependent.

Foot Seat

The origin of the internal coordinates of the foot is the ankle, x being the lateral coordinate axis, y the longitudinal one (heel to forefoot) and z the vertical one.

The location of the foot in space is given by the vertical position of the ankle and the rotations around the ankle, as shown in Fig. (). The equation of the “foot seat” FS plan is:

The terms respectively describe the lateral and longitudinal rotations around the ankle, whereas c is the vertical displacement of the ankle.

The considered rotation angles are small. In the approximation of small angles and

Internal Foot Shape

The Internal foot Shape, IS, is a surface profile describing the non-compressive material mainly due to the skeleton.

The internal foot shape is described by the distance of each foot point to the foot seat plane. Since the rotation angles are small, the distance has its main component, , along the z axis. The components along the and x axes due to the rotations are neglected.

IS can be precisely obtained by “classical” 3D measurement techniques, such as MRI or scanning. But it can also be approached by simple plantar pressure measurements.

For a bare foot on flat ground, the foot shape is similar to a 3D profile of the underfoot, such as the one provided by a 3D scanner. The internal foot shape consists in the foot skeleton covered by the other rigid media, such as the ligaments and the muscles, as shown in Figs. ( and ).

Elastic Medium

The internal foot shape (skeleton and rigid medium) is covered by soft tissues, as shown in Fig. (). The thickness and the elasticity of the soft tissues depends on the underfoot location. The compression of the tissues induces a pressure.

To simplify the modeling, it is assumed that the soft tissues can be described realistically by an equivalent uniform elastic medium (EM). The elastic medium is uniform in terms of thickness and elastic properties, Fig. (). For a bare foot, the Elastic Medium, EM, describes directly the soft tissues, which cover the internal foot shape.

The relation between pressure and compression of the EM should verify the following properties:

The EM does not apply any force when the foot is not in contact with the floor.

The EM cannot be compressed beyond its thickness.

The intensity of the force depends on the material stiffness.

For the same compression, the force depends on thickness and stiffness of the EM.

The vertical compression of the EM is directly related to the FS and IS functions. If there is a contact, that is to say if , then .

The usual linear spring compression models, issued from the small variation hypothesis, cannot be applied in this case. A simple function that can account for the above limitations is an exponential function, which links the EM compression to the applied force depending on the characteristics of the EM. It is given by:

or the inverse function for or where (mm) is the compression of the EM, thick (mm) the thickness, stiff (N / cm²) the stiffness, and F (N / cm²) the pressure.

This type of law can be verified through the experimental compression curve of an 8 mm thick elastomer foam with a 25 cm² area is shown on Fig. (). The dots are the measured points, the line is the fitted exponential with a stiffness of 0.2 daN/25cm². This type of material is used for insoles.

The 25 cm² area compares with the minimum support surface of the foot, when standing on the forefoot or the heel.

Pressure Map

The problem is spatially discretised. The size of a pixel corresponds to the F-SCAN one, , which is being used for the experiments.

The pressure cartography, , is directly derived from the foot seat , for a given internal foot shape and a given elastic medium. With, stiff in daN per pixel, the force per pixel is expressed by: with

Direct simulation of the plantar pressure distribution as a function of the foot seat FS can be achieved with arbitrary IS and EM, but for medical applications, the simulation requires the actual internal foot shape of the patient and the elastic medium characteristics.

Our purpose is to identify FS and EM from experimental data.

Internal Foot Shape Extraction

The pressure measurements are made with F-SCAN, at a frequency of 50 scans per second.

From these measurements one can obtain a set of pressure maps for different foot seats.

At each time step, , for each pixel of the foot the force is measured and stored. Note that each time step corresponds to a different foot seat.

and the compression are linked by the bijective relation:

The compression depends on the seat plan and the internal foot shape :

The relation can be written:

The unknowns are the thickness, thick, the stiffness, stiff, the 3 foot seat parameters for each step i and the internal foot shape for each pixel k. But the relation shows that the thickness thick appears as a proportional factor for the other unknown spatial parameters . There is no other relation, so that the thickness cannot be determined independently of . It is fixed to an arbitrary realistic value acting as a proportional factor on the foot seat and the internal foot shape.

For the foot seat parameters are unknown, along with the profile for .

The number of data is equal to . The number of unknown parameters are the N values of , and the values of the foot seat parameters .

The total number of data , is larger than the total unknown values .

After eliminating the unnecessary pixels,N is of the order of 500, and M is of the order of 1000 or more, depending on the duration of the measurement.

Most of the unknown parameters can be obtained by a linear system, which minimizes the estimation error.

Let be E the summation of the square of the errors made on the compression, which can be written:

where the are the experimental data.

The thickness is fixed. The values of the other unknowns can be computed by expressing that the partial derivatives of with respect to each unknown are null at the minimum of E (least squares method).

The partial derivatives, , lead to a linear equation system in . But, do not lead to a linear form. A first order development could have been used. We preferred to compute the stiffness which minimises E with method close to a dichotomy method. In the present case, the variation domain of the unknown is divided in 3 parts instead of 2. At first, the domain of variation of the stiffness is divided in 3 intervals, which correspond to 4 values of the stiffness. For each value of the stiffness, the least squares method is used to determine the parameters. The E value is then computed. The two neighbour values surrounding the stiffness, for which the E value is the lowest, define the new boundaries of the variation domain, which is in turn divided in 3 intervals. This iterative process is repeated until the required accuracy of the stiffness is reached.

The numerical method has been developed with MATLAB software on an ordinary desktop personal computer. The computing time is about 10 seconds.

In the following, we present a measurement of the plantar pressure distribution, on a 70 kg male subject with a foot length of 26 cm, bare-footed, in a standing position. The F-SCAN pressure maps are recorded at a frequency of 50 frames/second for 1 minute.

The thickness of the elastic medium has been fixed at 5 mm, a realistic value that corresponds to the mean footpad thickness under the heel and the forefoot. The calculated internal foot shape is shown in Fig. (). The height varies between 1 and 4.5 mm. A variation of 2 mm can be observed between the mid-foot and the forefoot or the heel.

Fig. () presents the standard deviation for each pixel k, around the average value computed from the different postures i. The amplitude is less than 0.15 mm for the main support zones. The standard deviation is the result of implicit errors due to the simplified model hypothesis and measurement uncertainty.

The low values of standard deviation, compared to the height of Internal Shape, validate the hypothesis and the model. Higher variation amplitude points are located on the toe phalanges, which are articulated with the metatarsus and may move independently to the foot seat.

Other high variation points may appear in the zones where only a few pressure data are available. In any case the precision is less than a few percent.

SUMMARY

After identification of the parameters for each pixel k at each step i the seat and the compression of the elastic medium are known together with the internal foot shape, which is independent of the step.

Then it is possible to calculate , the force applied to the pixel k for any seat defined by

with

RESULTS

Plantar Pressure Extrapolation from 3 Sensors

Since the foot seat is defined by the equation of a plane, three sensors are sufficient to measure it.

The concept has first been tested using an arbitrary IS and fixed EM, varying the foot seat in order to obtain the necessary database. A study has been carried out to evaluate the performances of the method independently of experimental measurement incertitude.

It allowed us to verify that the use of 4 or 5 sensors rather than 3 did not significantly change the precision of the results.

The 3 sensors are fixed on an insole. They are arranged under the heel, the first and fifth metatarsal respectively, as shown in Fig. (). The F-SCAN insole is placed above the sensor insole, in order to compare the plantar pressures and the sensor measurement.

Sensor Parameter Extraction

The effective stiffness value for the sensors has to be determined. The stiffness depends on the surface of the sensor. Presently, the actual surface and response of the sensors are the stiffness value.

The internal foot shape and the seat have been determined for each step . The value of the stiffness of the sensors is calculated using the same method as EM stiffness with the same thickness and from the same data base.

Moreover the exact positions of the sensors with regard to the F-SCAN insole are unknown. Practically the two insoles have no reason to correspond perfectly.

The relative distances of the 3 sensors from each other, on the insole, are known. The shift (rotation and translation) of the sensor coordinate, with regard to the F-SCAN one, is numerically determined simultaneously along with the extraction of the stiffness.

Foot Seat Extraction from 3 Sensors

At each time step i, the force measured by the sensors provides the compression of the elastic material above the sensors. Once the foot shape above each sensor is known, the foot seat extraction can be resumed by solving the equations to obtain the 3 unknowns .

After extracting the FS, and knowing the IS, the plantar pressure distribution map can be calculated.

Plantar Pressure Simulation

For each posture, the foot seat is determined from the 3 sensors responses. Using the foot seat, the calculated foot shape and the stiffness of the elastic medium the pressure maps are then real time simulated.

We now compare the simulated maps to the measured F-Scan pressure maps. Fig. () shows the plantar pressure measured for one of the standing positions that the extrapolated plantar pressure from the model highlights the anatomically support zones under the foot, such as the heel and the forefoot.

One can notice that the peak pressure value under the heel for the left and right foot is about 0.45 daN/(pixel area) for both measured and extrapolated pressures. We see that the pressure under the toes does not appear in this posture.

The average error measured at each pixel is 0.01 daN, with a standard deviation of 0.005 daN. This model therefore offers good prediction characteristics. It will be verified on a wider sample of subjects in further studies.

Foot in a Shoe

Previously we were considering a bare foot. The elastic medium EM was the plantar soft tissues.

If the foot lies on an insole, the EM is composed of several layers of different thicknesses and stiffnesses. Usually, the stiffness of soft tissues and insoles may be comparable. When identifying the stiffness of the EM from pressures data it is impossible to separate the contributions of the soft tissues and the insole.

In this case, we show in Fig. () that a single exponential satisfies the description of a multilayer system for the domain of pressure between 0.15 and 1 daN/pixel, which corresponds to the pressure of greater numerical weight in the linear regression method. It is impossible to determine the zero value of the functions, which implies that the compression can be defined with an additive constant. For example, Fig. () shows the summation, in circular forms, of two 5 mm layers with respective stiffness of 0.1 daN/pixel and 0.4 daN/pixel. The continuous black curve is the fitting exponential of the sum for a thickness of 7.6 mm and a stiffness of 0.25 daN/pixel.

So the elastic material is considered as a mean equivalent medium describing the multilayer behaviour in the main domain of pressures. The model and the method are still valid.

When the foot is inside a shoe, the shape of the shoe insole changes the problem. In such a case the compression is the distance between the internal foot shape and the shoe insole.

The shape of the shoe insole need not be known. The shoe shape will be included in the internal foot shape when identifying the parameters.

Where the new internal foot shape function is

The internal foot shape has to be considered as a numerical function, which may look different from a bare foot shape. But it does not change the capability of the method to extrapolate the foot pressures from 3 pressure sensors.

For illustration purposes, Fig. () presents a comparison between measured and simulated plantar pressure maps during walking for the same subject.

The simulated map for the left foot describes the plantar pressure evolution during heel strike, mid stance, forefoot.

DISCUSSION

The extrapolations of the plantar pressure maps, from 3 sensors per foot, are made by using the foot model developed in this study. This method depends on the foot parameters. For a bare foot on a planar ground, the IS and the stiffness of the EM are determined from the plantar pressure maps issued from F-Scan.

The extraction of the IS from the plantar pressure maps while the subject is standing, exhibits a variance about a few percent under the heel and the mid foot. A more important standard deviation (12 % of the measured height) can be seen under the toes. This variation is due to the relatively low pressure measured under the toes during calibration. A possible solution would be the application of a specific pressure under the toes when making measurements.

Starting from these barefoot measurements, the height of the extracted foot shape is extracted for a selected thickness of 5 mm. This value is comparable to the thickness of the plantar soft tissue. Otherwise, when the patient is standing on an elastomeric insole during measurements, the total thickness value has to be adjusted, depending on the stiffness of each medium and mechanical combination law.

For a foot in shoe, it should be noted that using a compensated insole during the measure of the plantar pressure changes the IS. In this case, the calculated IS include the form of the compensated insoles and the IS calculated when bare foot.

Thanks to the foot model developed in this study, the plantar pressure maps are calculated from the measures of 3 pressure sensors per foot. In the standing position, the average difference on the force map per pixel, compared to F-Scan one is 0.01 daN/Pixel, for a pressure variation between 0 and 0.45 daN/Pixel.

We also verified the validity of the method during walking. The plantar pressures maps have been also extrapolated for the same subject. Results show a high similarity in the pressure maps compared to F-Scan one. Therefore, the plantar pressure extrapolated during walking permits studies on the gait.

Moreover, the results confirm the good approximation of a uniform elastic medium with exponential elasticity properties. This assumption limits the complexity of modeling the plantar soft tissues, and foot insoles.

For simplicity, the model sets that the IS remains rigid when changing the FS. But as the toe phalanges are articulated to the metatarsus, they may move differently from the FS. Some error may occur on the under toe pressures, but they remain small with regards to heel or forefoot pressures.

The results confirm the hypotheses suggested, simplifying the foot model.

Note that the model validation has been carried out on a healthy foot, with “standard” pressure areas. Specific sensor number, location and size may be necessary in other cases. The model aims to provide simplified practical tool to generate pressure maps, inducing a cost/benefits balance concerning sensor number and complexity.

CONCLUSION

We have proposed a new foot model based on simplified anatomical modeling. It provides a practical tool to simulate the plantar pressure distribution for podiatrists.

Based on the foot model, the plantar pressure distribution can be generated using only three pressure sensors per foot. For this, the internal foot shape is first calculated from plantar pressure cartographies, issued from usual matrix systems, such as F-SCAN or Pedar, avoiding the use of expensive equipment.

The sensors are to be used in a portable ambulatory monitoring system, with low power consumption and wireless transmission (Wi-Fi, Bluetooth…). The reduced number of data would be transmitted to podiatrist, or locally processed within a smart phone processor. Plantar pressures cartographies simulations are real time proceeded.

This system could be suitable for the balance study for the elderly or preventing ulcers formation in diabetic foot. These applications will be presented elsewhere.

In the future, this model could be further developed to extrapolate the action of added insoles on the under foot pressure when assigning thickness and stiffness to each area of the elastic medium. It could be an insole design tool for diabetic patients.