A Novel Sensorised Insole for Sensing Feet Pressure Distributions

Wearable sensors are gaining in popularity because they enable outdoor experimental monitoring. This paper presents a cost-effective sensorised insole based on a mesh of tactile capacitive sensors. Each sensor’s spatial resolution is about 4 taxels/cm2 in order to have an accurate reconstruction of the contact pressure distribution. As a consequence, the insole provides information such as contact forces, moments, and centre of pressure. To retrieve this information, a calibration technique that fuses measurements from a vacuum chamber and shoes equipped with force/torque sensors is proposed. The validation analysis shows that the best performance achieved a root mean square error (RMSE) of about 7 N for the contact forces and 2 Nm for the contact moments when using the force/torque shoe data as ground truth. Thus, the insole may be an alternative to force/torque sensors for certain applications, with a considerably more cost-effective and less invasive hardware.


Introduction
The sense of touch is pivotal in many fields of engineering, such as robotics and wearable technologies [1]. Tactile sensing allows one to characterise physical contacts in terms of contact forces and moments between bodies. In human-centred technologies, instead, tactile sensors are embedded into wearable devices to provide feedback information on human movements. This paper contributes towards human-centred wearable technologies by proposing a novel capacity-based sensorised insole that can measure foot pressure distributions.
A popular application of human-centred wearable devices concerns gait analysis, that is, the systematic study of human walking. This is used in the medical field for the assessment of gait pathologies and for rehabilitation [2][3][4][5]. Gait analysis is also used for sport purposes, that is, to help athletes improve performances while minimising the risk of injuries at the joint level [6][7][8]. Different solutions for gait event detection and plantar pressure monitoring have been proposed, including force platforms, as well as the moments about horizontal axes with an accuracy comparable to commercial FT sensors. In the best case, the performances achieved by our insole, expressed in root mean square error (RMSE), were around 7 N for the contact forces and 2 N m for the contact moments when using force/torque sensors as ground truth. Instead, in most dynamic cases, higher RMSE could be measured, up to 96 N for the forces and 12 N m for the moments.
The paper is organised as follows: Section 2 describes the technology used for building the sensorised skin of the iCub humanoid robot, the mathematical model of taxels and the methods already developed for calibrating capacitive sensors. Section 3 focuses on the insole prototype description and the calibration procedure. Section 4 presents data extracted for the calibration experiments and describes validation experiments and results. This section also describes the tool developed for real-time visualisation of the quantities estimated by the insole. Section 5 contains a summary of the paper and the main advantages of using our insole. It also suggests some future perspectives from the point of view of the hardware, the mathematical model and the calibration technique.

Mathematical Notation
• Let R be the set of real numbers.
• Let x ∈ R n denote an n-dimensional column vector, while x denotes a scalar quantity.
• Let ||x|| be the Euclidean norm of the vector x.
• Let 0 m×n be the zero matrix ∈ R m×n . • Let I n be the identity matrix ∈ R n×n . • Let f z , m x and m y denote the vertical ground reaction force and the moments about the horizontal axes, respectively. • Let x i ∈ R denote a variable referred to the i-th sensor.

Tactile Sensors
The insole proposed in this paper is based on the same technology used for the skin of the iCub humanoid robot [24]. The iCub skin (Figure 1a) is an array of triangular modules of flexible printed circuit boards (PCBs), each containing 12 capacitive sensors with a circular shape (named taxels, Figure 1b): 10 measuring the change in capacitance induced by pressure and 2 embedded in the PCB used for the temperature compensation. Each triangle has an area of 3.24 cm 2 and has a unique 2-bit address. The PCB provides a spatial resolution of approximately 4 taxels/cm 2 (the area of each taxel is approximately 19 mm 2 ). The triangles also host an integrated capacitance to digital converter (CDC) that encodes measurements of each taxel into 16-bit words with a 0.12 fF resolution. A patch of PCB is composed of many interconnected triangles where each connection includes a dedicated Inter-Integrated Circuit (I2C) bus with four synchronised serial data lines and a shared serial clock. This allows up to 16 triangles to be connected to the same MTB. The MTB is a micro-controller-based board that interfaces the I2C bus and the CAN bus ( Figure 1d). The micro-controller recognises the devices connected to the I2C bus (i.e., the CDCs), reads measurements sent by the CDC of each triangle on the bus and performs an on-board temperature and an offset compensation. Then, the MTB cuts the 16-bit words into 8-bit words suitable for communication through CAN bus. The measurements are sent to a host through the CAN bus at a rate of 50 Hz. A deformable dielectric layer, coupled with a ground plane, is placed on top of the sensors (Figure 1c). The dielectric layer consists of a nylon mesh with good mechanical properties in terms of hysteresis and deformation response. The mesh is provided by commercial suppliers and is cost-effective. When it is pressed, the distance between sensors and the surface of the dielectric decreases, causing a change in capacitance. That is, the higher the pressure on the dielectric, the higher the capacitance variation, as shown in Figure 2. The user can set a bit-shift parameter that selects which 8 bits are considered of the 16. The settable value is in the range from 0 (the MTB selects the least significant byte) to 7 (the MTB selects the most significant byte). With a higher bit-shift value, the 8-bit word maps a larger range of analog capacitance values but reduces the sensibility to very small variations. The choice of this value depends on the specific application and on the specific stiffness of the used dielectric material.

PCB Dielectric Ground Plane
Object pressing the dielectric < l a t e x i t s h a 1 _ b a s e 6 4 = " A W 0 l 5 + / x W n + n y P

A Mathematical Model for Taxels
The behaviour of each taxel can be approximated as a parallel plate capacitor described by [30] where C a is the analog capacitance, 0 is the vacuum permittivity, r is the relative permittivity of the dielectric that covers the PCB, A is the area of the taxels (equal for each taxel) and d is the distance between the ground place and the taxels ( Figure 2). The mechanical behaviour of the skin can be approximated as a set of parallel springs connecting the ground plane and the sensors plane. The behaviour can be described with Hooke's law [31] as: where k is the stiffness of the skin and d 0 is the distance between the skin and the taxel, at rest conditions. The stiffness is approximated to a constant, thus, the model (2) is linear. The insole has a smooth, approximately flat shape. We assume the ground contact surface to be smooth as well, such that the insole deformation is locally flat around each taxel. Consequently, we consider only the force perpendicular to the sensors. The force can also be expressed in terms of differential pressure relative to the atmospheric pressure, such that when no external force is applied to the skin (lifted foot) the differential pressure is zero. By assuming a uniform pressure P over the taxel area A, the resulting force F is given by: By substituting (2) and (3) in (1), the capacitance measured by each taxel is given by: By substituting (1) and (2) in (3) and solving for P, we write the pressure estimated for each taxel as: In (5), the pressure is linear with respect to C −1 a but the CDC circuits measure the analog capacitance C a , convert and stream the capacitance as a raw value that we designate as digital capacitance. The relationship between C a and C d is expressed by a linear function: By considering the variation of both capacitances C a and C d , the linear relationship is still valid, as follows: where The variables ∆C d,max and ∆C a,max are the maximum variation of the digital and analog capacitance, respectively. The substitution of (4), (8a) and (8b) in (7) yields the pressure exerted on each taxel as: Nevertheless, the model in (9) can be applied only in ideal conditions. The variables d 0 , k and r are typically obtained via a calibration procedure which is affected by several limitations: r changes with the pressure because the quantity of air contained between the ground plane and the sensors decreases to zero (we can consider that beyond this point we are saturating the sensors); -k changes with the pressure and with time; -d 0 is not constant due to the hysteresis in the dielectric mechanical elasticity.
To overcome some of these drawbacks, a polynomial model can be used to describe the relationship between capacitance and pressure, thus overcoming some of the limitations of the model (9). Hysteresis phenomena, however, will not be taken into account.

Recall on Skin Calibration
The objective of this section is to present two different methods used in previous work for calibrating the iCub skin [24]. Both methods choose a fifth-order polynomial model and use setups based on a similar approach. They exploit an air pump and flexible material to uniformly distribute pressure on the surface of the skin.
• The first setup ( Figure 3a) consists of inserting the skin inside a Ziploc bag and pumping out air.
This decreases the pressure in the bag and creates a uniform pressure distribution over the skin surface. The flow rate is regulated by using valves to ensure a slow and steady pace in order to avoid the dynamic effects of the dielectric material. The pressure reached on the skin is measured with a pressure sensor that sends data to a computer through a CAN network [28]. • The second setup uses an isolation chamber ( Figure 3b). Contrary to the vacuum bag experiment explained above, the pressure is now increased in the external environment relative to the skin. The setup consists of an air compressor, the mentioned chamber, a pressure regulator and a micro-controller for pushing air and increasing the pressure inside the chamber [29].
In both cases, the capacitances measured by each taxel are related to the applied pressure by using a least-square fitting to find the polynomial coefficients.
In this paper, we use the first setup to take measurements then fused with data coming from sensorised shoes (Section 2.3).

