CN107464241B - Analysis method based on non-diagnosis and non-treatment purpose plantar pressure - Google Patents
Analysis method based on non-diagnosis and non-treatment purpose plantar pressure Download PDFInfo
- Publication number
- CN107464241B CN107464241B CN201710701325.6A CN201710701325A CN107464241B CN 107464241 B CN107464241 B CN 107464241B CN 201710701325 A CN201710701325 A CN 201710701325A CN 107464241 B CN107464241 B CN 107464241B
- Authority
- CN
- China
- Prior art keywords
- foot
- pressure
- data
- partition
- plantar
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 30
- 238000003745 diagnosis Methods 0.000 title claims abstract description 12
- 238000005192 partition Methods 0.000 claims abstract description 41
- 238000012545 processing Methods 0.000 claims abstract description 37
- 230000005484 gravity Effects 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 16
- 238000000638 solvent extraction Methods 0.000 claims description 9
- 230000007797 corrosion Effects 0.000 claims description 6
- 238000005260 corrosion Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 241001227561 Valgus Species 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 2
- 230000001225 therapeutic effect Effects 0.000 claims 2
- 238000009499 grossing Methods 0.000 claims 1
- 210000002683 foot Anatomy 0.000 description 123
- 238000010586 diagram Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000036544 posture Effects 0.000 description 3
- 230000037396 body weight Effects 0.000 description 2
- 210000004744 fore-foot Anatomy 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010205 computational analysis Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 210000003041 ligament Anatomy 0.000 description 1
- 210000003141 lower extremity Anatomy 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
The invention relates to a method for analyzing plantar pressure based on non-diagnosis and non-treatment purposes, which comprises the following steps: obtaining foot pressure data, and carrying out early-stage processing on the foot pressure data to obtain foot sole pressure indexes (maximum foot pressure, average foot pressure, variance foot pressure, standard deviation foot pressure, pressure intensity ratio and the like); further processing the foot pressure round edge in the foot pressure data after the previous processing to obtain a smooth foot pressure image meeting the preset requirement; and (4) carrying out partition processing on the smooth foot pressure image to obtain foot indexes (foot contact area, foot gravity center, arch length, foot length and foot width), and analyzing the standing posture, the eversion index and the arch index of the testee. The analysis method of the plantar pressure can comprehensively analyze the distribution condition of the plantar pressure, realize the further analysis of the plantar pressure data, optimize the analysis result, reduce the analysis error and accurately obtain various plantar pressure indexes.
Description
Technical Field
The invention relates to the field of analysis and processing of measurement data, in particular to an analysis method based on non-diagnosis and non-treatment target plantar pressure.
Background
The human foot is a complex result of layered communication as a mesh by 26 bones, 33 joints and 126 ligaments, muscles and nerves. Its basic function is mainly to support the body weight, buffer and absorb the impact force; generate forward thrust and help to adjust and maintain the balance of the human body. When a human body stands still or walks dynamically, the sole of a foot is subjected to a reaction force of the ground in the vertical direction under the action of the gravity of the human body, and the force is the sole pressure. The foot of the human body is a relatively small part of the whole body, but the pressure on the sole of the foot is enormous with each foot walk, which is approximately over 50% of the body weight. In terms of daily life, on average, each person walks 8000 to 10000 steps on both feet in about 4 hours per day. This means that a person's sole is subjected to a cumulative pressure of several hundred tons per day, which is a surprising figure. When some pathological changes or dysfunctions occur to the foot structure of the human body and the motion state of the human body is changed, the pressure on the sole and the pressure distribution are correspondingly changed. Different plantar pressure distribution characteristics and modes are revealed by analyzing and researching plantar pressure parameters, and the disease foot cause, the disease course derivation and the foot function evaluation can be analyzed by comparing and researching the plantar pressure parameters of a normal foot and a pathological foot; meanwhile, the analysis of the sole pressure also plays an important role in the design of shoes, the application of intelligent lower limb artificial limbs, the recovery of leg disease patients and the like.
The existing precise measuring equipment is generally only used by professional institutions due to the characteristics of more measuring points, complex control, high price and the like, and on the other hand, the existing simple measuring device cannot be used as the basis of professional analysis due to the characteristics of less measuring points, inaccurate positioning and the like. Moreover, because the tested person has different standing postures each time, the obtained pressure indexes have some differences; the positions of the testees standing each time may be different, and the obtained analysis results have some differences; the sole portion that does not contact the pressure plate has no pressure value and therefore a very complete foot shape cannot be obtained.
Therefore, it is necessary to develop an analysis method that can comprehensively reflect the pressure distribution of the sole of a foot, obtain a complete foot shape, and accurately calculate various foot indexes.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a non-diagnosis and non-treatment purpose-based plantar pressure analysis method, which can comprehensively analyze the distribution condition of plantar pressure according to pressure data, further analyze the plantar pressure data, optimize the analysis result and accurately obtain various plantar pressure indexes.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a non-diagnosis and non-treatment purpose-based plantar pressure analysis method, which comprises the following steps:
step (1), foot pressure data are obtained, and the foot pressure data are subjected to early-stage processing to obtain a sole pressure index;
step (2), further processing the foot pressure round edge in the foot pressure data after the preliminary treatment in the step (1) to obtain a smooth foot pressure image meeting the preset requirement;
and (3) carrying out partition processing on the smooth foot pressure image obtained in the step (2) to obtain a foot index, and analyzing the eversion index and the arch index.
In order to further optimize the above technical solution, the technical measures taken by the present invention further include:
further, the plantar pressure indexes comprise a maximum foot pressure value, a mean foot pressure value, a variance foot pressure value, a standard deviation foot pressure value and a pressure intensity ratio; the foot indexes comprise foot contact area, foot gravity center, arch length, foot length and foot width.
Further, the specific implementation process of the step 1) comprises the following steps:
a) acquiring foot pressure data by using a sensor;
b) calculating the average value, variance and standard deviation of the pressure data according to the foot pressure data obtained in the step a); a pressure center of gravity; the mean, variance, standard deviation and center of gravity of the left and right feet;
c) sequentially carrying out interpolation processing and binarization processing on the pressure data processed in the step a) and analyzing to obtain a connected domain of each coordinate point.
Furthermore, the interpolation processing is specifically to interpolate data increased by 2 × 2 times once, interpolate three times in total, and increase the data amount by 64 times, and it can be understood that the interpolation times can be adjusted accordingly here.
Further, the binarization processing is specifically as follows: the pressure value is set to 1, and the no pressure value is set to 0.
Further, the analyzing step may specifically obtain a connected domain of each coordinate point as follows: the length of the connected domain of each coordinate point in the x-axis direction and the y-axis direction is len _ x and len _ y respectively; wherein if one condition of (len _ x + len _ y) <20, len _ y/len _ x <0.05, len _ y/len _ x >20 is satisfied, the pressure value of the coordinate point is set to 0.
Further, the specific implementation process of step 2) includes:
A) carrying out corrosion expansion treatment on the pressure data subjected to the preliminary treatment in the step 1) by adopting ellipse checking;
B) and D), averagely dividing the pressure data subjected to corrosion expansion treatment in the step A) into two parts from the middle according to an x coordinate, respectively corresponding to the positions of the left foot and the right foot, and storing the corresponding pressure data into the data of the left foot and the right foot.
Further, the size of the elliptical kernel in step a) is adjusted according to a predetermined required degree of edge smoothness.
Further, the specific implementation process of the step 3) includes:
I) partitioning the left foot pressure image and the right foot pressure image processed in the step (2) according to a preset partitioning rule;
II) carrying out inverse interpolation processing on the pressure data partitioned in the step I) to restore the pressure data to the size of the acquired data, so as to obtain a partition effect graph;
III) calculating the inward and outward turning index and the arch index according to the partition effect graph in the step II).
Further, the partition rule in step I) is specifically as follows:
trisecting the left and right foot pressure images processed in the step (2) according to the length, dividing the left and right foot pressure images into an upper part, a middle part and a lower part, respectively calculating the gravity centers of the upper and lower parts of the left foot pressure image and the right foot pressure image, and respectively connecting the gravity centers of the upper and lower parts of the left foot pressure image and the right foot pressure image, wherein the connecting lines are respective L6 foot axes;
setting L1-L5 partition lines which are vertical to the axis of the L6 foot, wherein the L1 partition line and the L5 partition line are respectively tangent to the foot type, and the division of the L2 partition line, the L3 partition line and the L4 partition line respectively reflects the front toe part, the front sole part, the arch part and the heel part according to statistical data;
setting L7-L11 partition lines parallel to the axes of the L6 feet, wherein the L7 partition line and the L11 partition line are tangent to the front half part of the foot shape, the L8 partition line bisects the area between the axes of the L6 feet and the L7 partition line, and the L9 partition line and the L10 partition line trisect the area between the axes of the L6 feet and the L11 partition line;
the foot axis line L6 and the segmentation lines L1 to L11 segment the left and right foot shapes of the foot pressure image into 10 regions (region 1, region 2, region 3, region 4, region 5, region 6, region 7, region 8, region 9, and region 10) in total for each foot shape.
It is understood that more or fewer partition lines may be used for the area division based on the analysis purpose, i.e. the above partition rules may be adjusted appropriately.
Further, the inverse interpolation processing in step II) is specifically: the data amount of the inverse interpolation one time is changed into the original data amount (1/4), the data amount is changed into the original data amount 1/64 by the inverse interpolation three times; it is understood that the number of inverse interpolations can be adjusted appropriately, but the number should be equal to the number of interpolations.
Further, the valgus index is calculated by the formula:
the calculation formula of the arch index is as follows:
where S3 is defined as the area of region 3, S4 is defined as the area of region 4, S5 is defined as the area of region 5, S6 is defined as the area of region 6, S7 is defined as the area of region 7, S8 is defined as the area of region 8, S9 is defined as the area of region 9, and S10 is defined as the area of region 10.
The technical scheme of the invention comprises the following specific implementation processes:
firstly, preprocessing data to obtain some basic foot pressure indexes:
1) calculating the average value, variance and standard deviation of the pressure data according to the foot pressure data;
2) calculating the pressure gravity center:
3) dividing the pressure data into two parts from the middle according to x coordinates, respectively corresponding to the positions of the left foot and the right foot, and respectively calculating the average value, the variance, the standard deviation and the gravity center of the left foot and the right foot;
4) the pressure data is interpolated by 2 × 2 times, three times, and 64 times.
5) Carrying out binarization processing on the pressure data: namely, the pressure value is set to be 1, and the no pressure value is set to be 0;
6) and analyzing to obtain a connected domain of each coordinate point, wherein the lengths of the connected domain of each coordinate point in the directions of the x axis and the y axis are len _ x and len _ y respectively. If one of the conditions of (len _ x + len _ y) <20, len _ y/len _ x <0.05, len _ y/len _ x >20 is satisfied, the pressure value at the coordinate position is set to 0.
Secondly, further processing the edges of the foot pressure graph to obtain a smoother foot shape:
1) performing corrosion expansion processing on the pressure data by using an elliptical kernel, aiming at enabling the foot-shaped edge to be smoother, wherein the size of the elliptical kernel can be adjusted according to the requirement of the smoothness degree of the edge;
2) and dividing the processed pressure data into two parts from the middle according to the x coordinate, respectively corresponding to the positions of the left foot and the right foot, and storing the two parts into the data of the left foot and the right foot.
Thirdly, the foot pressure image is subjected to partition processing, and a foot pressure index is calculated:
1) dividing the left foot pressure image and the right foot pressure image into an upper part, a middle part and a lower part according to the length, calculating the gravity centers of the upper part and the lower part, wherein the connecting line of the two points is the axis of the foot shape and represents the direction of the tiptoes of different people when the people stand;
2) the upper and lower intersection points of the axis and the foot pressure image respectively represent the lengths of the foot shapes of the left foot and the right foot;
3) translating the axis to two sides until no intersection point exists between the axis and the front half part foot shape, wherein the distance between the two parallel lines is the foot width;
4) the left and right foot types are divided into 10 areas (as shown in fig. 2) and 11 dividing lines, and the dividing rule is as follows: l6 is the foot axis. L1, L2, L3, L4, L5 are all perpendicular to the axis of the foot, wherein L1 and L5 are tangent to the foot shape, respectively, and the divisions of L2, L3, L4 are derived from statistical data, reflecting the forefoot, ball, arch and heel portions, respectively. L7, L8, L9, L10, L11 are all parallel to the foot axis, L7 and L11 are tangent to the anterior half of the foot shape, L8 bisects the region between L6 and L7, L9 and L10 trisect the region between L6 and L11;
5) carrying out inverse interpolation processing on the pressure data, wherein the data amount of the inverse interpolation once is changed into the original value (1/4), the inverse interpolation is carried out three times, the data amount is changed into the original value 1/64, and the size of the acquired data is recovered; partitioning according to the partitioning rule to obtain a partitioning effect graph (shown in FIG. 3);
6) calculating the valgus index according to the region division: the sum of the areas 3, 4 and 9 minus the sum of the areas 5, 6, 7 and 10 is the ratio of the sum of the areas 3, 4 and 9;
7) calculating the arch index: the ratio of the area of region 8 to the sum of the areas of regions 3, 4, 5, 6, 7, 8, 9, 10.
Compared with the prior art, the invention has the following beneficial effects:
the analysis method analyzes the plantar pressure indexes, including the maximum value of the plantar pressure, the mean value of the plantar pressure, the variance of the plantar pressure, the standard deviation of the plantar pressure and the pressure intensity ratio; calculating foot indexes including foot contact area, foot gravity center, arch length, foot length and foot width; analyzing the standing posture of the testee, including the body gravity center and the axis direction of the feet; analyzing the type of the arch, analyzing the arch index and the turnover degree according to the sole pressure and the foot index data, and judging the type of the arch. The analysis method of the plantar pressure can comprehensively analyze the distribution condition of the plantar pressure according to the pressure data, further analyze the plantar pressure data, optimize the analysis result, reduce the analysis error, accurately obtain various plantar pressure indexes, obtain complete and accurate data of the plantar pressure under the condition of having a few measurement point measurement values, and save the hardware cost of the sensor.
Drawings
Fig. 1 is a block flow diagram of a method for analyzing plantar pressure in an embodiment of the present invention;
FIG. 2 is a region division diagram obtained according to a division rule in the analysis method of plantar pressure in an embodiment of the present invention;
fig. 3 is a diagram of a zoning effect obtained by performing inverse interpolation processing on pressure data in the analysis method of plantar pressure in an embodiment of the present invention.
Detailed Description
The invention provides a non-diagnosis and non-treatment purpose-based plantar pressure analysis method, which comprises the following steps of: obtaining foot pressure data, and carrying out early-stage processing on the foot pressure data to obtain a sole pressure index; further processing the foot pressure round edge in the foot pressure data after the previous processing to obtain a smooth foot pressure image meeting the preset requirement; and carrying out partition processing on the smooth foot pressure image to obtain foot indexes, and analyzing the eversion index and the arch index.
The following description of the embodiments of the present invention will be made with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
Example 1
As shown in fig. 1, the analysis method of plantar pressure in this embodiment includes the steps of preprocessing pressure data, performing erosion and dilation, partitioning foot pressure, calculating indexes, and the like, and specifically includes the following steps:
1) reading the raw pressure data into m _ dataDrawVector;
2) determining the number of rows m _ nRow and the number of columns m _ nCol according to the specification of the sensor;
3) calculating the mean/variance and standard deviation of the pressure data;
4) calculating the pressure gravity center:
x-axis coordinate of center of gravity: summing products of x coordinates corresponding to all pressure data and the pressure values, and dividing the products by the sum of all the pressure values;
y-axis coordinate of center of gravity: summing products of y coordinates corresponding to all pressure data and the pressure values, and dividing the products by the sum of all the pressure values;
5) dividing the pressure data into two parts from the middle according to x coordinates, and respectively corresponding to the positions of the left foot and the right foot;
6) respectively calculating the average value, the variance, the standard deviation and the gravity center of the left foot and the right foot;
7) interpolation processing is carried out on the pressure data, the interpolation is carried out once for increasing the data by 2 × 2 times, the difference value is three times, and the data amount is increased by 64 times;
8) carrying out binarization processing on the pressure data: namely, the pressure value is set to be 1, and the no pressure value is set to be 0;
9) respectively calculating the length len _ x and len _ y of the communicated area of each coordinate point in the x-axis direction and the y-axis direction; if one of the conditions of (len _ x + len _ y) <20, len _ y/len _ x <0.05, len _ y/len _ x >20 is satisfied, the pressure value at the coordinate position is set to 0;
10) performing corrosion expansion processing on the pressure data by using an elliptical kernel, aiming at enabling the foot-shaped edge to be smoother, wherein the size of the elliptical kernel can be adjusted according to the requirement of the smoothness degree of the edge;
11) dividing the processed pressure data into two parts from the middle according to x coordinates, respectively corresponding to the positions of the left foot and the right foot, and storing the two parts into the data of the left foot and the right foot;
12) trisecting the left foot and the right foot according to the length respectively, dividing the trisections into an upper part, a middle part and a lower part, calculating the gravity centers of the upper part and the lower part, and connecting lines of the two parts are the axes of the foot shapes and represent the standing postures of different people;
13) the length of the left foot and the length of the right foot are represented by the two intersection points of the axis and the foot shape respectively;
14) translating the axis to two sides until no intersection point exists between the axis and the front half part foot shape, wherein the distance between the two parallel lines is the foot width;
15) the left and right foot types are divided into 10 areas (as shown in fig. 2) and 11 dividing lines, and the dividing rule is as follows: l6 is the foot axis. L1, L2, L3, L4, L5 are all perpendicular to the axis of the foot, wherein L1 and L5 are tangent to the foot shape, respectively, and the divisions of L2, L3, L4 are derived from statistical data, reflecting the forefoot, ball, arch and heel portions, respectively. L7, L8, L9, L10, L11 are all parallel to the foot axis, L7 and L11 are tangent to the anterior half of the foot shape, L8 bisects the region between L6 and L7, and L9 and L10 trisect the region between L6 and L11.
16) Carrying out inverse interpolation processing on the pressure data, wherein the data amount of the inverse interpolation once is changed into the original value (1/4), the inverse interpolation is carried out three times, the data amount is changed into the original value 1/64, and the data obtained by recovery acquisition is large; the cells are partitioned according to a partitioning rule;
17) calculating the valgus index according to the region division: the sum of the areas 3, 4 and 9 minus the sum of the areas 5, 6, 7 and 10 is the ratio of the sum of the areas 3, 4 and 9;
18) calculating the arch index: the ratio of region 8 to the sum of regions 3, 4, 5, 6, 7, 8, 9, 10.
The analysis method of the present embodiment uses interpolation to obtain a data set in which the amount of pressure information data is increased by 64 times, and obtains a data set containing a large amount of plantar pressure information data by combining the sampling values and the estimated values. On one hand, more interpolation data are provided for basic data, so that the graph fitting of the basic information has scientific samples, and the graph of the basic information is closer to actual data; on the other hand, enough sole data are provided for calculating more sole data, so that the calculation of other indirect information has data basis, not only can the data information provided in the calculation method be calculated, but also other user-defined personalized data information can be calculated, thereby realizing the comprehensive analysis of the distribution condition of sole pressure, realizing the further analysis of sole pressure data, optimizing the analysis result and accurately obtaining various sole pressure indexes.
Those skilled in the art will appreciate that all or part of the computational analysis steps of the methods in the above embodiments may be implemented by associated hardware as an instruction of an associated program, and the program may be stored in a computer-readable storage medium, and the storage medium may include: read-only memory, random access memory, magnetic or optical disk, and the like.
The embodiments of the present invention have been described in detail, but the embodiments are merely examples, and the present invention is not limited to the embodiments described above. Any equivalent modifications or alterations to this practice will occur to those skilled in the art and are intended to be within the scope of this invention. Accordingly, equivalent changes and modifications made without departing from the spirit and scope of the present invention should be covered by the present invention.
Claims (8)
1. A method for analyzing plantar pressure based on non-diagnostic and non-therapeutic purposes, comprising the steps of:
step (1), foot pressure data are obtained, and the foot pressure data are subjected to early-stage processing to obtain a sole pressure index;
step (2), further processing the foot pressure round edge in the foot pressure data after the preliminary treatment in the step (1) to obtain a smooth foot pressure image meeting the preset requirement;
step (3), the smooth foot pressure image obtained in the step (2) is subjected to partition processing to obtain foot indexes, and the eversion index and the arch index are analyzed;
wherein, the specific implementation process of the step (1) comprises the following steps:
a) acquiring foot pressure data by using a sensor;
b) calculating the average value, variance and standard deviation of the pressure data according to the foot pressure data obtained in the step a); a pressure center of gravity; the mean, variance, standard deviation and center of gravity of the left and right feet;
c) sequentially carrying out interpolation processing and binarization processing on the pressure data processed in the step a) and analyzing to obtain a connected domain of each coordinate point;
the specific implementation process of the step (2) comprises the following steps:
A) performing corrosion expansion treatment on the pressure data subjected to the preliminary treatment in the step (1) by ellipse checking;
B) and D), averagely dividing the pressure data subjected to corrosion expansion treatment in the step A) into two parts from the middle according to an x coordinate, respectively corresponding to the positions of the left foot and the right foot, and storing the corresponding pressure data into the data of the left foot and the right foot.
2. The method for analyzing plantar pressure based on non-diagnosis and non-treatment purposes as claimed in claim 1, wherein the interpolation processing is specifically to interpolate data increased by 2 × 2 times once, interpolate three times in total and increase the data amount by 64 times.
3. The analysis method based on non-diagnosis and non-treatment purpose plantar pressure according to claim 1, characterized in that the binarization processing specifically comprises: setting a pressure value to be 1 and setting a no-pressure value to be 0; the connected domain of each coordinate point obtained by the analysis is specifically as follows: the length of the connected domain of each coordinate point in the x-axis direction and the y-axis direction is len _ x and len _ y respectively; wherein if one of conditions (len _ x + len _ y) <20, (len _ y/len _ x) <0.05, (len _ y/len _ x) >20 is satisfied, the pressure value of the coordinate point is set to 0.
4. The method for analyzing plantar pressure based on non-diagnostic and non-therapeutic purposes according to claim 1, wherein the size of the elliptical nucleus in step A) is adjusted according to a predetermined required degree of edge smoothing.
5. The method for analyzing plantar pressure based on non-diagnosis and non-treatment purposes according to claim 1, wherein the step (3) is implemented by the following steps:
I) partitioning the left foot pressure image and the right foot pressure image processed in the step (2) according to a preset partitioning rule;
II) carrying out inverse interpolation processing on the pressure data in the step I) to restore the pressure data to the size of the acquired data, and obtaining a partition effect graph;
III) calculating the inward and outward turning index and the arch index according to the partition effect graph in the step II).
6. The method for analyzing plantar pressure based on non-diagnosis and non-treatment purposes according to claim 5, wherein the partition rule in the step I) is as follows:
trisecting the left and right foot pressure images processed in the step (2) according to the length, dividing the left and right foot pressure images into an upper part, a middle part and a lower part, respectively calculating the gravity centers of the upper and lower parts of the left foot pressure image and the right foot pressure image, and respectively connecting the gravity centers of the upper and lower parts of the left foot pressure image and the right foot pressure image, wherein the connecting lines are respective L6 foot axes;
setting L1-L5 partition lines which are vertical to the axis of the L6 foot, wherein the L1 partition line and the L5 partition line are respectively tangent to the foot type, and the division of the L2 partition line, the L3 partition line and the L4 partition line respectively reflects the front toe part, the front sole part, the arch part and the heel part according to statistical data;
setting L7-L11 partition lines parallel to the axes of the L6 feet, wherein the L7 partition line and the L11 partition line are tangent to the front half part of the foot shape, the L8 partition line bisects the area between the axes of the L6 feet and the L7 partition line, and the L9 partition line and the L10 partition line trisect the area between the axes of the L6 feet and the L11 partition line;
the foot axis line L6 and the segmentation lines L1 to L11 segment the left and right foot shapes of the foot pressure image, respectively, and each foot shape is divided into 10 regions in total.
7. The method for analyzing plantar pressure based on non-diagnosis and non-treatment purposes according to claim 5, wherein the inverse interpolation processing in the step II) is specifically as follows: the data amount is inverted once (1/4) and is inverted three times, so that the data amount is changed to 1/64.
8. The method for analyzing plantar pressure according to claim 6, wherein the plantar pressure is measured by a computer,
the calculation formula of the valgus index is as follows:
the calculation formula of the arch index is as follows:
where S3 is defined as the area of region 3, S4 is defined as the area of region 4, S5 is defined as the area of region 5, S6 is defined as the area of region 6, S7 is defined as the area of region 7, S8 is defined as the area of region 8, S9 is defined as the area of region 9, and S10 is defined as the area of region 10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710701325.6A CN107464241B (en) | 2017-08-16 | 2017-08-16 | Analysis method based on non-diagnosis and non-treatment purpose plantar pressure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710701325.6A CN107464241B (en) | 2017-08-16 | 2017-08-16 | Analysis method based on non-diagnosis and non-treatment purpose plantar pressure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107464241A CN107464241A (en) | 2017-12-12 |
| CN107464241B true CN107464241B (en) | 2020-06-26 |
Family
ID=60549899
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710701325.6A Active CN107464241B (en) | 2017-08-16 | 2017-08-16 | Analysis method based on non-diagnosis and non-treatment purpose plantar pressure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107464241B (en) |
Families Citing this family (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107967931B (en) * | 2017-12-29 | 2024-01-26 | 新绎健康科技有限公司 | Body balance data collector |
| CN108836337A (en) * | 2018-05-04 | 2018-11-20 | 福建省莆田市双驰智能信息技术有限公司 | A method of personalized sufficient type health detection is carried out by foot movement state |
| CN108630316B (en) * | 2018-05-08 | 2021-07-16 | 四川大学 | Fitting method and system for diabetic foot protection shoes |
| CN109141714A (en) * | 2018-08-09 | 2019-01-04 | 重庆邮电大学 | A kind of infant foot bottom pressure DATA REASONING and processing method |
| CN110123330A (en) * | 2019-05-20 | 2019-08-16 | 浙江大学 | Foot type parameter measurement and pressure cloud atlas generation method based on plantar pressure data |
| CN110097551A (en) * | 2019-05-23 | 2019-08-06 | 福州市翼谷科技有限公司 | The operation method of flat foot detection device and its data acquisition based on image analysis |
| CN110037704A (en) * | 2019-05-23 | 2019-07-23 | 福州市翼谷科技有限公司 | Flat foot detection device and its operation method of data acquisition |
| CN110292384A (en) * | 2019-06-26 | 2019-10-01 | 浙江大学 | A kind of intelligent foot arch index measurement method based on plantar pressure data distribution |
| CN110664409B (en) * | 2019-10-24 | 2022-03-25 | 瑞昌芯迈科技有限公司 | Arch type identification method based on pressure acquisition |
| CN111358471B (en) * | 2020-04-15 | 2023-04-28 | 青岛一小步科技有限公司 | Body posture detection device and detection method |
| CN111358469B (en) * | 2020-04-15 | 2024-01-02 | 青岛一小步科技有限公司 | Plantar pressure detection device, foot health detection device and foot health detection method |
| CN111759310A (en) * | 2020-07-16 | 2020-10-13 | 深圳国微云技术有限公司 | Foot detection equipment based on 3D pressure data acquisition |
| CN113100751A (en) * | 2021-04-08 | 2021-07-13 | 宜步医疗科技(江苏)有限公司 | Judgment index and judgment basis of biological force line |
| CN114010181B (en) * | 2021-10-15 | 2024-01-16 | 成都识康堂医疗科技有限公司 | Method and system for converting resistance value into plantar pressure value |
| CN113907745B (en) * | 2021-10-15 | 2024-01-16 | 成都识康堂医疗科技有限公司 | System for converting resistance value into plantar pressure value |
| CN114469067A (en) * | 2022-01-27 | 2022-05-13 | 合肥工业大学 | Foot structure and stress characteristic measuring method and device |
| CN115797499A (en) * | 2022-11-24 | 2023-03-14 | 北京华益精点生物技术有限公司 | Gravity center trajectory drawing method and device, electronic equipment and storage medium |
| CN117100250B (en) * | 2023-07-10 | 2024-05-28 | 北京华益精点生物技术有限公司 | Toe position determining method and device, electronic equipment and storage medium |
| CN117958801B (en) * | 2024-04-02 | 2024-06-11 | 西安福音康复健康产业有限公司 | Calculation method and system for arch coefficient |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140009875A (en) * | 2012-07-13 | 2014-01-23 | 주식회사 학산 | Foot pressure measurement and insol pad selection system and thereof method |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101502419B (en) * | 2009-02-23 | 2013-02-13 | 上海体育学院 | Method and device for measuring arch of foot |
| CN101912271A (en) * | 2010-07-19 | 2010-12-15 | 中国科学院深圳先进技术研究院 | System and method for displaying plantar pressure in real time |
| CN202376277U (en) * | 2011-05-25 | 2012-08-15 | 功泰实业股份有限公司 | Integratedly-formed arch pad improved structure |
| CN102670207B (en) * | 2012-05-15 | 2013-09-04 | 北京大学 | Gait analysis method based on plantar pressure |
| CN104366906B (en) * | 2014-11-11 | 2017-02-22 | 广州励途健康科技有限公司 | arch contour detection method and system |
| CN106096033A (en) * | 2016-06-27 | 2016-11-09 | 北京大学深圳研究生院 | Plantar pressure information analysis method, Apparatus and system |
| CN106805980A (en) * | 2017-01-24 | 2017-06-09 | 重庆小爱科技有限公司 | A kind of gait analysis system and analysis method |
-
2017
- 2017-08-16 CN CN201710701325.6A patent/CN107464241B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20140009875A (en) * | 2012-07-13 | 2014-01-23 | 주식회사 학산 | Foot pressure measurement and insol pad selection system and thereof method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107464241A (en) | 2017-12-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107464241B (en) | Analysis method based on non-diagnosis and non-treatment purpose plantar pressure | |
| CN107123112B (en) | Blood flow state analysis system and method | |
| WO2021000423A1 (en) | Pig weight measurement method and apparatus | |
| Wang et al. | Automatic evaluation of the degree of facial nerve paralysis | |
| CN111358471A (en) | Body posture detection device and detection method | |
| CN111588353A (en) | Body temperature measuring method | |
| CN111383221A (en) | Method for generating scoliosis detection model and computer equipment | |
| CN108159576A (en) | Human body chest and abdomen surface region respiratory movement predicting method in a kind of radiotherapy | |
| Hatamzadeh et al. | A kinematic-geometric model based on ankles’ depth trajectory in frontal plane for gait analysis using a single RGB-D camera | |
| CN105389476B (en) | The interpolation algorithm of IMRT intended dose data based on Gradient Features | |
| Montefiori et al. | Variations of lower-limb joint kinematics associated with the use of different ankle joint models | |
| Bakke et al. | Shape-model scaled gait models can neglect segment markers without consequential change to inverse kinematics results | |
| Jimeno-Morenilla et al. | GNG based foot reconstruction for custom footwear manufacturing | |
| Price et al. | A model‐based motion capture marker location refinement approach using inverse kinematics from dynamic trials | |
| CN115375787A (en) | Artifact correction method, computer device and readable storage medium | |
| Milligan et al. | Within-participant variance in multiplanar breast kinematics during 5 km treadmill running | |
| CN113012812B (en) | Knee, ankle joint and plantar pressure data integration method and system | |
| CN111938655B (en) | Orbit soft tissue form evaluation method, system and equipment based on key point information | |
| KR102565867B1 (en) | Apparatus and method for measuring ground reaction force | |
| CN110755084B (en) | Motion function assessment method and device based on active-passive and staged actions | |
| CN112998695B (en) | Calcaneus angle quantification and pressure treadmill data integration method and system | |
| Hannah et al. | Relationship between sagittal plane kinematics, foot morphology and vertical forces applied to three regions of the foot | |
| Rungruangbaiyok et al. | Shoulder Angle Measurement (SAM) system for home-based rehabilitation using computer vision with a web camera. | |
| CN116959742B (en) | Blood glucose data processing method and system based on spherical coordinate kernel principal component analysis | |
| CN112053330B (en) | Diaphragm prediction system and method based on PCA and TSSM models |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |