Three-dimensional quantitative analysis of healthy foot shape: a proof of concept study

A hybrid statistical morphometry free-form deformation approach to 3D personalized foot-ankle models

Foot and ankle joint models are widely used in the biomechanics community for musculoskeletal and finite element analysis. However, personalizing a foot and ankle joint model is highly time-consuming in terms of medical image collection and data processing. This study aims to develop and evaluate a framework for constructing a comprehensive 3D foot model that integrates statistical shape modeling (SSM) with free-form deformation (FFD) of internal bones. The SSM component is derived from external foot surface scans (skin measurements) of 50 participants, utilizing principal component analysis (PCA) to capture the variance in foot shapes. The derived surface shapes from SSM then guide the FFD process to accurately reconstruct the internal bone structures. The workflow accuracy was established by comparing three model-generated foot models against corresponding skin and bone geometries manually segmented and not part of the original training set. We used the top ten principal components representing 85 % of the population variation to create the model. For prediction validation, the average Dice similarity coefficient, Hausdorff distance error, and root mean square error were 0.92 ± 0.01, 2.2 ± 0.19 mm, and 2.95 ± 0.23 mm for soft tissues, and 0.84 ± 0.03, 1.83 ± 0.1 mm, and 2.36 ± 0.12 mm for bones, respectively. This study presents an efficient approach for 3D personalized foot model reconstruction via SSM generation of the foot surface that informs bone reconstruction based on FFD. The proposed workflow is part of the open-source Musculoskeletal Atlas Project linked to OpenSim and makes it feasible to accurately generate foot models informed by population anatomy, and suitable for rigid body analysis and finite element simulation.

Direct visualization and measurement of the plantar aponeurosis behavior in foot arch deformation via the windlass mechanism

The plantar aponeurosis (PA) is an elastic longitudinal band that contributes to the generation of a propulsive force in the push-off phase during walking and running through the windlass mechanism. However, the dynamic behavior of the PA remains unclear owing to the lack of direct measurement of the strain it generates. Therefore, this study aimed to visualize and quantify the PA behavior during two distinct foot postures: (i) neutral posture and (ii) windlass posture with midtarsal joint plantarflexion and metatarsophalangeal joint dorsiflexion, using computed tomography scans. Six healthy adult males participated in the experiment, and three-dimensional reconstruction of the PA was conducted to calculate its path length, width, thickness, and cross-sectional area. This study successfully visualized and quantified the morphological changes in the PA induced by the windlass mechanism, providing a precise reference for biomechanical modeling. This study also highlighted the interindividual variability in the PA morphology and stretching patterns. Although the windlass posture was not identical to that observed in the push-off phase during walking, the observed PA behavior provides valuable insights into its mechanics and potential implications for foot disorders.

Technical note: A volumetric method for measuring the longitudinal arch of human tracks and feet

Fossil footprints (i.e., tracks) were believed to document arch anatomical evolution, although our recent work has shown that track arches record foot kinematics instead. Analyses of track arches can thereby inform the evolution of human locomotion, although quantifying this 3-D aspect of track morphology is difficult. Here, we present a volumetric method for measuring the arches of 3-D models of human tracks and feet, using both Autodesk Maya and Blender software. The method involves generation of a 3-D object that represents the space beneath the longitudinal arch, and measurement of that arch object's geometry and spatial orientation. We provide relevant tools and guidance for users to apply this technique to their own data. We present three case studies to demonstrate potential applications. These include, (1) measuring the arches of static and dynamic human feet, (2) comparing the arches of human tracks with the arches of the feet that made them, and (3) direct comparisons of human track and foot arch morphology throughout simulated track formation. The volumetric measurement tool proved robust for measuring 3-D models of human tracks and feet, in static and dynamic contexts. This tool enables researchers to quantitatively compare arches of fossil hominin tracks, in order to derive biomechanical interpretations from them, and/or offers a different approach for quantifying foot morphology in living humans.

Human foot form and function: variable and versatile, yet sufficiently related to predict function from form

The human foot is a complex structure that plays an important role in our capacity for upright locomotion. Comparisons of our feet with those of our closest extinct and extant relatives have linked shape features (e.g. the longitudinal and transverse arches, heel size and toe length) to specific mechanical functions. However, foot shape varies widely across the human population, so it remains unclear if and how specific shape variants are related to locomotor mechanics. Here we constructed a statistical shape-function model (SFM) from 100 healthy participants to directly explore the relationship between the shape and function of our feet. We also examined if we could predict the joint motion and moments occurring within a person's foot during locomotion based purely on shape features. The SFM revealed that the longitudinal and transverse arches, relative foot proportions and toe shape along with their associated joint mechanics were most variable. However, each of these only accounted for small proportions of the overall variation in shape, deformation and joint mechanics, most likely owing to the high structural complexity of the foot. Nevertheless, a leave-one-out analysis showed that the SFM can accurately predict joint mechanics of a novel foot, based on its shape and deformation.

Correlation analysis between body mass index and foot length in Chinese adolescents: a regional study

This study aims to analyse the relationship between body mass index and foot length in Chinese adolescents and to provide theoretical guidance for preventing a flat foot in Chinese adolescents. This study recruited 1477 students aged 14-23 years. The participants' height, weight, and body mass index were measured, as well as baseline data, including age, gender and foot length. Differences in foot length (bilateral) and flat foot distribution were statistically significant except for the normal foot and high arch foot distribution based on different body mass index groups. Linear correlation analysis demonstrated that body height, weight and body mass index were positively correlated with bilateral foot length regardless of gender. Body mass index acted as a risk factor for flat foot (bilateral) through disordered multi-classification logistic regression analysis. Body mass index was positively correlated with left and right foot length regardless of gender and acted as a risk factor for a flat foot in Chinese adolescents. Practitioner summary: Significant differences exist in the anthropometric data of various races and ethnic groups. The study was investigated in the form of a cross-sectional study. BMI was positively correlated with bilateral foot length and acted as a risk factor for a flat foot in Chinese adolescents.

Methodological and statistical approaches for the assessment of foot shape using three-dimensional foot scanning: a scoping review

OBJECTIVE

The objectives of this study were to: (i) review and provide a narrative synthesis of three-dimensional (3D) foot surface scanning methodological and statistical analysis protocols, and (ii) develop a set of recommendations for standardising the reporting of 3D foot scanning approaches.

METHODS

A systematic search of the SCOPUS, ProQuest, and Web of Science databases were conducted to identify papers reporting 3D foot scanning protocols and analysis techniques. To be included, studies were required to be published in English, have more than ten participants, and involve the use of static 3D surface scans of the foot. Papers were excluded if they reported two-dimensional footprints only, 3D scans that did not include the medial arch, dynamic scans, or derived foot data from a full body scan.

RESULTS

The search yielded 78 relevant studies from 17 different countries. The available evidence showed a large variation in scanning protocols. The subcategories displaying the most variation included scanner specifications (model, type, accuracy, resolution, capture duration), scanning conditions (markers, weightbearing, number of scans), foot measurements and definitions used, and statistical analysis approaches. A 16-item checklist was developed to improve the consistency of reporting of future 3D scanning studies.

CONCLUSION

3D foot scanning methodological and statistical analysis protocol consistency and reporting has been lacking in the literature to date. Improved reporting of the included subcategories could assist in data pooling and facilitate collaboration between researchers. As a result, larger sample sizes and diversification of population groups could be obtained to vastly improve the quantification of foot shape and inform the development of orthotic and footwear interventions and products.

Foot shape is related to load-induced shape deformations, but neither are good predictors of plantar soft tissue stiffness

Modern human feet are considered unique among primates in their capacity to transmit propulsive forces and re-use elastic energy. Considered central to both these capabilities are their arched configuration and the plantar aponeurosis (PA). However, recent evidence has shown that their interactions are not as simple as proposed by the theoretical and mechanical models that established their significance. Using three-dimensional foot scans and statistical shape and deformation modelling, we show that the shape of the longitudinal and transverse arches varies widely among the healthy adult population, and that the former is subject to load-induced arch flattening, whereas the latter is not. However, longitudinal arch shape and flattening are only one of the various foot shape-deformation relationships. PA stiffness was also found to vary widely. Yet only a small amount of this variability (approx. 10-18%) was explained by variations in foot shape, deformation and their combination. These findings add to the mounting evidence showing that foot mechanics are complex and cannot be accurately represented by simple models. Especially the interactions between longitudinal arch and PA appear to be far less constrained than originally proposed, most likely due to the many degrees of freedom provided by the structural complexity of our feet.

Toward improved understanding of foot shape, foot posture, and foot biomechanics during running: A narrative review

