JP2015165895A - Gait analysis device and gait analysis program - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a gait analysis device that can estimate a walking state in a short time.SOLUTION: A gait analysis device has a determination unit 42, a feature value calculation unit 460 and an estimation unit 48. The determination unit 42 determines a walking start time when a subject starts walking, based on the acceleration measured by a measuring unit 40 worn by the subject. The feature value calculation unit 460 calculates at least any one of acceleration feature values of: the acceleration in the vertical direction of the subject measured by the measuring unit 40 during a predetermined period from the walking start time as a period when the movement in the walking of the subject does not become routine, when the walking start time is determined; and the acceleration in the horizontal direction substantially perpendicular to the traveling direction of the subject. The estimation unit 48 identifies the walking state of the subject, based on the feature value calculated by the feature value calculation unit 460.

Description

  Embodiments described herein relate generally to a gait analyzer and a gait analysis program.

  2. Description of the Related Art Conventionally, apparatuses for evaluating walking motion using various sensor devices have been developed for the purpose of observing disease progress, preventing falls, and the like. Moreover, in the medical field, the risk of falls is determined by having the subject perform an action (such as standing on one foot) that requires balance ability and observing the behavior.

JP2011-92696A

  However, in the prior art, measurement and estimation of the degree of disease are made on the assumption that the subject is continuously walking, and estimating the walking state and the possibility of falling in a short time from the start of walking is could not. Moreover, in the technique implemented in the medical field, there was a problem that when the subject was allowed to perform an action that required balance ability, the subject himself / herself was at risk and burden.

  The walking analysis device of the embodiment includes a determination unit, a feature amount calculation unit, and an estimation unit. The determination unit determines a walking start time at which the subject starts walking based on the acceleration measured by the measurement unit worn by the subject. When the start time of walking is determined, the feature amount calculation unit determines the vertical direction of the target person measured by the measurement unit during a predetermined period from the start time of walking as the period when the movement of the target person is not steady. A feature quantity of at least one of acceleration and acceleration in a horizontal direction substantially orthogonal to the traveling direction of the subject is calculated. The estimation unit identifies the walking state of the subject based on the feature amount calculated by the feature amount calculation unit.