Contribution
This section describes our insole prototype by detailing the changes to the capacity-based technology recalled in Section 2.2. It also presents a new setup for the insole calibration as well as the mathematical theory developed for the parameter identification for the taxels' polynomial model.

Insole Prototype
The insole prototype is based on the tactile capacitive technology presented in Section 2.2 ( Figure 4a). The insole is mainly composed of four elements: a support, the electronics, a skin dielectric and the MTBs. The support is an element not present in the iCub skin, as the electronics are directly glued on the iCub cover. In our case, the electronic components are not directly glued on the internal base of the shoe, but on a dedicated support. This choice comes from the fact that the insole is continuously subjected to stress due to foot steps and high pressures, but at the same time the electronic components are very fragile and need to be protected. The support is a 3-mm-thick 3D-printed piece of plastic (Figure 4c) with gaps to accommodate the CDC (Figure 4b). The electronics of each insole are composed of 28 triangles, divided into three patches highlighted in Figure 4b: front in green, middle in red and rear in blue. Each patch is connected to a different MTB. The MTBs are connected in series and they send information to a host computer connected to the insoles through a CAN bus interface. A total of 280 capacitive pressure sensors and 56 capacitive temperature sensors are available on each insole. The electronics are covered by the same cost-effective dielectric material used for the iCub skin, but made thicker to sustain the human weight range. The dielectric layer is covered by conductive lycra and is glued to the PCB (Figure 4a). The total weight of each insole is about 100 g. The cost of a pair of insoles is not higher than 500 euros.

New Calibration Method
As done in previous work [32], we are interested in estimating the pressure distribution, contact forces and moments using only our insole prototype. The relationship between these quantities and the digital capacitance is obtained through calibration. Nevertheless, the calibration methods described in Section 2.4 cannot be used for the calibration of our insole prototype owing to the following drawbacks: • The vacuum bag method can reach a maximum pressure of 1 Atm, which is not high enough for our application. We need higher pressures during the calibration in order to take data that is consistent with real working conditions. For instance, if we consider a person weighing 80 kg s standing still on a single foot heel, a vertical force f z of about 800 N is applied on a small area of the insole. This force is spread over 4~5 triangles. The area of each triangle is about A = 3.24 × 10 −4 m 2 . Thus, recalling (3) and considering four triangles, the pressure on each sensor will be: • The calibration setup based on the isolating vessel is not big enough to hold the proposed insole and can only reach a maximum pressure of 3 Atm.
Because of these drawbacks, we propose a new calibration setup and method. We combine data of two different setups for the insole calibration. The first setup is the vacuum bag. Although it does not allow pressures higher than 1 Atm to be reached, it allows to: (i) have a uniform pressure distribution over the insole; (ii) control the flow rate for increasing the pressure; (iii) measure the applied pressure as described in Section 2.4. These characteristics are important for avoiding the dynamic effects of the dielectric and thus characterising only the static behaviour of the skin. The second setup is meant to overcome the limits of the vacuum bag method. The idea is to compare reliable and accurate force/moment measurements from third-party sensors with the respective quantities estimated through the insoles. Among those sensors there are sensorised treadmills, force plates and force/torque sensors. By placing the insole on one of the mentioned sensors, they can be used simultaneously and data from the insole and from force/moments sensors can be acquired in a synchronised way. This kind of setup allows higher pressures to be reached (depending on the weight of the person) and allows the taxels to be excited with different pressure distributions. The quantities that will be compared are: pressure, vertical ground reaction force and moments about horizontal axes.