The current narrative review has explored known associations between foot shape, foot posture, and foot conditions during running. The artificial intelligence was found to be a useful metric of foot posture but was less useful in developing and obese individuals. Care should be taken when using the foot posture index to associate pronation with injury risk, and the Achilles tendon and longitudinal arch angles are required to elucidate the risk. The statistical shape modeling (SSM) may derive learnt information from population-based inference and fill in missing data from personalized information. Bone shapes and tissue morphology have been associated with pathology, gender, age, and height and may develop rapid population-specific foot classifiers. Based on this review, future studies are suggested for 1) tracking the internal multi-segmental foot motion and mapping the biplanar 2D motion to 3D shape motion using the SSM; 2) implementing multivariate machine learning or convolutional neural network to address nonlinear correlations in foot mechanics with shape or posture; 3) standardizing wearable data for rapid prediction of instant mechanics, load accumulation, injury risks and adaptation in foot tissue and bones, and correlation with shapes; 4) analyzing dynamic shape and posture via marker-less and real-time techniques under real-life scenarios for precise evaluation of clinical foot conditions and performance-fit footwear development.

New Distinct Component Patterns for Plantar Pressure Variables by Using Principal Component Analysis

BACKGROUND

It is important to determine the plantar pressure distribution of schoolchildren by applying static and dynamic foot analyses using a pedobarography device. However, it is difficult to obtain clear interpretations from results that can be explained by a large number of plantar pressure variables. The aim of this study was to use principal component analysis (PCA) to predict the main components for reducing the size of big data sets, provide a practical overview, and minimize information loss on the subject of plantar pressure assessment in youths.

METHODS

In total, 112 schoolchildren were included in the study (mean ± SD: age, 10.58 ± 1.27 years; body mass index, 18.86 ± 4.33). During the research, a pedobarography device was used to obtain plantar pressure data. Each foot was divided into six anatomical regions and evaluated. Global and regional plantar pressure distributions, load and surface areas, pressure-time integrals, weight ratios, and geometric foot properties were calculated.

RESULTS

The PCA yielded ten principal components that together account for 81.88% of the variation in the data set and represent new and distinct patterns. Thus, 137 variables affecting the subject were reduced to ten components.

CONCLUSIONS

The numerous variables that affect static and dynamic plantar pressure distributions can be reduced to ten components by PCA, making the research results more concise and understandable.

Foot Orthosis and Sensorized House Slipper by 3D Printing

BACKGROUND

In clinical practice, specific customization is needed to address foot pathology, which must be disease and patient-specific. To date, the traditional methods for manufacturing custom functional Foot Orthoses (FO) are based on plaster casting and manual manufacturing, hence orthotic therapy depends entirely on the skills and expertise of individual practitioners. This makes the procedures difficult to standardize and replicate, as well as expensive, time-consuming and material-wasting, as well as difficult to standardize and replicate. 3D printing offers new perspectives in the development of patient-specific orthoses, as it permits addressing all the limitations of currently available technologies, but has been so far scarcely explored for the podiatric field, so many aspects remain unmet, especially for what regards customization, which requires the definition of a protocol that entails all stages from patient scanning to manufacturing.

METHODS

A feasibility study was carried out involving interdisciplinary cooperation between industrial engineers and podiatrists. To that end: (i) For patient-specific data acquisition, 3D scanning of the foot is compared to traditional casting. (ii) a modelling GD workflow is first created to design a process permitting easy creations of customized shapes, enabling the end user (the podiatrist) to interactively customize the orthoses. Then, (iii) a comparison is made between different printing materials, in order to reproduce the same mechanical behavior shown by standard orthoses. To do this, the mechanical properties of standard materials (Polycarbonate sheets), cut and hand-shaped, are compared with four groups of 3D printed samples: poly(ethylene glycol) (PETG), poly(acrylonitrile-butadiene.styrene) (ABS), polycarbonate (PC) and poly(lactic acid) (PLA) obtained by Fused Filament Fabrication (FFF).

RESULTS

Differences found between the foot plaster model obtained with the plaster slipper cast in a neutral position and the model of the real foot obtained with 3D scanning in the same position can be ascribed to the non-stationarity of the patient during the acquisition process, and were limited by a locking system with which no substantial differences in the almost entire sole of the foot scan were observed.

CONCLUSIONS

Using the designed GD workflow, podiatrists with limited CAD skills can easily design and interactively customize foot orthoses to adapt them to the patients' clinical needs. 3D printing enables the complex shape of the orthoses to be reproduced easily and quickly. Compared to Polycarbonate sheets (gold standard), all the printed materials were less deformable and reached lower yield stress for comparable deformation. No modifications in any of the materials as a result of printing process were observed.