The hardware block diagram of the gait analyzer concerning an embodiment. The schematic diagram which shows the mounting state of the gait analyzer concerning embodiment. The functional block diagram of the walk analysis device concerning an embodiment. The graph which shows the relationship between an acceleration, a determination result, and an analysis object. The graph which shows the example of analysis of feature-value calculation. The chart which shows the label of three classes learned as an identification example. The conceptual diagram which shows the example of an output of an output part. The flowchart which shows operation | movement of the gait analyzer concerning embodiment.

  Embodiments of a gait analyzer will be described in detail below with reference to the accompanying drawings. FIG. 1 is a configuration diagram illustrating the configuration of a gait analyzer 1 according to an embodiment. As shown in FIG. 1, the gait analyzer 1 includes a main body 10 and a mounting part 30.

  The main body unit 10 includes a control unit 12, a storage unit 14, an input unit 16, an output unit 18, a detection unit 20, and a communication unit 22. The control unit 12 includes, for example, a CPU 120 and controls each unit constituting the main body unit 10. The storage unit 14 includes a ROM (Read Only Memory) and a RAM (Random Access Memory) not shown, and stores a program executed by the control unit 12 and data used when the control unit 12 executes the program. In addition, the main body 10 is provided with a storage medium 140 such as a memory card that can exchange programs and data with the storage unit 14 in a detachable manner.

  The input unit 16 is, for example, an input key or a switch, and accepts a user input to the main body unit 10. The output unit 18 includes, for example, a display unit 180 such as a liquid crystal panel, a speaker 182 that outputs sound, a vibrator (not shown), and the like, and displays the state of the main unit 10 and the results of processing operations of the main unit 10 on a screen. Outputs at least one of sound and vibration. Further, the input unit 16 and the display unit 180 may be integrated by a touch panel or the like.

  The detection unit 20 is, for example, a triaxial acceleration sensor with a sampling frequency of 128 Hz, and the acceleration measurement measurement range is, for example, ± 6 G or more. For example, when the gait analysis device 1 is attached to the subject of the gait analysis, the detection unit 20 is in a vertical direction, a horizontal direction (left and right direction) substantially orthogonal to the subject's traveling direction and the subject's traveling direction. Acceleration (acceleration in three directions) is detected.

  The communication unit 22 is a general-purpose interface that performs communication with the outside, and is configured to be connectable to, for example, wired communication, long-distance wireless communication, and proximity wireless communication.

  The mounting unit 30 is configured by, for example, a belt wound around the waist of the subject, and the body 10 is mounted near the waist of the subject, for example.

  FIG. 2 is a schematic view illustrating a state in which the gait analyzer 1 is mounted near the waist of the subject and the direction of acceleration measured by the gait analyzer 1. As shown in FIG. 2, when the gait analyzer 1 is mounted near the waist of the subject, the gait analyzer 1 is substantially in the vertical direction (Y direction), the direction of travel of the subject (Z direction), and the direction of travel of the subject. The acceleration in the orthogonal horizontal direction (X direction) is detected.

  Next, the function of the gait analyzer 1 will be described. FIG. 3 is a block diagram illustrating an overview of the functions of the gait analyzer 1. FIG. 4 is a graph illustrating the relationship between the acceleration detected by the detection unit 20, the determination result of the determination unit 42 described later, and the analysis target extracted by the extraction unit 44.

  As illustrated in FIG. 3, the gait analyzer 1 includes a measurement unit 40, a determination unit 42, an extraction unit 44, an analysis unit 46, an estimation unit 48, and an output unit 18. The output unit 18 illustrated in FIG. 3 corresponds to the output unit 18 illustrated in FIG.

  The measurement unit 40 detects the movement of the subject. Specifically, the measurement unit 40 includes the detection unit 20, and measures the motion (acceleration) of the subject by detecting, for example, acceleration in three directions that changes according to the motion of the subject. For example, when the power source of the main body unit 10 is turned on, the measurement unit 40 continuously measures acceleration in three directions in order to output a result of the gait analyzer 1 estimating the walking state of the subject. All acceleration data measured by the measurement unit 40 is stored until, for example, the storage unit 14 (FIG. 1) becomes unnecessary with time (or time).

  The determination unit 42 determines whether the subject has started walking based on the result measured by the measurement unit 40. Specifically, the determination unit 42 first acquires the acceleration data stored in the storage unit 14, and determines a predetermined time frame (first time interval a) at a predetermined time interval (first time interval a). The variance value of the acceleration data measured by the measuring unit 40 within one set period A) is calculated. For example, as shown in FIG. 4, the determination unit 42 newly measures the measurement unit 40 every 0.5 seconds (a = 0.5 seconds) and every 0.5 seconds (A = 0.5 seconds). A variance value of acceleration data is calculated. When the sampling frequency of the detection unit 20 is 128 Hz, the acceleration data measured by the measurement unit 40 within the first set period A (A = 0.5 seconds) is 64 per direction. The first time interval a is preferably 0.5 seconds or less. Further, the first setting period A may be set longer than the first time interval a.

  Next, the determination unit 42 (FIG. 3) determines that the variance value of the acceleration data is the threshold value after a period in which the variance value of the acceleration data is equal to or less than a predetermined threshold value (threshold value σ) continues for the second set period B or more. When the period exceeding σ continues for the third set period C or longer, the time when the variance value of the acceleration data starts to exceed the threshold value σ (or finally determined that the variance value of the acceleration data is equal to or less than the threshold value σ) Time) is determined to be the walking start time when the subject started walking. The threshold σ is set to 0.04, for example. The second set period B and the third set period C preferably take values in the range of 2 to 6 seconds.

  For example, as illustrated in FIG. 4, the determination unit 42 determines that the acceleration data variance value after the period in which the acceleration data variance value is equal to or less than the threshold σ continues for 6 seconds (second setting period B = 6 seconds) or longer. When the dispersion value of the acceleration data starts to exceed the threshold σ when the period during which the threshold value exceeds the threshold σ continues for 6 seconds (third set period C = 6 seconds) or longer (start of the third set period C) Time) is determined retroactively to be the time when the subject started walking. In the example illustrated in FIG. 4, the determination unit 42 determines that the measurement unit 40 is 214 seconds and the time 208 seconds is the walking start time. Thus, the determination unit 42 determines whether or not the subject has started walking while determining the walking start time. The determination unit 42 may be configured to calculate a standard deviation instead of the variance value.

  When the determination unit 42 determines that the subject has started walking, the extraction unit 44 (FIG. 3) extracts the analysis target from the result measured by the measurement unit 40. Specifically, when the determination unit 42 determines that the subject has started walking, the extraction unit 44 measures the extraction unit 44 during the extraction period from the start of walking until the fourth set period D elapses. The acceleration data is extracted from the acceleration data stored in the storage unit 14. The fourth setting period D (extraction period) is set to 3 seconds, for example.

  Furthermore, the extraction unit 44 is measured by the measurement unit 40 during a plurality of other extraction periods in which the start time of the extraction period is sequentially delayed from the walking start time at a predetermined time interval (second time interval b). A plurality of acceleration data is also extracted from the acceleration data stored in the storage unit 14.

  For example, as illustrated in FIG. 4, when the determination unit 42 determines that the subject has started walking, the extraction unit 44 passes 3 seconds (fourth set period D = 3 seconds) from the start of walking. Acceleration data measured by the measurement unit 40 is extracted during the extraction period D1 until. Further, the extraction unit 44 measures, for example, six other extraction periods D2 to D7 in which the start timing of the extraction period is sequentially delayed from the start of walking at an interval of 0.5 seconds (second time interval b = 0.5). Six acceleration data measured by the unit 40 are also extracted. Hereinafter, when any of the extraction periods D1 to D7 is indicated without being specified, it may be simply referred to as an extraction period D (or a fourth set period D). The acceleration data extracted by the extraction unit 44 is used as an analysis target by the analysis unit 46.

  The analysis unit 46 (FIG. 3) includes a feature amount calculation unit 460 and a variation value calculation unit 462. The feature amount calculation unit 460 calculates the analysis target feature amount extracted by the extraction unit 44. Specifically, the feature amount calculation unit 460 receives the analysis target extracted by the extraction unit 44, calculates, for example, an autocorrelation function of acceleration data in the vertical direction for each extraction period D as shown in FIG. The first peak value (or the second peak value) is taken as the maximum autocorrelation value. The feature amount calculation unit 460 may be configured to calculate an autocorrelation function after passing the analysis target through a low-pass filter. In FIG. 5, the extraction period D is simply denoted as D, and the extraction periods D1 to D7 are simply denoted as D1 to D7. The feature amount calculation unit 460 outputs each calculated maximum autocorrelation value to the variation value calculation unit 462 and the estimation unit 48 as a feature amount for each extraction period D.

  Note that the autocorrelation function calculated by the feature amount calculation unit 460 has a value other than 0 when the acceleration data has periodicity, and has a characteristic that the value increases as the amplitude of change in the acceleration data increases and the noise decreases. is there. For example, as shown in FIG. 5, in the extraction period D1, the maximum autocorrelation value is smaller than that in the extraction periods D2 to D7. That is, the extraction period D <b> 1 is a transition state between a state in which no motion is performed (such as a stationary state) and a movement becomes steady, and it is considered that the change amplitude of the acceleration data is small and the periodicity is low.

  In addition, the feature amount calculation unit 460 may be configured to calculate the left-right symmetry of the analysis target extracted by the extraction unit 44 as the feature amount with respect to the lateral acceleration detected by the detection unit 20. For example, the feature amount calculation unit 460 may calculate an average value of acceleration data as the feature amount with respect to the acceleration in the left-right direction (X direction) detected by the detection unit 20 in each extraction period D.

  The fluctuation value calculation unit 462 (FIG. 3) calculates a numerical value (fluctuation value) that serves as an index of fluctuation of the maximum autocorrelation value for each extraction period D calculated by the feature quantity calculation unit 460. For example, the fluctuation value calculation unit 462 may calculate a variance value, a standard deviation value of the maximum autocorrelation value (see the variation shown in FIG. 5) for each extraction period D, or the maximum autocorrelation value of the extraction period D as a reference and other values. A value using the sum of differences from the maximum autocorrelation value in the extraction period D is calculated as a variation value. Note that when the subject's movement becomes steady, the variation in the maximum autocorrelation value is considered to be small.

  Further, the fluctuation value calculation unit 462 may be configured to calculate the fluctuation value of the left-right symmetry of the analysis target extracted by the extraction unit 44 with respect to the lateral acceleration detected by the detection unit 20. For example, the fluctuation value calculation unit 462 may calculate the variance value of the average value of acceleration data in the horizontal direction (X direction) calculated by the feature amount calculation unit 460 in each extraction period D as the fluctuation value.

  Here, it is thought that the subject who is excellent in balance ability will be in a steady walking state earlier, so that the subject who is excellent in balance ability is also expected to reduce the fluctuation of the feature amount for each extraction period D early.

  The estimation unit 48 receives the feature amount calculated by the feature amount calculation unit 460 and the variation value calculated by the variation value calculation unit 462, and uses at least one of the received feature amount and variation value to determine the walking state of the subject person. presume.

  For example, the estimation unit 48 receives the maximum autocorrelation value in the extraction period D1 and the variance of the maximum autocorrelation values in the extraction periods D1 to D7 as inputs, and uses the discriminator to determine the walking state of the subject safely and carefully (required). Caution) ・ Classify into 3 groups (3 classes) of danger.

  The estimation unit 48 uses, for example, an SVM (Support Vector Machine) algorithm for the classifier. SVM is a two-group pattern identification method that can deal with nonlinear classification using a kernel function, but the estimation unit 48 uses a one-against-one method, one-against-all method, etc. for the purpose of multi-group identification. Using the extended method using a plurality of classifiers, classification into three groups is realized.

  For example, the estimation unit 48 classifies the walking state of the subject person into the three classes shown in FIG. Each class has a safety, caution and danger label (class label).

  Here, the label will be described in detail. As shown in FIG. 6, for example, each label is associated with a range of points indicating a walking state. The range of the score associated with each label corresponds to the score of BBS (Berg Balance Scale). BBS is a balance test that consists of 14 types of actions such as standing on one foot and changing direction, and the evaluation target's balance ability is distinguished by a score of 0 to 56 points (0 to 4 points for 14 types of actions). To do.

  The estimation unit 48 is trained based on data measured in advance. For the subjects including those who have a relatively high possibility of falling and healthy subjects in advance, the above-mentioned calculation of the feature amount and fluctuation value of acceleration during walking and the provision of points using BBS are prepared as sample data. The correspondence between the feature value and the fluctuation value and the range of the score indicating the walking state is learned.

  That is, the estimation unit 48 learns in advance about the relationship between the three labels for classifying the walking state and the above-described feature amount and variation value, and based on the feature amount and variation value of the subject who newly analyzes the walking state. Then, the walking state of the subject person is identified (estimated).

  Note that the estimation unit 48 may be configured to estimate the walking state of the subject using at least one of the feature amount and the variation value. Moreover, the estimation part 48 may be comprised so that a walking state may be estimated by setting only a threshold value in a feature-value and a fluctuation value, without utilizing a learning algorithm. The estimation unit 48 is configured to execute a pattern recognition algorithm after performing dimensional compression using a neural network as another algorithm or using a self-organizing map, kernel principal component analysis, or the like. Also good.

  The output unit 18 (FIG. 3) accepts and outputs the result estimated by the estimation unit 48. FIG. 7 is a conceptual diagram illustrating an output example of the output unit 18. As illustrated in FIG. 7, the output unit 18 outputs a label or the like indicating the result estimated by the estimation unit 48 by performing display and sound output. The output unit 18 may be configured to output at least one of the result calculated by the analysis unit 46 and the result estimated by the estimation unit 48.

  Next, the operation of the gait analyzer 1 will be described. FIG. 8 is a flowchart illustrating an operation performed by the gait analysis apparatus 1 according to the embodiment by executing a program corresponding to the function illustrated in FIG. 3. The gait analysis apparatus 1 starts gait analysis when the power is turned on, for example, in a state where no exercise such as walking is performed (such as a stationary state).

  As shown in FIG. 8, in step 100 (S100), in the gait analyzer 1, for example, when the power is turned on, the detection unit 20 starts detecting the acceleration of the subject.

  In step 102 (S102), the determination unit 42 determines whether or not the first time interval a has elapsed. When determining that the first time interval a has elapsed (S102: Yes), the determination unit 42 proceeds to the processing of S104, and when determining that the first time interval a has not elapsed (S102: No). ) Continues the process of S102.

  In step 104 (S104), the determination unit 42 calculates the variance value of the acceleration data measured by the measurement unit 40 within the first setting period A, for example.

  In step 106 (S106), the determination unit 42 determines whether or not the variance value of the acceleration data calculated in the process of S104 is equal to or less than a predetermined threshold (threshold σ). If the determination unit 42 determines that the variance value of the acceleration data is equal to or less than the threshold σ (S106: Yes), the process proceeds to S108, and determines that the variance value of the acceleration data is not equal to or less than the threshold σ (S106: In No), the process proceeds to S112. Here, the process may proceed to S102 to continue the processing.

  In step 108 (S108), the determination unit 42 determines whether or not the period in which the variance value of the acceleration data is equal to or less than the threshold σ continues for the second set period B or more. If the determination unit 42 determines that the period in which the variance value of the acceleration data is equal to or less than the threshold σ continues for the second set period B or more (S108: Yes), the determination unit 42 proceeds to the processing of S110, and the variance value of the acceleration data is When it is determined that the period equal to or less than the threshold σ has not continued for the second set period B or longer (S108: No), the process proceeds to S102.

  In step 110 (S110), the determination unit 42 determines whether or not the variance value of the acceleration data is greater than the threshold value σ. If the determination unit 42 determines that the variance value of the acceleration data is greater than the threshold σ (S110: Yes), the determination unit 42 proceeds to the process of S112, and determines that the variance value of the acceleration data is equal to or less than the threshold σ (S110). : No), the process proceeds to S102.

  In step 112 (S112), the determination unit 42 determines whether or not the period in which the variance value of the acceleration data is greater than the threshold σ continues for the third set period C or more. If the determination unit 42 determines that the period in which the variance value of the acceleration data is greater than the threshold σ continues for the third set period C or more (S112: Yes), the determination unit 42 proceeds to the process of S114, and the variance value of the acceleration data is When it is determined that the period longer than the threshold σ has not continued for the third set period C or longer (S112: No), the process proceeds to S102.

  In step 114 (S114), the determination unit 42 determines the walking start time.

  In step 116 (S116), the extraction unit 44 accepts the determination result of the determination unit 42 and extracts the analysis target.

  In step 118 (S118), the feature amount calculation unit 460 receives the analysis target extracted by the extraction unit 44 and calculates the feature amount.

  In step 120 (S120), the variation value calculation unit 462 accepts the feature amount calculated by the feature amount calculation unit 460 and calculates the variation value.

  In step 122 (S122), the estimation unit 48 accepts at least one of the feature amount calculated by the feature amount calculation unit 460 and the variation value calculated by the variation value calculation unit 462, and estimates the walking state of the subject.

  In step 124 (S124), the output unit 18 outputs the result estimated by the estimation unit 48.

  Note that the gait analyzer 1 may be configured to estimate the walking state of the subject using all of the accelerations in the three directions detected by the detection unit 20 (or arbitrarily combined), or in one direction. The walking state of the subject person may be estimated using at least one of the feature amount and the variation value calculated from the acceleration data. In the above-described embodiment, the case where the feature amount calculated by the feature amount calculation unit 460 and the variation value calculated by the variation value calculation unit 462 are distinguished from each other has been described as an example. The variation value calculated by the unit 462 may be configured to estimate the walking state of the subject person as one of the feature amounts. In other words, the walking analysis apparatus 1 may be configured to estimate the walking state of the subject, assuming that any value calculated by the analysis unit 46 is one of the feature amounts corresponding to the motion of the subject.

  When the walking analysis apparatus 1 estimates the walking state of the subject based on the acceleration data of the extraction period D1, the period set in advance as the period during which the subject's movement is not steady is the same as the extraction period D1. (3 seconds in the example shown in FIG. 4). In addition, when estimating the walking state of the subject based on the acceleration data of the extraction periods D1 to D7, a period that is predetermined as a period in which the movement of the subject is not steady is from the start of walking to the end of the extraction period D7. (6 seconds in the example shown in FIG. 4).

  In addition, the gait analyzer 1 is not limited to the configuration shown in the above embodiment, and the main body unit 10 includes, for example, the detection unit 20 (or the measurement unit 40) and the communication unit 22, and is connected via a network. A PC (Personal Computer) or the like may include a determination unit 42, an extraction unit 44, an analysis unit 46, an estimation unit 48, and an output unit 18.

  In addition, the main body 10 may be directly attached to the subject's body by the attachment part 30 made of an adhesive member, or may be attached to an exterior part such as a shoulder bag.

  The program executed by the gait analysis apparatus 1 of the present embodiment has a module configuration including the above-described units (determination unit 42, extraction unit 44, analysis unit 46, and estimation unit 48).

  According to the embodiment described above, even if one sensor is used, the walking state is estimated based on the feature amount at the start of walking of the subject, so that the walking state of the subject can be estimated in a short time. . Moreover, according to the said embodiment, even if it does not perform a complicated act with a high risk and burden for the subject person himself / herself, it is possible to estimate the possibility of falling.

  Moreover, although several embodiment of this invention was described by several combination, these embodiment is shown as an example and is not intending limiting the range of invention. These novel embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the spirit of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, and are included in the invention described in the claims and the equivalents thereof.