Model Parameters Identification from the New Calibration Method
The theoretical model of taxels described in Section 2.3 considers ideal behaviours and precise knowledge of quantities like the skin stiffness k and the dielectric constant r . In practice, these values are difficult to estimate. We resort to a polynomial function for describing taxels' calibration curves for two reasons: -The quantities d 0 and r to estimate in the model (9) are not linear with respect to the pressure and to find them we need to solve a non-linear optimisation problem; -Hypotheses for deriving the theoretical model are based on ideal behaviours, but the physical insole does not always respect the ideally hypotised case. Thus, a more general model helps to consider factors that we did not take into account in (9).
We use a third-degree polynomial, where the degree was chosen on the basis of the sensors' answer. We checked the shape of the taxel calibration curves relating the pressure P(C i ) and digital capacitances C i acquired by increasing the pressure with the vacuum bag setup. Indeed, the calibration curves often have two inflexions due to a nonlinearity at low differential pressure and near saturation of the sensor, and the two inflexions are typical of a cubic function. Thus, the model for the i-th taxel is defined as: The model may be different for each taxel, and the individual coefficients are estimated through the calibration procedure. It is possible to write each taxel model in a compact form: where contains the powers of the capacitance C i for the sensor i, and k i = a 0,i a 1,i a 2,i a 3,i T contains the respective polynomial coefficients to estimate. The resultant force acting on the insole is computed considering the force on each taxel: where f SK ∈ R 3 , and for each taxel i, A i is the taxel area on the ground and P i is the measured pressure. n is the unit vector normal to the ground and pointing upwards, since each taxel measures a force normal to the sensor plate and we consider both the insole and the ground contact to be flat. Another direct consequence is that it is not possible to compute the ground contact horizontal forces. In addition, all the taxels have the same known area denoted by A. For compactness in the formulation, we consider a reference frame W with its origin placed anywhere on the ground plane, the x and y axes parallel to the ground plane, the z axis pointing up. All the dynamic quantities and taxel positions are expressed in that frame. We haven = 0 0 1 T such that we can consider just the z component of the normal measured forces: being f SK z ∈ R. By defining ϕ f ∈ R 4n×1 as: where n is the number of taxels that compose the insole, and k ∈ R 4n×1 as: we can write the resultant force as: Let X ∈ R 4n×4n and Y ∈ R 4n×4n denote the block diagonal matrices containing respectively the x and y coordinates of each taxel expressed in the reference frame W, in the form: . .
We can then write the resultant moments as: The number of samples acquired is greater than the number of unknown coefficients. Thus, it is possible to define four over-constrained systems in the following form: where Φ ∈ R m×4n is the regressor containing the capacitances' powers measured by the taxels, k ∈ R 4n are the coefficients to estimate and b ∈ R m , with m 4n, containing the known quantities measured by the accurate sensors. In general, a problem in the form Φk = b can be solved through the least-square method. Nevertheless, in our case there are multiple variables to be compared, that is, vertical ground reaction force, moments about horizontal axes and pressure. The optimisation problem could be ill-posed for two reasons: -Variables to compare, force, moments and pressure have different order. Indeed, force can reach a value of about 1000 N, moments can reach a value of about 100 N m and pressure assumes values between 0 and 10 6 Pa. -The sample numbers of the two data sets (FT shoes data set and vacuum bag data set) is not the same because they come from different experiments.
Thus, we need to normalise the variables of the optimisation problem (forces, moments and pressures) by considering the order of magnitude and the sample numbers of each data set. The normalisation is obtained by adopting the weighted least-square method: where f FT z ∈ R l f , m FT x ∈ R l f , m FT y ∈ R l f are measured by the FT sensors, P ∈ R l P is measured by the pressure sensor and they correspond to b in (20). l f and l P are the respective numbers of samples. Instead, regressors are defined as below: The problem can also be written in a quadratic programming form: where The problem could be ill-posed, considering the nature of the problem and the fact that with a sensorised shoes setup we do not have control of the acquired data. To overcome this issue, the problem can be regularised with a Tikhonov regularisation [33] as follows: The above quadratic programming problem is solved using an open-source C++ software package qpOASES [34]. The main advantage of this tool is the possibility of adding constraints to our problem if required.