Foot-surface-structure analysis using a smartphone-based 3D foot scanner

BACKGROUND

A thorough understanding of the influence of the foot skeletal structure on hallux valgus (HV) is required for HV prevention. We developed a system using a 3D foot scanner on a smartphone to clarify the relationships between foot features and HV risk.

METHODS

Two-dimensional video images were recorded on a smartphone, sent to a computer or cloud server, and used to construct a 3D foot-feature model, considering 10 foot features associated with HV. The participants (419 individuals, aged 40-89 years) stood with their toes 12 cm apart and heels 8 cm apart during video recording. The height and weight were measured for body-mass index calculation.

RESULTS

Age-dependent foot-feature variations were observed slightly for males and distinctively for females. For females, the great toe-first metatarsal head-heel (GFH) angle associated with HV increased with age, i.e., the GFH angle increased with age, suggesting that HV increased with age. Multiple regression analysis revealed that the features determining the GFH angle are the second toe-heel-navicular angle, bone distance axis, and transverse arch length and height. The adjusted coefficients of determination were 0.54 and 0.52 for males and females, respectively.

CONCLUSION

This approach enables simple foot structure assessment for HV risk evaluation.

Reliability and quality of statistical shape and deformation models constructed from optical foot scans

The unique shape of modern human feet, and how they change shape when loaded are thought to be integral to effective upright gait. This unique shape, and the natural variations therein, have previously been analysed using a range of methods; from visual assessments, anthropometric measurements, and footprints, to x-ray, ultrasound and magnetic resonance images. However, these methods are often limited by their use of linear two-dimensional measures. Only recently have advances in three-dimensional (3D) scanning technology and statistical shape analysis been applied to studying 3D foot shape variations. Given their novelty, information regarding the reliability and repeatability of 3D foot scanning and shape modelling is lacking. To investigate whether repeated foot scans captured by two examiners give the same 3D shape and produce consistent statistical shape models, 17 healthy adults' left feet were scanned while bearing half and full bodyweight, as well as minimal weight. Surface to surface distances between corresponding foot meshes and differences between shape model quality criteria were both found to be small and insignificant. The only exception being the specificity criterion for minimally loaded foot scans. Furthermore, Euclidean vectors were used to model the magnitude and direction of deformation that feet undergo as a consequence of increased loading. The deformation models showed that loading a minimally loaded foot results in greater, but less consistent, shape changes than when increasing the load on an already loaded foot. These results show that the utilized methods offer a valuable, reliable and repeatable approach to analysing foot shape and deformation.

Modelling the shape of the pig scapula

BACKGROUND

The shape of pig scapula is complex and is important for sow robustness and health. To better understand the relationship between 3D shape of the scapula and functional traits, it is necessary to build a model that explains most of the morphological variation between animals. This requires point correspondence, i.e. a map that explains which points represent the same piece of tissue among individuals. The objective of this study was to further develop an automated computational pipeline for the segmentation of computed tomography (CT) scans to incorporate 3D modelling of the scapula, and to develop a genetic prediction model for 3D morphology.

RESULTS

The surface voxels of the scapula were identified on 2143 CT-scanned pigs, and point correspondence was established by predicting the coordinates of 1234 semi-landmarks on each animal, using the coherent point drift algorithm. A subsequent principal component analysis showed that the first 10 principal components covered more than 80% of the total variation in 3D shape of the scapula. Using principal component scores as phenotypes in a genetic model, estimates of heritability ranged from 0.4 to 0.8 (with standard errors from 0.07 to 0.08). To validate the entire computational pipeline, a statistical model was trained to predict scapula shape based on marker genotype data. The mean prediction reliability averaged over the whole scapula was equal to 0.18 (standard deviation = 0.05) with a higher reliability in convex than in concave regions.

CONCLUSIONS

Estimates of heritability of the principal components were high and indicated that the computational pipeline that processes CT data to principal component phenotypes was associated with little error. Furthermore, we showed that it is possible to predict the 3D shape of scapula based on marker genotype data. Taken together, these results show that the proposed computational pipeline closes the gap between a point cloud representing the shape of an animal and its underlying genetic components.

Morphological variations of the human talus investigated using three-dimensional geometric morphometrics

INTRODUCTION

The shape of the talus determines the positional and kinematic features of the subtalar, talonavicular, and talocrural joints during walking. Thus, detailed knowledge of the pattern of sexual dimorphism of the human talus may be useful for revealing the pathogenetic mechanism of foot and knee disorders, which are more prevalent in females. The aim of this study was to characterize and visualize the three-dimensional shape variations of the talus in relation to sex and age using geometric morphometrics.