DESCRIPTION OF SYMBOLS 1 Walk analyzer 10 Main body part 12 Control part 120 CPU
DESCRIPTION OF SYMBOLS 14 Storage part 140 Storage medium 16 Input part 18 Output part 180 Display part 182 Speaker 20 Detection part 22 Communication part 30 Wearing part 40 Measurement part 42 Judgment part 44 Extraction part 46 Analysis part 460 Feature-value calculation part 462 Fluctuation value calculation part 48 Estimation Part

Claims (11)

Based on the acceleration measured by the measurement unit worn by the subject, a determination unit that determines a walking start time at which the subject started walking; and
When the walking start time is determined, the subject's vertical acceleration measured by the measurement unit during a predetermined period from the walking start time as a period in which the movement of the target person is not steady And a feature amount calculation unit that calculates a feature amount of at least one of accelerations in a horizontal direction substantially orthogonal to the traveling direction of the subject,
An estimation unit for identifying the walking state of the subject based on the feature amount calculated by the feature amount calculation unit;
A gait analyzer. The estimation unit includes
The feature amount calculated by the feature amount calculation unit is input, and the walking state of the target person is identified by a pattern recognition algorithm that outputs any one of a plurality of class labels having different possibility of falls. The gait analyzer according to 1. The determination unit
After a period in which the variance value of the acceleration measured by the measurement unit within the first set period at a first time interval continues for a second set period or longer, the variance value of the acceleration reaches the threshold value. The gait analysis device according to claim 1, wherein when the exceeding period continues for a third set period or more, a time point at which an acceleration variance value starts to exceed the threshold is determined as the walking start time point. The feature amount calculation unit includes:
The gait analysis apparatus according to claim 3, wherein a feature amount of acceleration measured by the measurement unit is calculated from the start of walking until a fourth set period elapses. The feature amount calculation unit includes:
The gait analyzer according to claim 4, wherein at least one of a maximum autocorrelation value in an autocorrelation function of vertical acceleration measured by the measurement unit and an average value of acceleration in the horizontal direction is calculated as a feature amount. The feature amount calculation unit includes:
The gait analysis apparatus according to claim 4, wherein the feature quantity of the acceleration measured by the measurement unit is calculated during a plurality of periods obtained by sequentially delaying the fourth setting period at a second time interval from the start time of walking. The feature amount calculation unit includes:
Further calculating at least one of a plurality of maximum autocorrelation values in the autocorrelation function of vertical acceleration measured in the plurality of periods and a plurality of average values of accelerations in the horizontal direction as a feature amount,
A variation value calculating unit that calculates at least one of a variation value of the plurality of maximum autocorrelation values and a variation value of the plurality of average values;
The estimation unit includes
The gait analysis apparatus according to claim 6, wherein the walking state of the subject is identified based on at least one of a variation value calculated by the variation value calculation unit and a feature amount calculated by the feature amount calculation unit. The estimation unit includes
At least one of the feature amount calculated by the feature amount calculation unit and the variation value calculated by the variation value calculation unit is input, and any one of a plurality of class labels having different fall possibility levels is output. The gait analysis device according to claim 7, wherein the gait state of the subject is identified by a pattern recognition algorithm.   The gait analysis according to claim 8, further comprising: an output unit that outputs at least one of the feature amount calculated by the feature amount calculation unit, the variation value calculated by the variation value calculation unit, and the result identified by the estimation unit. apparatus. Based on the acceleration measured by the measurement unit worn by the subject, a determination unit that determines a walking start time at which the subject started walking; and
When the walking start time is determined, the subject's vertical acceleration measured by the measurement unit during a predetermined period from the walking start time as a period in which the movement of the target person is not steady And a feature amount calculation unit that calculates a feature amount of at least one of accelerations in a horizontal direction substantially orthogonal to the traveling direction of the subject,
Based on the feature amount calculated by the feature amount calculation unit, an estimation unit that classifies the walking state of the subject as safe, caution, and danger;
A gait analyzer. Based on the acceleration measured by the measurement unit worn by the subject, a step of determining a walking start time at which the subject started walking;
When the walking start time is determined, the vertical direction of the subject measured by the measurement unit during a predetermined period from the walking start time as a period during which the movement of the target person is not steady is determined. Calculating a feature quantity of at least one of acceleration and acceleration in a horizontal direction substantially orthogonal to the traveling direction of the subject;
Identifying the walking state of the subject based on the calculated feature amount;
Gait analysis program to run a computer.

Images