Experiments for Calibration and Validation
The objective of the experiments was (i) to collect data needed for the calibration of the taxels and (ii) to validate the effectiveness of the identification algorithm in (25). An experimental session was carried out at the Istituto Italiano di Tecnologia (IIT), Genova, Italy. One volunteer (mass: 63 kg) was recruited for the experiments to acquire the calibration dataset and a second volunteer (mass: 54 kg) was recruited for validating the estimation algorithm outcome.
For collecting data to use at the calibration phase, two kinds of data sets were taken.
• One data set was collected using the setup described in Section 2.4. The experiments with a vacuum bag ( Figure 5a) reached a pressure of 90 kPa and the pressure was uniformly distributed over the insole. The experiment consisted of three cycles of pressure reduction and increase ( Figure 6). In Figure 6 we can see that the pressure was negative, which can be explained by the choice of the reference frame W that has the z axis pointing up. Indeed, in order to have compression, the force applied on each taxel has to be negative and as consequence of (3) even the pressure has to be negative. We can even observe that when the pressure increased the measured capacitance decreased. This is because the chips do not stream analog capacitances, as explained in Section 2.3, but the analog capacitances are mapped to a digital value at an 8-bit resolution and the relation (6). The mapping was inverse, that is, values from 255 to 240 indicate that the taxel was not pressed and 0 means that the taxel was measuring the highest analog capacitance that is possible to map with 8 bits. • Another data set was taken using a pair of sensorised shoes developed at IIT (Figure 5b). Each shoe was equipped with two six-axis FT sensors by the iCub group at Istituto Italiano di Tecnologia (IIT) [35], able to measure the forces and moments. As mentioned in Section 3.2, the data sets had different features in order to excite the activation of the taxels as much as possible in order to obtain different combinations to cover the entire capacitance range of each sensor (Figure 7).  Several tasks were performed with the sensorised shoes, including standing still on both feet (T1), standing still on one foot (T2), performing slow movements on the feet (T3) (used for the calibration) and walking (T4-T6) (used for the validation). The pressure reached with these experiments was higher than the one applied on the insole using the vacuum bag (up to 1 MPa) and was not uniformly distributed ( Figure 8). The data sets used for calibration were pre-processed in order to filter the acquisition noise. The filter we used in pre-processing phase was the MATLAB implementation of the Savitzky-Golay filter [36]. The filter parameters were the order and the frame length and we set them to 5 and 41, respectively. Data were processed and analysed with MATLAB R . The results described in the next paragraph were obtained from the left insole. The data sets for the calibration phase were filtered in order to consider only slowly varying data. The reason for this is that, compared to the insole sensors, the FT sensors and the pressure sensor are more responsive to rapid force/moment/pressure changes, which leads to a mismatch in the compared data and might affect the calibration results.
Let n f and n P be the number of samples in the data set collected with the sensorised shoe and the number of samples collected with the vacuum, respectively. The weights of the least-square problem were manually tuned and chosen as w f = 1/n f , w m x = 33/n f , w m y = 10/n f and w P i = 10 −6 /n P , with i = 1, ..., n, where n is the number of taxels on a single insole. We obtained the best performances with the chosen weights. The regularisation coefficient λ was instead set to 10 −9 . The highest force and moments reached were | f z | = 650 N, |m x | = 25 N m and |m y | = 90 N m.

Validation
The validation process involved the comparison between the data measured by the FT sensors mounted on the sensorised shoes (in which the insoles were placed) and the data estimated with the calibrated insoles starting from the measured capacitances and using the relations (13) and (19). The validation was performed in MATLAB. The data sets for the validation phase were not filtered to discard capacitances with high variation, as done with the calibration data sets. Table 1 lists the maximum pressures sensed by the taxels in each task. It is worth specifying that the pressure was not uniformly distributed, and thus only some taxels sensed the maximum value. The accuracy of the calibration was evaluated for different tasks through the estimation of the root mean square error (RMSE) resulting from that comparison ( Table 1). The RMSE varied from 7 to 95 N for the vertical force f z , from 0.7 to 1 N m for the moment about the x axis m x and from 2 to 8 N m for the moment about the y axis. The RMSE computed for tasks T1 and T2 was obtained at steady state; indeed, the experiments consisted of staying still on the insoles without changing the foot position. We can see that the RMSE increased for tasks T3 to T6. The obtained results are justified by the introduction of dynamic movements on the insole, where the dynamics was faster for tasks from T4 to T6, which were the tasks that reached the largest error. In addition to the dynamics that involved transition phases, we needed to consider another factor influencing the results. Indeed, the insole was subjected to a different pressure intensity and pressure distribution than the static case. Tasks from T3 to T6 involved a greater vertical force on smaller areas corresponding to the forefoot and the heel. To overcome this problem, one solution could be the introduction of a dynamic sensors model in place of the static one adopted in this study. In conclusion, it is important to highlight that performances related to the steady state and the transition phase were different. For a better interpretation of the results, Table 2 shows the tracking accuracy of the variables f z , m x and m y [37]. As already commented, the accuracy was better for the static tasks. In Figure 9, we can see the estimated force and moments for the Task T3 and Task T4 data sets. The graph indicates that the estimated force f z and moments m x , m y closely followed the actual force and moments applied on the insole.   Furthermore, the calibration goodness could be analysed by comparing the CoP estimated by the FT sensors and the insole. The CoP estimated by the insole was computed as: where x i and y i represent the x and y components of the i-th taxel in the plane of the insole. Instead, the CoP estimated by the FT sensors was computed as: Figure 11 shows the components x and y of the CoP related to Tasks T3 and T4. It can be observed that the values related to the insole were close to those related to the FT sensors (i.e., the ground truth).   Table 3 shows the RMSE for several tasks. The estimation of the CoP varied between 11 and 38 mm for the x component and between 3 and 7 mm for the y component. Instead, Table 4 contains the tracking accuracy. It shows that, contrary to the force and moments tracking, the accuracy was higher for dynamic tasks. We developed a MATLAB-based tool for online visualisation of the estimated CoP, capacitances and pressure distribution, and estimated forces and moments as shown in Figure 12. Unfiltered data were streamed through the YARP protocol with a frequency of 50 Hz and read in MATLAB by opening a YARP port [38]. Data were received and processed in real-time. The cycle for reading, processing and visualising the results had a frequency of 27 Hz. This means that we processed one sample for every two sent, where the sample indicates the set of taxel measurements in a given instant.
On the left, the figure shows the comparison between the centre of pressure estimated from the FT sensor data and that estimated by the insole. In the middle, the graph contains the comparison between vertical contact force and horizontal moments estimated through the FT sensors and the insole, respectively. On the right, the capacitance distribution is given by the taxel measurements and the pressure distribution was computed using the coefficients estimated during the calibration phase. For visualisation purposes, the estimated forces and the pressure distribution are positive. The insole is very sensitive and thus the measured capacitance values were not equal to the rest value even when there was no weight on the foot. This is due to the contact between the foot and the insole, which becomes strong when the shoes are strapped. To avoid considering fake pressures, we used a threshold on the measured capacitance values set to 30. Furthermore, the pressure values were interpolated over the entire insole surface in order to obtain a more accurate estimation of the pressure distribution and the CoP. Additionally, the capacitance values were interpolated to visualise the capacitance distribution. A video of the online visualisation tool is provided to show the results of the validation in real time. On the left-hand side, the comparison between centres of pressure estimated by force/torque (FT) sensors and insole. In the middle, estimation of vertical contact forces f z (N) and moments m x , m y (N m) using the calibrated insole. On the right-hand side, the pressure distribution.

Conclusions
In this paper, we presented a novel insole prototype for the real-time monitoring of plantar pressure distributions. Compared to the existing devices, the measurements are not affected by temperature changes because of an internal compensation. The sensors cover almost the entire surface of the insole, having a high spatial resolution. We defined and identified a model for each single taxel through a proved calibration procedure, which is required for an accurate estimation of the pressure distribution, vertical forces and moments about the horizontal axes. The insole needs to be calibrated only once and provides a calibration model for each individual taxel. The choice of materials allows a cost-effective wearable sensor that is almost comparable in performance with other expensive solutions. We showed that the calibration of the insole gave good results and the estimated variables were quite reliable when compared to those measured and estimated with the FT sensors. We also described an example of a real scenario, demonstrating that it was possible to obtain good results for real-time applications.
Nevertheless, the proposed prototype can still be improved. A wireless module can be used to remove cables. The PCB and the support for the electronics can be more stretchable and soft to provide a more comfortable wearing experience for the user. Additionally, we can improve the calibration with the following additional mathematical considerations: • A weighted iterative least-square method can be used in order to avoid manually tuning the weights; • The formulation for describing the relationship between capacitance and pressure for each taxel could also take into account some hysteresis in the stiffness of the skin dielectric and the electric state of the sensor; • Constraints on the model can be added in the optimisation problem; • The accuracy and the repeatability of the data collection procedure can be improved.
As future work, we aim to redesign the calibration device mentioned in Section 2.4 in order to reach higher pressures and control the pressure variation applied on the insole. Finally, we wish to keep improving the insole design with the desirable hardware characteristics mentioned above.