MATERIALS AND METHODS

Computed tomography images of 56 feet without talar injuries or disorders were used in this study. Thirty-seven anatomical landmarks were identified on a bone model of the talus to calculate principal components (PCs) of shape variations among specimens. PC scores were compared between sexes, and their correlations with age were also investigated.

RESULTS

The female talus had a longer neck and narrower head width than the male talus. The superior trochlea was tilted more laterally in the frontal plane in females. Furthermore, the female talar head was more twisted and was more elongated in the dorsoplantar direction.

CONCLUSIONS

Morphological features of the talus in females could alter the subtalar and talonavicular joint kinematics during walking and could be a structural factor in the pathogenetic mechanism underlying foot and knee disorders. This study contributes to the comprehensive understanding of shape variations in the human talus.

Three-dimensional foot shape analysis in children: a pilot analysis using three-dimensional shape descriptors

BACKGROUND

Existing clinical measures to describe foot morphology are limited in that they are commonly two-dimensional, low in resolution and accuracy, and do not accurately represent the multi-planar and complex changes during development across childhood. Using three-dimensional (3D) scanner technology provides the opportunity to understand more about morphological changes throughout childhood with higher resolution and potentially more relevant 3D shape measures. This is important to advance the prevailing arguments about the typical development of children's feet and inform the development of appropriate clinical measures. 3D shape descriptors derived from 3D scanning can be used to quantify changes in shape at each point of the 3D surface. The aim of this study was to determine whether 3D shape descriptors derived from 3D scanning data can identify differences in foot morphology between children of different ages.

METHODS

Fifteen children were recruited from three age groups (2, 5, and 7 years of age). Both feet were scanned in bipedal stance, using the Artec Eva (Artec Group, Luxembourg, Luxembourg) hand-held scanner. Three dimensional shape descriptors were extracted from the 3D scans of the right foot, to create histograms for each age group and heat maps of representative participants for comparison.

RESULTS

There were changes to the dorsal, medial and lateral surfaces of the feet with age. The surfaces became less round along with an increase in indented areas. This is supported by the heat maps which demonstrated that the surfaces of the anatomical landmarks (e.g. the malleoli and navicular tuberosity) became more rounded and protruding, with indented surfaces appearing around these landmarks. On the plantar surface, the concavity of the midfoot was evident and this concavity extended into the midfoot from the medial aspect as age increased.

CONCLUSIONS

The findings of this study indicated that with increasing age the foot becomes thinner in 3D, with bony architecture emerging, and the medial longitudinal arch (MLA) increases in area and concavity. Three-dimensional shape descriptors have shown good potential for locating and quantifying changes in foot structure across childhood. Three-dimensional shape descriptor data will be beneficial for understanding more about foot development and quantifying changes over time.

Analysis of 1.2 million foot scans from North America, Europe and Asia

For decades, footwear brands have developed products using outdated methods and measurements, working with limited insight into the foot shapes and dimensions of their target customers. The integration of 3D scanning technology into footwear retail stores has made it possible for this research to analyze a database containing a large number of male and female 3D foot scans collected across North America, Europe, and Asia. Foot scans were classified into length classes with 5mm length increments; mean width, instep height, and heel width were calculated for each length class. This study confirms the existence of many statistically significant differences in mean foot measurements amongst the regions and between the sexes, and a large dispersion of foot measurements within each group of customers. Therefore, shoes should be developed separately for each group, region, and sex, and at least 3 shoe widths per length class are required to provide a proper fit for 90% of customers. Beyond this, our analysis asserts that a shoe designed for a single group will fit a different segment of the population in another group, and that existing last grading tables should be updated to reflect the foot dimensions of current consumers.

Sex Differences in the Footprint Analysis During the Entire Gait Cycle in a Functional Equinus Condition: Novel Cross Sectional Research

Some studies suggest that gender is related to gait. Females show significantly higher ankle motion and vertical ground reaction forces. Males have significantly larger plantar contact surface areas in all regions of the foot than females in most, but not all, prior studies. However, there is no research on sex differences in a functional equinus condition. In this study, 119 individuals, including 59 females (29.7 ± 5.15 years, 58.74 ± 6.66 kg, 163.65 ± 5.58 cm) and 60 males (31.22 ± 6.06 years, 75.67 ± 9.81 kg, 177.10 ± 6.16 cm), with a functional equinus condition walked onto a pressure platform. In two separate testing sessions, five trials of each foot were conducted for the first, second, and third steps. We measured the contact surface areas for each of the three phases of the stance phase. We computed the intraclass correlation coefficient and standard error of the mean to assess the reliability. We found significantly greater contact surface areas in males than females in the first, second, and third steps in all phases of the stance phase: heel strike, mid-stance, and take-off. This is important information for the design of footwear and orthotics and gender knowledge. In a functional equinus condition, males have registered greater contact surface areas than females in all phases of the dynamic footprint of the stance phase.

A fitting problem: Standardising shoe fit standards to reduce related diabetic foot ulcers

AIMS

Incorrectly fitting shoes are implicated in callus formation and a significant proportion of diabetic foot ulcers, yet remain surprisingly prevalent. We review the current shoe fit guidelines for consistency and discuss ways in which technology may assist us in standardising methods of footwear assessment.

METHODS

Narrative review.

RESULTS

Incorrectly fitted shoes are implicated in the development of some diabetic foot ulcers yet surprisingly there's no consensus on shoe fit, despite substantial spending on prescription footwear. Suggested toe gaps vary from 6 to 20 mm and measurement methods also vary from Brannock Devices and callipers to manual measurement.

CONCLUSIONS

To prevent fit-related foot ulceration, we need to standardise our biomechanical definition of fit. Future research should (1) evaluate the potential use of 3D scanning technology to provide a standardised means of capturing foot morphology; (2) develop a working biomechanical definition of fit, including toe gap through the identification of key physiological markers that capture and predict dynamic foot shape changes during different physical activities and body weight loading conditions; and (3) determine whether changes in dynamic foot shape of those with diabetes differs from those without, impacting on their shoe fitting needs, potentially necessitating specialist footwear at an earlier stage to avoid ulceration.

An assessment of the information lost when applying data reduction techniques to dynamic plantar pressure measurements

Data reduction techniques are commonly applied to dynamic plantar pressure measurements, often prior to the measurement's analysis. In performing these data reductions, information is discarded from the measurement before it can be evaluated, leading to unkonwn consequences. In this study, we aim to provide the first assessment of what impact data reduction techniques have on plantar pressure measurements. Specifically, we quantify the extent to which information of any kind is discarded when performing common data reductions. Plantar pressure measurements were collected from 33 healthy controls, 8 Hallux Valgus patients, and 10 Metatarsalgia patients. Eleven common data reductions were then applied to the measurements, and the resulting datasets were compared to the original measurement in three ways. First, information theory was used to estimate the information content present in the original and reduced datasets. Second, principal component analysis was used to estimate the number of intrinsic dimensions present. Finally, a permutational multivariate ANOVA was performed to evaluate the significance of group differences between the healthy controls, Hallux Valgus, and Metatarsalgia groups. The evaluated data reductions showed a minimum of 99.1% loss in information content and losses of dimensionality between 20.8% and 83.3%. Significant group differences were also lost after each of the 11 data reductions (α=0.05), but these results may differ for other patient groups (especially those with highly-deformed footprints) or other region of interest definitions. Nevertheless, the existence of these results suggest that the diagnostic content of dynamic plantar pressure measurements is yet to be fully exploited.

Effect of In-Shoe Foot Orthosis Contours on Heel Pain Due to Calcaneal Spurs

The objective of this study is to investigate the effect of contouring the shoe insole on calcaneal pressure and heel pain in calcaneal spur patients. Calcaneal pressure was measured using three force sensors from 13 patients including three males and 10 females. These patients have plantar heel pain due to calcaneal spurs, and we examined five customized contour insole foot areas (0–100%). Sensors were attached at the central heel (CH), lateral heel (LH) and medial heel (MH) of the foot. The pain was measured using an algometer and evaluated by the pain minimum compressive pressure (PMCP). In this study, it was observed that the calcaneal pressure decreased with increasing insole foot area. In addition, increasing the insole foot area from 25% to 50% can reduce the calcaneal pressure approximately 17.4% at the LH and 30.9% at the MH, which are smaller than the PMCP, while at the MH, pressure reduced 6.9%, which is greater than the PMCP. Therefore, to reduce pain, one can use 50% insole foot area, even though at MH it is still 19.3% greater than the PMCP. Excellent pain relief was observed when using 100% insole foot area, as the pressures in those three areas are lower than the PMCPs, but it is not recommended because it requires large production costs.