WO2014181602A1 - Walking posture meter and program - Google Patents
Walking posture meter and program Download PDFInfo
- Publication number
- WO2014181602A1 WO2014181602A1 PCT/JP2014/058910 JP2014058910W WO2014181602A1 WO 2014181602 A1 WO2014181602 A1 WO 2014181602A1 JP 2014058910 W JP2014058910 W JP 2014058910W WO 2014181602 A1 WO2014181602 A1 WO 2014181602A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- point
- value
- reference period
- temporal change
- walking
- 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.)
- Ceased
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6813—Specially adapted to be attached to a specific body part
- A61B5/6823—Trunk, e.g., chest, back, abdomen, hip
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
- A61B5/112—Gait analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
- A61B5/1121—Determining geometric values, e.g. centre of rotation or angular range of movement
- A61B5/1122—Determining geometric values, e.g. centre of rotation or angular range of movement of movement trajectories
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
Definitions
- the present invention relates to a walking posture meter, and more particularly to a walking posture meter that quantitatively evaluates whether or not a human walking posture is a correct posture.
- the present invention also relates to a program for causing a computer to execute a method for quantitatively evaluating whether or not a human walking posture is a correct posture.
- Patent Document 1 Japanese Patent Application Laid-Open No. 2011-078728. This device detects the direction of gravitational acceleration based on the output from an acceleration sensor attached to the human hip when the human takes a predetermined posture, so that the tilt angle of the hip with respect to the ground when the posture is taken To estimate the pelvic tilt.
- Patent Document 2 Japanese Patent Application Laid-Open No. 2011-2510173 calculates a movement amount from an output of an acceleration sensor attached to a human waist and acquires a walking locus based on the movement amount. To do.
- a walking posture meter that can easily evaluate whether or not a human walking posture is tilted forward or backward with respect to the traveling direction.
- a walking posture meter is a walking posture meter that evaluates the walking posture of the measurement subject, and an acceleration sensor that is mounted on the midline of the measurement subject's waist; Using the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output from the acceleration sensor, an amount corresponding to the forward / backward inclination of the posture of the person being measured is calculated. And an arithmetic unit.
- the “inclination in the front-rear direction” of the posture while walking means the inclination along the traveling direction of the upper body of the human being walking (the upper body from the pelvis).
- the inclination may be for an arbitrary posture (leg phase) that appears while the gait goes around.
- the inclination may be about the inclination of the posture of the person to be measured in the vicinity of the time when the heel of the front leg contacts the road surface.
- the inclination may be about the inclination of the posture of the measurement subject in the vicinity of the time point where the rear leg, which is a free leg, coincides with the front leg, which is a standing leg, along the traveling direction.
- the acceleration sensor is mounted on the midline of the waist.
- the arithmetic unit uses the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output from the acceleration sensor to respond to the degree of inclination of the posture of the subject during walking in the longitudinal direction.
- the amount to be calculated is calculated. Therefore, according to the amount, it is possible to evaluate the degree of inclination in the front-rear direction of the posture of the measurement subject during walking.
- this walking posture meter can perform the above-described evaluation based on the output of the acceleration sensor mounted on the midline of the subject's waist, so that it does not depend on large-scale equipment such as motion capture. It is possible to evaluate easily.
- the calculation unit derives a characteristic parameter indicating characteristics of the temporal change of the vertical axis acceleration and the temporal change of the longitudinal axis acceleration, and the inclination in the front-rear direction based on the characteristic parameter. An amount corresponding to is calculated.
- the calculation unit derives a feature parameter indicating the temporal change characteristics of the vertical axis acceleration and the longitudinal axis acceleration, and calculates an amount corresponding to the tilt in the front-rear direction based on the feature parameter. An amount corresponding to the degree of inclination of the posture in the front-rear direction can be quickly calculated.
- the calculation unit derives a plurality of types of feature parameters, and an amount corresponding to the inclination in the front-rear direction based on an amount obtained by weighted addition of the derived types of feature parameters. Is calculated.
- the calculation unit can evaluate the degree of inclination in the front-rear direction of the posture of the measurement subject while walking by weight-adding the plurality of derived characteristic parameters. Therefore, the calculation required for the evaluation is extremely simple, and quick evaluation is possible.
- the calculation unit estimates the step length of the subject based on a temporal change waveform of the vertical axis acceleration, and performs the weighted addition based on the estimated step length.
- the weighting coefficient to be used is changed.
- the calculation unit estimates the measured person's stride based on the temporal change waveform of the vertical axis acceleration, and changes the weighting coefficient used for the weighted addition based on the estimated stride. Appropriate evaluation is possible regardless of the characteristics of the direction (the width of the stride).
- the calculation unit is configured to calculate the walking cycle of the subject based on at least one of the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration.
- a reference period corresponding to one step in the walking cycle is obtained, and a weighting coefficient used for the weighted addition is changed based on the obtained walking cycle or reference period.
- the calculation unit obtains a reference period corresponding to one step of the walking cycle of the measured person or the walking cycle based on at least one of the temporal change waveform of the vertical axis acceleration and the longitudinal axis acceleration. Since the weighting coefficient used for the weighted addition is changed based on the walking cycle or the reference period, appropriate evaluation can be performed regardless of the characteristics of the person to be measured (walking cycle length).
- the amount corresponding to the inclination in the front-rear direction is an angle formed between the pelvis of the measurement subject and the horizontal direction during walking.
- the calculation unit calculates the tilt angle (angle formed with the horizontal) of the person being measured while walking, and evaluates the tilt of the posture of the person being measured in the front-rear direction using the calculated tilt angle. Therefore, this walking posture meter can perform evaluation with extremely high accuracy.
- a program according to another aspect of the present invention is a program for causing a computer to execute a method for evaluating a walking posture of a measured person, the method including an acceleration mounted on a midline of the waist of the measured person. Using the step of obtaining the sensor output and the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output by the acceleration sensor, Calculating an amount corresponding to the inclination of the direction.
- the computer By executing the program, the computer first acquires the output of the acceleration sensor mounted on the midline of the subject's waist. Then, using the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output by the acceleration sensor, it corresponds to the degree of the forward / backward inclination of the posture of the measurement subject during walking. Calculate the amount. Therefore, the inclination of the posture of the measurement subject during walking can be evaluated according to the amount.
- this program enables the above-mentioned evaluation based on the output of the acceleration sensor mounted on the midline of the subject's waist, so it is easy to use regardless of large equipment such as motion capture. It is possible to evaluate.
- the walking posture meter according to one aspect of the present invention, it is possible to easily evaluate the inclination in the front-rear direction of the posture of the measurement subject during walking.
- FIG. 4A is a diagram showing a mode in which the activity meter is attached to the measurement subject.
- FIG. 4B is a diagram for explaining the X-axis (front-rear axis), the Y-axis (left-right axis), and the Z-axis (up-down axis).
- 5A is a diagram for explaining the relationship between the forward inclination (backward inclination) of the walking posture and the angle between the pelvis and the horizontal.
- 5 (B), 5 (C), and 5 (D) are diagrams showing a posture of a human being walking (timing when a heel of a front leg is grounded).
- FIG. 5B is a schematic diagram when a human walking backward is viewed from the side.
- FIG. 5C is a schematic view of a human walking in a good posture that is neither backward nor forward tilted when viewed from the side.
- FIG. 5D is a schematic diagram when a human walking forward is viewed from the side.
- FIG. It is a figure (vertical axis acceleration time change waveform figure) explaining the characteristic parameter which shows the characteristic of the time change of the vertical axis acceleration which the said acceleration sensor outputs. It is a figure (vertical axis speed time change waveform figure) explaining the characteristic parameter which shows the characteristic of the time change of the vertical axis acceleration which the said acceleration sensor outputs. It is a figure (vertical axis orbit time change waveform figure) explaining the characteristic parameter which shows the characteristic of the time change of the vertical axis acceleration which the said acceleration sensor outputs.
- FIG. 1 shows a system configuration of a walking posture meter (the whole is denoted by reference numeral 1) according to an embodiment of the present invention.
- This walking posture meter 1 includes an activity meter 100 and a smartphone 200.
- the active mass meter 100 and the smartphone 200 can communicate with each other by BLE (Bluetooth low energy; specified in Bluetooth® Core Specification Ver. 4.0).
- BLE Bluetooth low energy
- the activity meter 100 includes a casing 100M, a control unit 110, an oscillation unit 111, an acceleration sensor 112, a memory 120, an operation unit 130, and a display mounted on the casing 100M.
- Unit 140 BLE communication unit 180, power supply unit 190, and reset unit 199.
- the casing 100M is formed in a size that fits in the palm of a human hand so that the activity meter 100 can be easily carried.
- the oscillating unit 111 includes a crystal resonator and generates a clock signal that serves as a reference for the operation timing of the activity meter 100.
- the oscillation unit 111 may be a module chip having a function as a clock generator.
- the acceleration sensor 112 detects the three-axis (three directions) acceleration received by the casing 100M and outputs it to the control unit 110.
- the acceleration sensor 112 may be a module chip of a triaxial acceleration sensor.
- the memory 120 includes a ROM (Read Only Memory) and a RAM (Random Access Memory).
- the ROM stores program data for controlling the activity meter 100.
- the RAM stores setting data for setting various functions of the activity meter 100, acceleration measurement results, calculation result data, and the like.
- the control unit 110 includes a CPU (Central Processing Unit) that operates based on the clock signal.
- the control unit 110 receives a signal from the acceleration sensor 112 according to a program for controlling the activity meter 100 stored in the memory 120. Based on the detection signal, each part of the activity meter 100 (including the memory 120, the display part 140, and the BLE communication part 180) is controlled.
- the control unit 110 includes at least a signal processing system that can process time-series data of vertical axis acceleration and longitudinal axis acceleration.
- the operation unit 130 includes a button switch in this example, and accepts appropriate operation inputs such as a power on / off switching operation and a display content switching operation.
- the display unit 140 includes a display screen made up of an LCD (liquid crystal display element) or an organic EL (electroluminescence) display, and displays predetermined information on the display screen according to a signal received from the control unit 110.
- the display unit 140 may be an LED (light emitting diode) that displays on / off of a power source, an operation state, and the like by turning on, turning off, flashing, and the like.
- the power supply unit 190 is composed of a button battery in this example, and supplies power to each part of the activity meter 100.
- the BLE communication unit 180 communicates with the smartphone 200 in real time. For example, information representing a measurement result is transmitted to the smartphone 200. In addition, an operation instruction is received from the smartphone 200.
- the BLE communication unit 180 may be a module chip having a BLE function.
- the reset unit 199 includes a switch, and resets and initializes the operation of the control unit 110 and the stored content of the memory 120.
- the smartphone 200 includes a main body 200M, a control unit 210, a memory 220, an operation unit 230, a display unit 240, a BLE communication unit 280, and network communication that are mounted on the main body 200M. Part 290.
- This smartphone 200 is obtained by installing application software (computer program) so that a commercially available smartphone can instruct the activity meter 100.
- the control unit 210 includes a CPU and its auxiliary circuit, controls each unit of the smartphone 200, and executes processing according to programs and data stored in the memory 220. That is, data input from the operation unit 230 and the communication units 280 and 290 is processed, and the processed data is stored in the memory 220, displayed on the display unit 240, or output from the communication units 280 and 290. Or
- the memory 220 includes a RAM used as a work area necessary for executing the program by the control unit 210 and a ROM for storing a basic program to be executed by the control unit 210. Further, as a storage medium of an auxiliary storage device for assisting the storage area of the memory 220, a semiconductor memory (memory card, SSD (Solid State Drive)) or the like may be used.
- a semiconductor memory memory card, SSD (Solid State Drive)
- the operation unit 230 includes a touch panel provided on the display unit 240 in this example.
- a keyboard or other hardware operation device may be included.
- the display unit 240 includes a display screen (for example, an LCD or an organic EL display).
- the display unit 240 is controlled by the control unit 210 to display a predetermined image on the display screen.
- BLE communication unit 280 communicates with activity meter 100 in real time. For example, an operation instruction is transmitted to the activity meter 100. In addition, information representing a measurement result is received from the activity meter 100.
- the network communication unit 290 transmits information from the control unit 210 to other devices via the network 900, and receives information transmitted from other devices via the network 900 and passes the information to the control unit 210. be able to.
- this walking posture meter 1 when this walking posture meter 1 is used by a person 90 to be measured as a user, for example, the activity meter 100 is attached by a wearing clip 100C (shown in FIG. 1). It is attached to the back side of the waist on 90 midline 91.
- the longitudinal direction is the X axis
- the horizontal direction is the Y axis
- the vertical direction is the Z axis.
- the acceleration sensor 112 of the activity meter 100 includes an X-axis (front-rear axis) acceleration, a Y-axis (left-right axis) acceleration, a Z-axis (which is measured by the casing 100M as the subject 90 walks forward.
- the acceleration (vertical axis) is output.
- the person under measurement 90 turns on the power of the activity meter 100 and the smartphone 200.
- the application software of the smartphone 200 is activated, and the activity meter 100 is instructed to start measurement via the operation unit 230 and the BLE communication unit 280.
- the person to be measured 90 walks straight forward by a predetermined number of steps, in this example, 10 steps. Then, the person under measurement 90 instructs the activity meter 100 to output the calculation and the calculation result via the operation unit 230 and the BLE communication unit 280 of the smartphone 200.
- control part 110 of the active mass meter 100 works as a calculating part, and performs the calculation mentioned later. Then, information representing the calculation result is transmitted to the smartphone 200 via the BLE communication unit 180.
- FIG. 12 shows an operation flow by the control unit 110 of the activity meter 100.
- the control unit 110 of the activity meter 100 waits for an instruction to start measurement from the smartphone 200 as shown in step S1.
- the control unit 110 acquires the output of the triaxial acceleration by the acceleration sensor 112 as shown in step S2. Acquisition of the output of the acceleration sensor 112 is performed only for a predetermined period (for example, 14 seconds) as a period including acceleration time-series data for 10 steps in this example.
- the acquired acceleration time series data is temporarily stored in the memory 120.
- step S ⁇ b> 3 the control unit 110 waits for a measurement start instruction from the smartphone 200.
- step S4 the control unit 110 calculates an amount corresponding to the inclination angle of the pelvis. And as shown to step S5, the control part 110 works as an evaluation part, and uses the calculation result (pelvic inclination angle estimation result) to measure the degree of inclination in the front-rear direction of the posture of the person being measured while walking. To evaluate. Then, as shown in step S6, the evaluation result is output (transmitted) to the smartphone 200. Note that the control unit 110 may execute the process of step S4 as soon as acceleration time-series data for at least one step is obtained. In that case, the step of determination shown as step S3 may be omitted.
- FIG. 5 (A) is a diagram for explaining the definition of the pelvic inclination angle.
- This figure is a schematic view of the human waist when viewed from the side, and the pelvis PV is represented by a triangle.
- the pelvic inclination angle ⁇ is a line from the superior anterior iliac spine ASIS (ASIS: Anterior Superior Iliac Spine) to the superior posterior iliac spine PSIS (PSIS: Posterior Superior Iliac Spine) with the clockwise rotation direction being positive. The angle between the minute and the horizontal H.
- ASIS Anterior Superior Iliac Spine
- PSIS Posterior Superior Iliac Spine
- a pelvic inclination angle ⁇ of less than +3 degrees is associated with the walking posture “backward inclination”
- a pelvic inclination angle ⁇ of +3 degrees or more and less than +12 degrees is associated with the walking posture “middle”, and +12 degrees.
- the above pelvic tilt angle ⁇ is associated with the walking posture “forward tilt”.
- the boundary value is an example, and other boundary values may be used.
- the posture is divided into three stages of “forward leaning”, “middle”, and “backward leaning”, but other stages may be used.
- FIGS. 5B to 5D are schematic diagrams showing examples of human walking postures “backward tilt”, “middle”, and “forward tilt”, respectively.
- a human who has a tendency to walk in the walking posture “backward tilting” walks in a posture that turns upside down.
- FIG. 5C is an example of the walking posture “middle”. In this way, a human walking in the walking posture “middle” has a straight back and is able to walk correctly (beautifully).
- FIG. 5D a human who has a tendency to walk in a walking posture “forward leaning” walks in a leaning posture with forward leaning. (In addition, some people who are “forward leaning” walk with their stomachs protruding.)
- the vertical axis acceleration changes from negative to positive through the zero cross point.
- the timing at which the stance (right foot in the figure) and the free leg (left leg in the figure) substantially coincide with each other in the traveling direction corresponds to the vicinity of the timing at which the third peak P3 appears.
- a waveform as illustrated and described appears during one step of human walking.
- a period (Step T) from the timing when the forefoot heel contacts the ground (heel contact timing) to the next heel contact timing is defined as a reference period.
- the period from the left foot heel contact timing to the right foot heel contact timing is referred to as the left foot reference period, and from the right foot heel contact timing to the left foot heel contact timing only when there is a particular need to distinguish in the following description.
- the right foot reference period the reference periods of one step with the left foot and one step with the right foot are distinguished.
- the control unit 110 described above uses the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration over at least one reference period in the process of step S4 in FIG. 12 (pelvic inclination angle estimation process). Then, an amount corresponding to the inclination in the front-rear direction of the posture of the measurement subject while walking is calculated.
- the amount corresponding to the degree of inclination in the front-rear direction is not limited to the average pelvic inclination angle ⁇ .
- step S4 in FIG. 12 will be described in detail with reference to the graphs in FIGS. 7 to 11, the flowchart in FIG. 13, and the reference table in FIG.
- Each step of the flowchart of FIG. 13 may be executed by the control unit 110 described above.
- the execution subject of each step may be the control unit 210 instead of the control unit 110.
- step S41 the control unit 110 calculates the vertical axis (Z-axis) acceleration time series data (FIG. 7) and the longitudinal axis (FIG. 7) from the time change waveform of the triaxial acceleration acquired in step S2.
- Control unit 110 similarly detects the next zero-crossing point that changes from negative to positive, and specifies the timing as the start of the next reference period (the period of the current reference period). Thereby, the reference period (Step T) is determined.
- step S42 the controller 110 calculates vertical axis velocity time series data (FIG. 8) and vertical axis trajectory time series data (FIG. 9) from the vertical axis acceleration time series data (FIG. 7).
- step S43 the control unit 110 calculates longitudinal axis velocity time series data (FIG. 11) from longitudinal axis acceleration time series data (FIG. 10).
- step S44 the control unit 110 derives a characteristic parameter indicating characteristics of the temporal change in the vertical axis acceleration and the temporal change in the longitudinal axis acceleration.
- the control unit 110 uses the time-series data of the vertical axis acceleration, the time-series data of the vertical axis velocity, and the time-series data of the vertical axis trajectory to indicate the characteristic of the temporal change of the vertical axis acceleration. Deriving parameters.
- the control unit 110 derives a characteristic parameter indicating a characteristic of the temporal change in the longitudinal axis acceleration using at least one of the time series data of the longitudinal axis acceleration and the time series data of the longitudinal axis velocity.
- the target to be derived in step S44 is a feature amount that represents the temporal change feature of each axis acceleration.
- the control unit 110 derives an amount (feature parameter) representing a characteristic of the temporal change in the acceleration of each axis using the time series data of the velocity of each axis and the time series data of the trajectory.
- characteristic parameters PRM1 to PRM33 there are 33 types of characteristic parameters (characteristic parameters PRM1 to PRM33, that is, PRMi (i: an integer of 1 to 33)).
- PRMi an integer of 1 to 33
- each feature parameter will be described in detail. Note that the number of types of feature parameters and the definition of each feature parameter are examples, and different types of features and definitions of feature parameters may be used.
- the characteristic parameter PRM1 is a parameter related to the appearance timing of the maximum point that appears first when the vertical peak acceleration is based on the first peak in the reference period, that is, the zero cross point that changes from negative to positive.
- PRM1 is a dimensionless quantity T ZAP1 / StepT obtained by normalizing the timing T ZAP1 with a reference period StepT.
- the characteristic parameter PRM2 is a parameter related to the appearance timing of the second maximum peak that appears when the vertical peak acceleration is based on the second peak in the reference period, that is, the zero cross point that changes from negative to positive.
- PRM2 is a dimensionless quantity T ZAP2 / StepT obtained by normalizing the timing T ZAP2 with a reference period Step T.
- the characteristic parameter PRM3 is a parameter relating to the appearance timing of the third maximum point that appears when the vertical peak acceleration is based on the third peak in the reference period, that is, the zero cross point that changes from negative to positive.
- PRM3 is a dimensionless amount T ZAP3 / StepT obtained by normalizing the timing T ZAP3 with a reference period StepT.
- the characteristic parameter PRM4 exists between the second maximum point that appears second and the maximum point that appears third when the vertical valley acceleration is based on the second valley in the reference period, that is, the zero cross point that changes from negative to positive. This is a parameter related to the appearance timing of the minimum point.
- PRM4 is a dimensionless amount T ZAV2 / StepT obtained by normalizing the timing T ZAV2 with the reference period Step T.
- the characteristic parameter PRM5 is a parameter related to the size of the maximum point that appears first when the vertical peak acceleration is based on the first peak in the reference period, that is, the zero cross point that changes from negative to positive.
- PRM5 is a dimensionless quantity ZAP1 / ZAP2 obtained by normalizing the value ZAP1 of the point with the maximum value ZAP2 that appears second.
- the characteristic parameter PRM6 is a parameter relating to the size of the maximum point that appears third when the vertical axis acceleration is based on the third peak in the reference period, that is, the zero cross point that changes from negative to positive.
- PRM6 is a dimensionless amount ZAP3 / ZAP2 obtained by normalizing the value ZAP3 of the point with the maximum value ZAP2 that appears second.
- the characteristic parameter PRM7 exists between the second maximum and the third maximum appearing when the vertical axis acceleration is based on the second valley in the reference period, that is, the zero cross point that changes from negative to positive. This is a parameter related to the size of the local minimum point.
- PRM7 is a dimensionless quantity ZAV2 / ZAP2 obtained by normalizing the value ZAV2 of the point with the value ZAP2 of the maximum point that appears second.
- the characteristic parameter PRM8 is a parameter related to the appearance timing of the minimum point in the reference period for the vertical axis acceleration.
- PRM8 is a dimensionless amount T ZAMN / StepT obtained by normalizing the timing T ZAMN with a reference period Step T.
- the characteristic parameter PRM9 is a parameter related to the size of the minimum point in the reference period for the vertical axis acceleration.
- PRM8 is the value ZAMN (unit: [m / sec 2 ]) of the point.
- the characteristic parameter PRM10 relates to the difference between the first peak in the reference period with respect to the vertical axis acceleration, that is, the size of the maximum point that appears first when the zero cross point that changes from negative to positive is used as a reference, and the size of the minimum point in the same reference period. It is a parameter.
- PRM10 is ZAP1-ZAMN (unit: [m / sec 2 ]).
- the characteristic parameter PRM11 is the second peak in the reference period with respect to the vertical axis acceleration, that is, the size of the maximum point that appears second when the zero cross point that changes from negative to positive is used as the reference, and the size of the minimum point in the same reference period. It is a parameter related to the difference.
- PRM11 is ZAP2-ZAMN (unit: [m / sec 2 ]).
- the characteristic parameter PRM12 is the third peak in the reference period with respect to the vertical axis acceleration, that is, the size of the maximum point that appears third when the zero cross point that changes from negative to positive is used as the reference, and the size of the minimum point in the same reference period. It is a parameter related to the difference.
- the PRM 12 is ZAP3-ZAMN (unit: [m / sec 2 ]).
- the characteristic parameter PRM13 is a parameter relating to the appearance timing of the maximum point in the reference period for the first-order integral (vertical axis velocity) of the vertical axis acceleration.
- the PRM 13 is a dimensionless amount T ZVMX / StepT obtained by normalizing the timing T ZVMX with the reference period StepT.
- the characteristic parameter PRM14 is a parameter related to the timing at which 90% of the maximum point is reached for the first time in the reference period with respect to the first-order integral (vertical axis velocity) of the vertical axis acceleration.
- the PRM 14 is a dimensionless quantity T ZVMX90 / StepT obtained by normalizing the timing T ZVMX90 with the reference period StepT.
- the characteristic parameter PRM15 is a parameter related to the timing at which 80% of the maximum point is reached for the first time in the reference period for the first-order integral (vertical axis velocity) of the vertical axis acceleration.
- the PRM 15 is a dimensionless amount T ZVMX80 / StepT obtained by normalizing the timing T ZVMX80 with the reference period StepT.
- the characteristic parameter PRM16 is a parameter related to the timing at which 50% of the maximum point is reached for the first time in the reference period with respect to the first-order integral (vertical axis velocity) of the vertical axis acceleration.
- the PRM 16 is a dimensionless amount T ZVMX50 / StepT obtained by normalizing the timing TZVMX50 with the reference period StepT.
- the characteristic parameter PRM17 is a parameter related to the appearance timing of the maximum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration.
- the PRM 17 is a dimensionless amount TZTMX / StepT obtained by normalizing the timing TZTMX with the reference period StepT.
- the characteristic parameter PRM18 is a parameter related to the appearance timing of the minimum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration.
- the PRM 18 is a dimensionless amount T ZTMN / StepT obtained by normalizing the timing T ZTMN with a reference period Step T.
- the characteristic parameter PRM19 is a parameter related to the appearance timing of the zero point (orbit neutral point) in the reference period (not including the period boundary) for the second-order integral (vertical axis trajectory) of the vertical axis acceleration.
- the PRM 19 is a dimensionless amount T ZTNTL / Step T obtained by normalizing the timing T ZTNTL with a reference period Step T.
- the characteristic parameter PRM20 is a parameter related to the size of the maximum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration.
- the PRM 20 is the value ZTMX (unit: [m]) of the point.
- the characteristic parameter PRM21 is a parameter related to the size of the minimum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration.
- PRM21 is the value ZTMN (unit: [m]) of the point.
- the characteristic parameter PRM22 is a parameter related to the difference between the maximum point size in the reference period and the minimum point size in the same reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration.
- the PRM 22 is ZTMX-ZTMN (unit: [m]).
- the first feature parameter PRM1 to the twenty-second feature parameter PRM22 are feature parameters indicating the temporal change feature of the vertical axis acceleration.
- the characteristic parameter which shows the characteristic of the time change of a longitudinal-axis acceleration is demonstrated.
- the characteristic parameter PRM23 is a parameter related to the appearance timing of the maximum point that appears first when the first peak that appears after the minimum point in the reference period with respect to the longitudinal acceleration, that is, the zero cross point that turns from the minimum point negative to positive, is used as a reference.
- the PRM 23 is a dimensionless amount T XAP1 / StepT obtained by normalizing the timing T XAP1 with the reference period StepT.
- the characteristic parameter PRM24 is a parameter relating to the appearance timing of the maximum point appearing second when the peak appears second after the minimum point in the reference period with respect to the longitudinal acceleration, that is, when the zero cross point that turns from the minimum point negative to positive is used as a reference. It is.
- the PRM 24 is a dimensionless amount T XAP2 / StepT obtained by normalizing the timing T XAP2 with the reference period StepT.
- the characteristic parameter PRM25 is a parameter related to the appearance timing of the minimum point that appears in the reference period with respect to the longitudinal acceleration.
- the PRM 25 is a dimensionless amount T XAMN / StepT obtained by normalizing the timing T XAMN with the reference period StepT.
- the characteristic parameter PRM26 is a parameter relating to the size of the maximum point that appears first when the first peak that appears after the minimum point in the reference period with respect to the longitudinal acceleration, that is, the zero cross point that changes from the minimum point negative to positive, is used as a reference.
- the PRM 26 is the value XAP1 (unit: [m / sec 2 ]) of the point.
- the characteristic parameter PRM27 is a parameter relating to the size of the maximum point that appears second when the peak appears second after the minimum point in the reference period with respect to the longitudinal acceleration, that is, the zero cross point that turns from the minimum point negative to positive. It is.
- the PRM 27 is a dimensionless amount XAP2 / XAP1 obtained by normalizing the value XAP2 of the point with the value of the first maximum point.
- the characteristic parameter PRM28 is a parameter relating to the size of the minimum point in the reference period for the longitudinal acceleration.
- the PRM 28 is the value XAMN (unit: [m / sec 2 ]) of the point.
- the characteristic parameter PRM29 is a parameter related to the appearance timing of the maximum point in the reference period for the first-order integration (front-rear axis velocity) of the front-rear axis acceleration.
- the PRM 29 is a dimensionless amount T ZVMX / StepT obtained by normalizing the timing T XVMX with the reference period StepT.
- the characteristic parameter PRM30 is a parameter related to the timing at which 90% of the maximum point is reached for the first time in the reference period with respect to the first order integration (front-rear axis velocity) of the front-rear axis acceleration.
- the PRM 30 is a dimensionless amount T XVMX90 / StepT obtained by normalizing the timing T XVMX90 with the reference period StepT.
- the characteristic parameter PRM31 is a parameter related to the timing at which 80% of the maximum point is reached for the first time in the reference period with respect to the first order integral (front-rear axis velocity) of the front-rear axis acceleration.
- the PRM 31 is a dimensionless quantity T XVMX80 / StepT obtained by normalizing the timing T XVMX80 with the reference period StepT.
- the characteristic parameter PRM32 is a parameter related to the appearance timing of the minimum point in the reference period for the first order integral (front-rear axis velocity) of the front-rear axis acceleration.
- the PRM 32 is a dimensionless amount T XVMN / StepT obtained by normalizing the timing T XVMN with the reference period StepT.
- the characteristic parameter PRM33 is a parameter related to the timing at which 50% of the minimum point is reached for the first time in the reference period with respect to the first order integration (front-rear axis velocity) of the front-rear axis acceleration.
- the PRM 33 is a dimensionless amount T XVMN50 / StepT obtained by normalizing the timing T XVMN50 with the reference period StepT.
- step S46 the 33 types of feature parameters derived by the control unit 110 in the present embodiment.
- These feature parameters PRM1 to PRM33 are weighted and added in step S46 (FIG. 13) as will be described later.
- weighted addition for the characteristic parameters PRM1 to PRM33 represented by the following equations:
- the weighting factors (k 1 , k 2 ,..., K 32 , k 33 ) and constant terms k 0 (k 0 to k 33 are real constants including zero) used for weighted addition will be described later.
- steps S45a to S45c it is appropriately determined according to the stride and / or the reference period of the measurement subject.
- step S45a the control unit 110 uses the vertical axis trajectory time series data (FIG. 9) to evaluate the step width of the measurement subject. For example, the control unit 110 estimates the step length of the person to be measured based on the magnitude of the difference (ZTPP in FIG. 9) between the maximum value and the minimum value of the vertical axis trajectory in one reference period. The estimation may be performed using a conversion table (conversion formula) between the difference ZTPP between the maximum value and the minimum value of the vertical axis trajectory experimentally obtained in advance and the stride.
- conversion table conversion formula
- step S45b the control unit 110 obtains a reference period corresponding to one step in the walking cycle of the subject using the vertical axis acceleration time-series data (FIG. 7).
- the control part 110 should just obtain
- step S45c the control unit 110 determines a weighting coefficient used for calculating the average pelvic inclination angle ⁇ according to the above equation 1.
- FIG. 14 is a reference table used for determining the weight coefficient set.
- the control unit 110 determines whether the step length of the measurement subject is “small”, “medium”, or “large” based on the step width of the measurement subject evaluated in step S45a.
- the control unit 110 further determines whether the reference period is “short”, “medium”, or “long”, thereby using the weight coefficient set K mn to be used. (M: 1, 2, 3, n: 1, 2, 3) is determined.
- Each of the weighting coefficient sets K mn is a set including weighting coefficients (k 1 , k 2 ,..., K 32 , k 33 ) and a constant term k 0 used in Equation 1.
- the weighting coefficient set K 11 are stride relatively small and reference period the average pelvic tilt angle for the Formula 1 the group consisting of a number of the subject having the common feature that relatively short ⁇ can be calculated accurately Thus, it is a set of values of weight coefficients and constant terms obtained experimentally.
- the weight coefficient set K 33 is obtained from a large number of subjects having a common feature that the stride is relatively large and the reference period is relatively long so that the average pelvic inclination angle ⁇ can be accurately calculated by Equation 1. It is a set of values of weighting factors and constant terms obtained experimentally for a group.
- Each of the other weight coefficient sets is a set of weight coefficient values and constant term values obtained experimentally in the same manner.
- step S45a to S45c is processing performed to select a weighting factor used in the evaluation of Equation 1 when a plurality of weighting factor sets are provided according to the step length and / or the reference period of the measurement subject. . Therefore, when there is one weighting coefficient set used by the control unit 110, the processing from steps S45a to S45c may be skipped.
- step S46 the control unit 110 uses the feature parameters PRM1 to PRM33 calculated in step S44 and the weighting coefficient set K mn determined in step S45c to
- PRM i i: an integer from 1 to 33
- k mn, 0 and k mn, i i: an integer from 1 to 33
- K It is a constant term and a weighting factor included in mn .
- step S5 the control unit 110, based on the average pelvic inclination angle ⁇ obtained in step S4, indicates that the posture of the person to be measured is “forward tilt” or “middle”. Or “backward tilt” (see FIG. 5A).
- control unit 110 outputs (transmits) the evaluation result obtained in step S5 to the smartphone 200 (step S6).
- the smartphone 200 When the smartphone 200 receives the information from the activity meter 100, the smartphone 200 displays the evaluation result together with the total score on the display unit 240. For example, a message such as “Your walking posture has a tendency to lean forward” is displayed on the display unit 240 of the smartphone 200. It should be noted that the display unit 240 may display such that the forward / backward inclination can be intuitively understood using illustrations, animations, and the like. Further, the display unit 240 may display the average pelvic inclination angle ⁇ itself obtained in step S46 or an evaluation numerical value obtained by appropriately processing the value of ⁇ .
- the user can know whether the posture during walking is forward leaning or backward leaning.
- the inventor has obtained one or a plurality of feature quantities (feature parameters) that capture the characteristics of the waveform shape of the time-varying acceleration output from the acceleration sensor attached to the waist of the subject. It was experimentally found that the result of the weighted addition of has a good correlation with the inclination angle of the pelvis when the subject walks. According to this walking posture meter 1, it is possible to appropriately evaluate whether or not the posture of the measurement subject during walking is forward / backward. Moreover, since this walking posture meter 1 performs the evaluation based on the output of the acceleration sensor 112, it can be easily evaluated regardless of a large facility such as motion capture.
- 33 types of characteristic parameters PRM1 to PRM33 are calculated, but the present invention is not limited to this.
- the present invention is not limited to this.
- instead of calculating all 33 types of feature parameters PRM1 to PRM33 only a part is calculated, and only the feature parameters are used to evaluate the tendency of the subject to lean forward / backward during walking. May be.
- the weight coefficient set K mn is set in accordance with the step length and / or the reference period of the measurement subject, but the present invention is not limited to this.
- a walking cycle (corresponding to a period of two steps including the left foot reference period and the right foot reference period) is used, and the weighting coefficient sets K mn are respectively set according to the step length and / or the walking cycle of the measurement subject. May be set.
- the activity meter 100 and the smartphone 200 communicate with each other by BLE communication, but the present invention is not limited to this.
- the activity meter 100 and the smartphone 200 may communicate with each other when the smartphone 200 and the activity meter 100 approach each other by NFC (Near Field Communication).
- the walking posture meter of the present invention is configured as a system including the activity meter 100 and the smartphone 200.
- the present invention is not limited to this.
- the walking posture meter of the present invention may be configured with only the smartphone 200.
- the smartphone 200 includes an acceleration sensor.
- the memory 220 of the smartphone 200 has a program for quantitatively evaluating whether or not the human walking posture is a correct posture in the control unit 210, more specifically, the forward / backward tilt of the posture while walking. Install a program to assess the degree.
- the walking posture meter of this invention can be comprised compactly and compactly.
- the program can be recorded as application software on a recording medium such as a CD, a DVD, or a flash memory.
- the application software recorded on this recording medium is installed in a substantial computer device such as a smartphone, personal computer, PDA (Personal Digital Assistance), etc., so that the human walking posture is correct in these computer devices. It is possible to execute a method for quantitatively evaluating whether or not.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Geometry (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
Description
この発明は歩行姿勢計に関し、より詳しくは、ヒトの歩行姿勢が正しい姿勢であるか否かを定量的に評価する歩行姿勢計に関する。 The present invention relates to a walking posture meter, and more particularly to a walking posture meter that quantitatively evaluates whether or not a human walking posture is a correct posture.
また、この発明はヒトの歩行姿勢が正しい姿勢であるか否かを定量的に評価する方法をコンピュータに実行させるためのプログラムに関する。 The present invention also relates to a program for causing a computer to execute a method for quantitatively evaluating whether or not a human walking posture is a correct posture.
従来、ヒトの歩行中の姿勢を評価する装置としては、例えば特許文献1(特開2011-078728号公報)に示される装置がある。この装置は、ヒトが所定の姿勢をとったときにヒトの腰部に装着された加速度センサから出力にもとづき重力加速度の方向を検知することにより、当該姿勢をとったときの腰部の地面に対する傾き角を求め、これより骨盤の傾きを推定する。 Conventionally, as an apparatus for evaluating the posture of a human during walking, for example, there is an apparatus disclosed in Patent Document 1 (Japanese Patent Application Laid-Open No. 2011-078728). This device detects the direction of gravitational acceleration based on the output from an acceleration sensor attached to the human hip when the human takes a predetermined posture, so that the tilt angle of the hip with respect to the ground when the posture is taken To estimate the pelvic tilt.
また、例えば特許文献2(特開2011-251013号公報)に示される機器は、ヒトの腰部に装着された加速度センサの出力から移動量を算出し、当該移動量にもとづいて歩行の軌跡を取得する。 Further, for example, the device disclosed in Patent Document 2 (Japanese Patent Application Laid-Open No. 2011-251013) calculates a movement amount from an output of an acceleration sensor attached to a human waist and acquires a walking locus based on the movement amount. To do.
ヒトの歩行姿勢の正しさ(あるいは美しさ)を評価するポイントとして、背中がまっすぐ天に向かい、進行方向に向かって前のめり(前傾)でも及び腰(後傾)でもない中立的な姿勢で歩くことができているか否かという観点がある。 As a point to evaluate the correctness (or beauty) of a human walking posture, walk in a neutral posture with your back straight to the sky and neither forward leaning (forward leaning) nor waist (back leaning) in the direction of travel. There is a viewpoint of whether or not.
しかしながら、従来は、歩いているヒトの姿勢が上述した前傾姿勢か、後傾姿勢であるか、そのいずれでもない中立的な姿勢であるか、を手軽に定量的に評価する手段がなかった。 However, conventionally, there has been no means for easily and quantitatively evaluating whether the posture of the walking human is the above-described forward leaning posture, the backward leaning posture, or the neutral posture that is neither of them. .
以上を鑑み、本発明の一態様により、ヒトの歩行姿勢が進行方向に対して前傾・後傾しているか否かを手軽に評価できる歩行姿勢計が提供される。 In view of the above, according to one aspect of the present invention, there is provided a walking posture meter that can easily evaluate whether or not a human walking posture is tilted forward or backward with respect to the traveling direction.
また、本発明の別の一態様により、ヒトの歩行姿勢が進行方向に対して前傾・後傾しているか否かを手軽に評価する方法をコンピュータに実行させるためのプログラムが提供される。 Further, according to another aspect of the present invention, there is provided a program for causing a computer to execute a method for easily evaluating whether or not a human walking posture is tilted forward / backward with respect to the traveling direction.
上記課題を解決するため、本発明の一態様による歩行姿勢計は、被測定者の歩行姿勢を評価する歩行姿勢計であって、被測定者の腰の正中線上に装着される加速度センサと、前記加速度センサが出力する上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、前記被測定者の歩行中の姿勢の前後方向の傾きに対応する量を算出する演算部と、を備える。 In order to solve the above problems, a walking posture meter according to one aspect of the present invention is a walking posture meter that evaluates the walking posture of the measurement subject, and an acceleration sensor that is mounted on the midline of the measurement subject's waist; Using the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output from the acceleration sensor, an amount corresponding to the forward / backward inclination of the posture of the person being measured is calculated. And an arithmetic unit.
本明細書で、歩行中の姿勢の「前後方向の傾き」とは、歩行中のヒトの上半身(骨盤から上の半身)の進行方向に沿った傾きを意味する。当該傾きは、歩容が一巡する間に現れる任意の姿勢(脚位相)についてのものでよい。例えば、傾きは、前脚の踵が路面に接地する時点近傍における被測定者の姿勢の傾きについてのものでよい。また例えば、傾きは、遊脚である後脚が立脚である前脚に進行方向に沿って一致する時点近傍における被測定者の姿勢の傾きについてのものでよい。 In this specification, the “inclination in the front-rear direction” of the posture while walking means the inclination along the traveling direction of the upper body of the human being walking (the upper body from the pelvis). The inclination may be for an arbitrary posture (leg phase) that appears while the gait goes around. For example, the inclination may be about the inclination of the posture of the person to be measured in the vicinity of the time when the heel of the front leg contacts the road surface. Further, for example, the inclination may be about the inclination of the posture of the measurement subject in the vicinity of the time point where the rear leg, which is a free leg, coincides with the front leg, which is a standing leg, along the traveling direction.
本発明の一態様による歩行姿勢計では、加速度センサが腰の正中線上に装着される。演算部は、加速度センサが出力する上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、被測定者の歩行中の姿勢の前後方向の傾きの程度に対応する量を算出する。したがって、当該量に応じて、被測定者の歩行中の姿勢の前後方向の傾きの程度を評価することができる。また、この歩行姿勢計では、被測定者の腰の正中線上に装着された加速度センサの出力に基づいて上記したような評価が可能であるため、モーションキャプチャのような大がかりな設備によらずに手軽に評価することが可能である。 In the walking posture meter according to one aspect of the present invention, the acceleration sensor is mounted on the midline of the waist. The arithmetic unit uses the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output from the acceleration sensor to respond to the degree of inclination of the posture of the subject during walking in the longitudinal direction. The amount to be calculated is calculated. Therefore, according to the amount, it is possible to evaluate the degree of inclination in the front-rear direction of the posture of the measurement subject during walking. In addition, this walking posture meter can perform the above-described evaluation based on the output of the acceleration sensor mounted on the midline of the subject's waist, so that it does not depend on large-scale equipment such as motion capture. It is possible to evaluate easily.
一実施形態による歩行姿勢計では、前記演算部は、前記上下軸加速度の時間変化および前記前後軸加速度の時間変化の特徴を示す特徴パラメータを導出し、当該特徴パラメータに基づいて前記前後方向の傾きに対応する量を算出する、ことを特徴とする。 In the walking posture meter according to the embodiment, the calculation unit derives a characteristic parameter indicating characteristics of the temporal change of the vertical axis acceleration and the temporal change of the longitudinal axis acceleration, and the inclination in the front-rear direction based on the characteristic parameter. An amount corresponding to is calculated.
演算部は、上下軸加速度および前後軸加速度の時間変化の特徴を示す特徴パラメータを導出し、当該特徴パラメータに基づいて前後方向の傾きに対応する量を算出するので、被測定者の歩行中の姿勢の前後方向の傾きの程度に対応する量を迅速に算出することができる。 The calculation unit derives a feature parameter indicating the temporal change characteristics of the vertical axis acceleration and the longitudinal axis acceleration, and calculates an amount corresponding to the tilt in the front-rear direction based on the feature parameter. An amount corresponding to the degree of inclination of the posture in the front-rear direction can be quickly calculated.
一実施形態による歩行姿勢計では、前記演算部は、複数種類の特徴パラメータを導出し、導出した複数種類の特徴パラメータを重み付け加算して得られる量に基づいて前記前後方向の傾きに対応する量を算出する、ことを特徴とする。 In the walking posture meter according to one embodiment, the calculation unit derives a plurality of types of feature parameters, and an amount corresponding to the inclination in the front-rear direction based on an amount obtained by weighted addition of the derived types of feature parameters. Is calculated.
演算部は、導出した複数種類の特徴パラメータを重み付け加算することにより被測定者の歩行中の姿勢の前後方向の傾きの程度を評価することができる。したがって、評価に要する計算は至極単純であり、迅速な評価が可能となっている。 The calculation unit can evaluate the degree of inclination in the front-rear direction of the posture of the measurement subject while walking by weight-adding the plurality of derived characteristic parameters. Therefore, the calculation required for the evaluation is extremely simple, and quick evaluation is possible.
一実施形態による歩行姿勢計では、前記演算部は、前記上下軸加速度の時間的な変化波形に基づいて前記被測定者の歩幅を推定し、当該推定された歩幅に基づいて、前記重み付け加算に用いる重み係数を変化させる、ことを特徴とする。 In the walking posture meter according to one embodiment, the calculation unit estimates the step length of the subject based on a temporal change waveform of the vertical axis acceleration, and performs the weighted addition based on the estimated step length. The weighting coefficient to be used is changed.
演算部は、上下軸加速度の時間的な変化波形に基づいて被測定者の歩幅を推定し、当該推定された歩幅に基づいて、重み付け加算に用いる重み係数を変化させるため、被測定者の歩き方の特徴(歩幅の広狭)に因らず適切な評価ができる。 The calculation unit estimates the measured person's stride based on the temporal change waveform of the vertical axis acceleration, and changes the weighting coefficient used for the weighted addition based on the estimated stride. Appropriate evaluation is possible regardless of the characteristics of the direction (the width of the stride).
一実施形態による歩行姿勢計では、前記演算部は、前記上下軸加速度の時間的な変化波形および前記前後軸加速度の時間的な変化波形の少なくともいずれか一方に基づいて前記被測定者の歩行周期またはその歩行周期のうち1歩分に相当する基準期間を求め、当該求めた歩行周期または基準期間に基づいて、前記重み付け加算に用いる重み係数を変化させる、ことを特徴とする。 In the walking posture meter according to an embodiment, the calculation unit is configured to calculate the walking cycle of the subject based on at least one of the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration. Alternatively, a reference period corresponding to one step in the walking cycle is obtained, and a weighting coefficient used for the weighted addition is changed based on the obtained walking cycle or reference period.
演算部は、上下軸加速度および前後軸加速度の時間的な変化波形の少なくともいずれか一方に基づいて被測定者の歩行周期またはその歩行周期のうち1歩分に相当する基準期間を求め、求めた歩行周期または基準期間に基づいて、重み付け加算に用いる重み係数を変化させるため、被測定者の歩き方の特徴(歩行周期の長短)に因らず適切な評価ができる。 The calculation unit obtains a reference period corresponding to one step of the walking cycle of the measured person or the walking cycle based on at least one of the temporal change waveform of the vertical axis acceleration and the longitudinal axis acceleration. Since the weighting coefficient used for the weighted addition is changed based on the walking cycle or the reference period, appropriate evaluation can be performed regardless of the characteristics of the person to be measured (walking cycle length).
一実施形態による歩行姿勢計では、前記前後方向の傾きに対応する量は、歩行中の前記被測定者の骨盤と水平とがなす角度である、ことを特徴とする。 In the walking posture meter according to an embodiment, the amount corresponding to the inclination in the front-rear direction is an angle formed between the pelvis of the measurement subject and the horizontal direction during walking.
演算部は、歩行中の被測定者の骨盤の傾斜角(水平となす角度)を算出し、算出した傾斜角を用いて被測定者の歩行中の姿勢の前後方向の傾きを評価する。したがって、本歩行姿勢計は、極めて確度の高い評価をすることができる。 The calculation unit calculates the tilt angle (angle formed with the horizontal) of the person being measured while walking, and evaluates the tilt of the posture of the person being measured in the front-rear direction using the calculated tilt angle. Therefore, this walking posture meter can perform evaluation with extremely high accuracy.
本発明の別の一態様によるプログラムは、被測定者の歩行姿勢を評価する方法をコンピュータに実行させるためのプログラムであって、前記方法は、被測定者の腰の正中線上に装着された加速度センサの出力を取得するステップと、前記加速度センサが出力する上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、前記被測定者の歩行中の姿勢の前後方向の傾きに対応する量を算出するステップと、を備えたことを特徴とする。 A program according to another aspect of the present invention is a program for causing a computer to execute a method for evaluating a walking posture of a measured person, the method including an acceleration mounted on a midline of the waist of the measured person. Using the step of obtaining the sensor output and the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output by the acceleration sensor, Calculating an amount corresponding to the inclination of the direction.
当該プログラムを実行させることにより、コンピュータは、まず被測定者の腰の正中線上に装着された加速度センサの出力を取得する。そして、記加速度センサが出力する上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、被測定者の歩行中の姿勢の前後方向の傾きの程度に対応する量を算出する。したがって、当該量に応じて、歩行中の被測定者の姿勢の前後方向の傾きを評価することができる。また、このプログラムでは、被測定者の腰の正中線上に装着された加速度センサの出力に基づいて上記したような評価が可能であるため、モーションキャプチャのような大がかりな設備によらずに手軽に評価することが可能である。 By executing the program, the computer first acquires the output of the acceleration sensor mounted on the midline of the subject's waist. Then, using the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output by the acceleration sensor, it corresponds to the degree of the forward / backward inclination of the posture of the measurement subject during walking. Calculate the amount. Therefore, the inclination of the posture of the measurement subject during walking can be evaluated according to the amount. In addition, this program enables the above-mentioned evaluation based on the output of the acceleration sensor mounted on the midline of the subject's waist, so it is easy to use regardless of large equipment such as motion capture. It is possible to evaluate.
以上より明らかなように、この発明の一態様による歩行姿勢計によれば、被測定者の歩行中の姿勢の前後方向の傾きを手軽に評価することができる。 As is clear from the above, according to the walking posture meter according to one aspect of the present invention, it is possible to easily evaluate the inclination in the front-rear direction of the posture of the measurement subject during walking.
また、この発明の別の一態様によるプログラムをコンピュータに実行させることにより、被測定者の歩行中の姿勢の前後方向の傾きを手軽に評価することができる。 In addition, by causing a computer to execute a program according to another aspect of the present invention, it is possible to easily evaluate the tilt in the front-rear direction of the posture of the person being measured while walking.
以下、この発明の実施の形態を、図面を参照しながら詳細に説明する。 Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
図1は、この発明の一実施形態の歩行姿勢計(全体を符号1で示す。)のシステム構成を示している。この歩行姿勢計1は、活動量計100と、スマートフォン200とを含んでいる。活動量計100とスマートフォン200とは、この例ではBLE(Bluetooth low energy;低消費電力Bluetooth、Bluetooth Core Specification Ver. 4.0において規定。)通信によって互いに通信可能になっている。
FIG. 1 shows a system configuration of a walking posture meter (the whole is denoted by reference numeral 1) according to an embodiment of the present invention. This walking
図2に示すように、活動量計100は、ケーシング100Mと、このケーシング100Mに搭載された、制御部110と、発振部111と、加速度センサ112と、メモリ120と、操作部130と、表示部140と、BLE通信部180と、電源部190と、リセット部199とを含む。
As shown in FIG. 2, the
ケーシング100Mは、この活動量計100を携帯し易いように、ヒトの手のひらに収まる程度の大きさに形成されている。
The
発振部111は、水晶振動子を含み、この活動量計100の動作タイミングの基準となるクロック信号を発生する。発振部111は、クロックジェネレータとしての機能を有するモジュールチップでよい。
The
加速度センサ112は、ケーシング100Mが受ける3軸(3方向)の加速度をそれぞれ検出して、制御部110へ出力する。加速度センサ112は、3軸加速度センサのモジュールチップでよい。
The
メモリ120は、ROM(Read Only Memory)とRAM(Random Access Memory)とを含む。ROMは、この活動量計100を制御するためのプログラムのデータを記憶する。また、RAMは、この活動量計100の各種機能を設定するための設定データ、加速度測定結果および演算結果のデータなどを記憶する。
The
制御部110は、上記クロック信号に基づいて動作するCPU(Central Processing Unit;中央演算処理装置)を含み、メモリ120に記憶された活動量計100を制御するためのプログラムに従って、加速度センサ112からの検知信号に基づいて、この活動量計100の各部(メモリ120、表示部140およびBLE通信部180を含む。)を制御する。この制御部110は、少なくとも、上下軸加速度および前後軸加速度の時系列データを処理することができる信号処理系を含む。
The
操作部130は、この例ではボタンスイッチからなり、電源オン・オフ切り替えの操作、表示内容切り替えの操作など、適宜の操作入力を受け付ける。
The
表示部140は、この例ではLCD(液晶表示素子)または有機EL(エレクトロルミネッセンス)ディスプレイからなる表示画面を含み、この表示画面に制御部110から受けた信号に従って所定の情報を表示する。表示部140は、電源のオン/オフ、動作状態等を点灯、消灯、点滅、等で表示するLED(発光ダイオード)であってもよい。
In this example, the
電源部190は、この例ではボタン電池からなり、この活動量計100の各部へ電力を供給する。
The
BLE通信部180は、スマートフォン200との間でリアルタイムで通信を行う。例えば、スマートフォン200へ測定結果を表す情報などを送信する。また、スマートフォン200から操作指示を受信する。BLE通信部180は、BLE機能を備えたモジュールチップでよい。
The
リセット部199は、スイッチからなり、制御部110の動作やメモリ120の記憶内容をリセットして初期化する。
The
図3に示すように、スマートフォン200は、本体200Mと、この本体200Mに搭載された、制御部210と、メモリ220と、操作部230と、表示部240と、BLE通信部280と、ネットワーク通信部290とを含む。このスマートフォン200は、市販のスマートフォンに、活動量計100への指示を行わせるようにアプリケーションソフトウェア(コンピュータプログラム)をインストールしたものである。
As illustrated in FIG. 3, the
制御部210は、CPUおよびその補助回路を含み、スマートフォン200の各部を制御し、メモリ220に記憶されたプログラムおよびデータに従って処理を実行する。すなわち、操作部230、および、通信部280,290から入力されたデータを処理し、処理したデータを、メモリ220に記憶させたり、表示部240で表示させたり、通信部280,290から出力させたりする。
The
メモリ220は、制御部210でプログラムを実行するために必要な作業領域として用いられるRAMと、制御部210で実行するための基本的なプログラムを記憶するためのROMとを含む。また、メモリ220の記憶領域を補助するための補助記憶装置の記憶媒体として、半導体メモリ(メモリカード、SSD(Solid State Drive))などが用いられてもよい。
The
操作部230は、この例では、表示部240上に設けられたタッチパネルからなっている。なお、キーボードその他のハードウェア操作デバイスを含んでいてもよい。
The
表示部240は、表示画面(例えばLCDまたは有機ELディスプレイからなる)を含む。表示部240は、制御部210によって制御されて、所定の画像を表示画面に表示させる。
The
BLE通信部280は、活動量計100との間でリアルタイムで通信を行う。例えば、活動量計100へ操作指示を送信する。また、活動量計100から測定結果を表す情報などを受信する。
ネットワーク通信部290は、制御部210からの情報をネットワーク900を介して他の装置へ送信するとともに、他の装置からネットワーク900を介して送信されてきた情報を受信して制御部210に受け渡すことができる。
The
例えば図4(A)に示すように、この歩行姿勢計1が例えばユーザとしての被測定者90によって使用される場合、活動量計100が装着クリップ100C(図1中に示す)によって被測定者90の正中線91上の腰の背面側に装着される。
For example, as shown in FIG. 4 (A), when this walking
この例では、図4(B)に示すように、被測定者90にとって前後方向をX軸、左右方向をY軸、上下方向をZ軸とする。そして、活動量計100の加速度センサ112は、被測定者90が前方へ歩行するのに伴ってケーシング100Mが受けるX軸(前後軸)の加速度、Y軸(左右軸)の加速度、Z軸(上下軸)の加速度をそれぞれ出力するものとする。
In this example, as shown in FIG. 4B, for the
この歩行姿勢計1によって測定を行う場合、被測定者90は、活動量計100とスマートフォン200の電源をオンする。それとともに、スマートフォン200のアプリケーションソフトウェアを起動して、操作部230、BLE通信部280を介して、活動量計100へ測定スタートを指示する。
When measuring with this walking
その状態で、被測定者90は前方へ真っ直ぐに所定の歩数、この例では10歩だけ歩行する。そして、被測定者90は、スマートフォン200の操作部230、BLE通信部280を介して、活動量計100へ演算および演算結果の出力を指示する。
In this state, the person to be measured 90 walks straight forward by a predetermined number of steps, in this example, 10 steps. Then, the person under
すると、活動量計100の制御部110は演算部として働いて、後述する演算を行う。そして、その演算結果を表す情報をBLE通信部180を介して、スマートフォン200へ送信する。
Then, the
図12は、活動量計100の制御部110による動作フローを示している。活動量計100の制御部110は、電源がオンされると、ステップS1に示すように、スマートフォン200からの測定スタートの指示を待つ。スマートフォン200からの測定スタートの指示を受信すると(ステップS1でYES)、ステップS2に示すように、制御部110は、加速度センサ112による3軸方向加速度の出力を取得する。加速度センサ112の出力の取得は、この例では10歩分の加速度時系列データを含む期間として、予め定められた期間(例えば14秒間)だけ行われる。取得された加速度の時系列データは、メモリ120に一旦記憶される。次に、制御部110は、ステップS3に示すように、スマートフォン200からの測定スタートの指示を待つ。スマートフォン200からの演算の指示を受信すると(ステップS3でYES)、ステップS4に示すように、制御部110は、骨盤の傾斜角に対応する量の算出を行う。そして、ステップS5に示すように、制御部110は評価部として働いて、その演算結果(骨盤傾斜角推定結果)を用いて、被測定者の歩行中の姿勢の前後方向の傾きの程度を多段的に評価する。その後、ステップS6に示すように、その評価の結果をスマートフォン200へ出力(送信)する。なお、制御部110は、少なくとも一歩分の加速度時系列データが得られ次第、ステップS4の処理を実行してもよい。その場合、ステップS3として示した判断のステップは省略されてもよい。
FIG. 12 shows an operation flow by the
図5(A)は、骨盤傾斜角の定義を説明する図である。本図は、ヒトの腰部を横から見たときの模式図であって、骨盤PVが、三角形で表現されている。本明細書において、骨盤傾斜角θは、時計回転方向を正として、上前腸骨棘ASIS(ASIS: Anterior Superior Iliac Spine)から上後腸骨棘PSIS(PSIS: Posterior Superior Iliac Spine)に向かう線分が水平Hとなす角である。また、以下の説明では、+3度未満の骨盤傾斜角θを歩行姿勢「後傾」と対応付け、+3度以上+12度未満の骨盤傾斜角θを、歩行姿勢「真中」と対応付け、+12度以上の骨盤傾斜角θを、歩行姿勢「前傾」と対応付けすることとする。なお、上記境界値は一例であって、その他の境界値であってもよい。また、上記例では、姿勢を「前傾」「真中」「後傾」の三段階に区分けしているが、その他の段数であってもよい。 FIG. 5 (A) is a diagram for explaining the definition of the pelvic inclination angle. This figure is a schematic view of the human waist when viewed from the side, and the pelvis PV is represented by a triangle. In this specification, the pelvic inclination angle θ is a line from the superior anterior iliac spine ASIS (ASIS: Anterior Superior Iliac Spine) to the superior posterior iliac spine PSIS (PSIS: Posterior Superior Iliac Spine) with the clockwise rotation direction being positive. The angle between the minute and the horizontal H. In the following description, a pelvic inclination angle θ of less than +3 degrees is associated with the walking posture “backward inclination”, and a pelvic inclination angle θ of +3 degrees or more and less than +12 degrees is associated with the walking posture “middle”, and +12 degrees. The above pelvic tilt angle θ is associated with the walking posture “forward tilt”. The boundary value is an example, and other boundary values may be used. In the above example, the posture is divided into three stages of “forward leaning”, “middle”, and “backward leaning”, but other stages may be used.
図5(B)~(D)は、それぞれ、ヒトの歩行姿勢「後傾」、「真中」、「前傾」の例を示す模式図である。図5(B)を参照すれば判るように、歩行姿勢「後傾」で歩く傾向を有するヒトは、反っくり返ったような姿勢で歩行する。(ほかに、「後傾」のヒトには、極めて強い猫背姿勢で歩行するヒトもいる。)図5(C)は、歩行姿勢「真中」の例図である。このように、歩行姿勢「真中」で歩くヒトは、背中がまっすぐに上方に伸びており、正しく(美しく)歩くことができている。図5(D)を参照すれば、歩行姿勢「前傾」で歩く傾向を有するヒトは、前かがみでうつむき気味の姿勢で歩行する。(ほかに、「前傾」のヒトには、おなかを突き出した姿勢で歩行するヒトもいる。) FIGS. 5B to 5D are schematic diagrams showing examples of human walking postures “backward tilt”, “middle”, and “forward tilt”, respectively. As can be seen from FIG. 5B, a human who has a tendency to walk in the walking posture “backward tilting” walks in a posture that turns upside down. (In addition, some people who are “backward tilted” walk with a very stooped posture.) FIG. 5C is an example of the walking posture “middle”. In this way, a human walking in the walking posture “middle” has a straight back and is able to walk correctly (beautifully). Referring to FIG. 5D, a human who has a tendency to walk in a walking posture “forward leaning” walks in a leaning posture with forward leaning. (In addition, some people who are “forward leaning” walk with their stomachs protruding.)
図6は、ヒトの歩容と、歩行周期のうち1歩分に相当する基準期間(図中のT7(=StepT))の間に、腰に装着された活動量計100の加速度センサ112から出力される上下軸加速度(鉛直上方を正とするZ軸方向加速度)の時間変化波形の典型例との関係を示す図である。
FIG. 6 illustrates the
繰り出された前足(図では右足)の踵が移動面と接地するタイミング(踵接地タイミング)の近傍において、上下軸加速度は、ゼロクロス点を通過して負から正へ転じる。 In the vicinity of the timing when the heel of the extended forefoot (right foot in the figure) comes into contact with the moving surface (the heel contact timing), the vertical axis acceleration changes from negative to positive through the zero cross point.
その後、上下軸加速度には、3つのピーク(極大点)(P1(時間t=T1)、P2(時間t=T3)、P3(時間t=T5))およびその間の谷(極小点)(V1(時間t=T2)、V2(時間t=T4))が現れる。歩容における、立脚(図では右足)と遊脚(図では左脚)とが進行方向に関して略一致するタイミング(立脚中期タイミング)は、第3のピークP3が現れたタイミング近傍に対応する。 Thereafter, the vertical axis acceleration includes three peaks (maximum points) (P 1 (time t = T 1 ), P 2 (time t = T 3 ), P 3 (time t = T 5 )) and a valley between them. (Minimum point) (V 1 (time t = T 2 ), V 2 (time t = T 4 )) appears. In the gait, the timing at which the stance (right foot in the figure) and the free leg (left leg in the figure) substantially coincide with each other in the traveling direction (mid-stance stance timing) corresponds to the vicinity of the timing at which the third peak P3 appears.
歩容における立脚中期タイミングを超えると、上下軸加速度は、再びゼロクロス点を通過して正から負に転じ、最小点(V3(時間t=T6))を通過し、やがて時間t=T7において再度ゼロクロス点(時間t=T7)を通過して負から正へ転じる。時間t=T7におけるゼロクロス点は、(図では左足を前足とする)次の一歩の踵接地タイミングである。 When the mid-stance timing of the gait in the gait is exceeded, the vertical axis acceleration passes through the zero cross point again, changes from positive to negative, passes through the minimum point (V 3 (time t = T 6 )), and eventually time t = T. 7 again passes the zero cross point (time t = T 7 ) and turns from negative to positive. The zero-cross point at time t = T 7 is the next one-step heel contact timing (with the left foot as the front foot in the figure).
このように、上下軸加速度には、ヒトの歩行の一歩の間に図示して説明したような波形が現れる。本明細書では、前足の踵が接地したタイミング(踵接地タイミング)から次の踵接地タイミングまでの期間(StepT)を、基準期間として規定する。なお、以下の説明において特に区別する必要がある場合に限り、左足の踵接地タイミングから右足の踵接地タイミングまでの期間を左足基準期間と称し、右足の踵接地タイミングから左足の踵接地タイミングまでの期間を右足基準期間と称することとして、左足による一歩と右足による一歩のそれぞれの基準期間を区別する。 Thus, in the vertical axis acceleration, a waveform as illustrated and described appears during one step of human walking. In this specification, a period (Step T) from the timing when the forefoot heel contacts the ground (heel contact timing) to the next heel contact timing is defined as a reference period. Note that the period from the left foot heel contact timing to the right foot heel contact timing is referred to as the left foot reference period, and from the right foot heel contact timing to the left foot heel contact timing only when there is a particular need to distinguish in the following description. By referring to the period as the right foot reference period, the reference periods of one step with the left foot and one step with the right foot are distinguished.
上方を正とする前記上下軸加速度の時間的な変化波形においては、加速度値が負から正へ変化するゼロクロス点の出現タイミングから次の負から正へ変化するゼロクロス点の出現タイミングまでの期間が一基準期間に相当する。 In the temporal change waveform of the vertical axis acceleration with the upper side being positive, there is a period from the appearance timing of the zero cross point at which the acceleration value changes from negative to positive to the appearance timing of the next zero cross point at which the negative value changes from positive to positive. It corresponds to one reference period.
上述の制御部110は、図12のステップS4の処理(骨盤傾斜角推定処理)において、少なくとも一基準期間にわたる上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、被測定者の歩行中の姿勢の前後方向の傾きに対応する量を算出する。ここでは、進行方向の前傾/後傾の程度に対応する量として、踵接地タイミングにおける被測定者の骨盤傾斜角θ1と、同一基準期間における立脚中期タイミングにおける被測定者の骨盤傾斜角θ2との平均骨盤傾斜角θ(θ=(θ1+θ2)/2)を用いる。なお、前後方向の傾きの程度に対応する量は、平均骨盤傾斜角θに限定されるものではない。
The
以下、図12のステップS4の処理(骨盤傾斜角推定処理)について、図7~図11のグラフ、図13のフローチャート、および、図14の参照表を参照して詳細に説明する。図13のフローチャートの各ステップは、上述の制御部110により実行されればよい。なお、当該各ステップの実行主体は、制御部110に代えて、制御部210であってもよい。
Hereinafter, the processing (pelvic tilt angle estimation processing) in step S4 in FIG. 12 will be described in detail with reference to the graphs in FIGS. 7 to 11, the flowchart in FIG. 13, and the reference table in FIG. Each step of the flowchart of FIG. 13 may be executed by the
図13のフローチャートを参照すれば、ステップS41において、制御部110は、ステップS2において取得した三軸加速度の時間変化波形から、上下軸(Z軸)加速度時系列データ(図7)および前後軸(X軸)加速度時系列データ(図10)を生成し、生成した上下軸加速度時系列データから負から正へ転じるゼロクロス点のタイミング(図7における下横軸の時間t=約6.2(上横軸における時間t=0))を検出し、(当該ゼロクロス点が踵接地タイミングと略一致するという実験的事実に基づいて)検出したタイミングを踵接地タイミングとして特定する。制御部110は、次の負から正へ転じるゼロクロス点を同様にして検出し、当該タイミングを次の一歩の基準期間の始期として(現基準期間の周期)特定する。これにより、基準期間(StepT)が確定される。
Referring to the flowchart of FIG. 13, in step S41, the
次に、ステップS42において、制御部110は、上下軸加速度時系列データ(図7)から上下軸速度時系列データ(図8)および上下軸軌道時系列データ(図9)を算出する。
Next, in step S42, the
次に、ステップS43において、制御部110は、前後軸加速度時系列データ(図10)から前後軸速度時系列データ(図11)を算出する。
Next, in step S43, the
次に、ステップS44において、制御部110は、上下軸加速度の時間変化および前後軸加速度の時間変化の特徴を示す特徴パラメータを導出する。制御部110は、上下軸加速度の時系列データ、上下軸速度の時系列データ、および、上下軸軌道の時系列データの少なくともいずれか1つを用いて上下軸加速度の時間変化の特徴を示す特徴パラメータを導出する。同様に、制御部110は、前後軸加速度の時系列データ、および、前後軸速度の時系列データの少なくともいずれか1つを用いて前後軸加速度の時間変化の特徴を示す特徴パラメータを導出する。ステップS44において導出すべき対象は、各軸加速度の時間変化の特徴を表す特徴量であるが、特徴量によっては、各軸速度、軌道の時間変化の時系列データを用いる方が導出容易である場合がある。そのため、制御部110は、各軸の速度の時系列データや軌道の時系列データを用いて、各軸の加速度の時間変化の特徴を表す量(特徴パラメータ)を導出する。
Next, in step S44, the
本実施の形態においては、特徴パラメータは33種(特徴パラメータPRM1~PRM33、即ち、PRMi(i:1~33の整数。))存在する。以下、各特徴パラメータについて詳細に説明する。なお、特徴パラメータの種類の数、および、各特徴パラメータの定義は、一例であって、異なる種類数および定義の特徴パラメータを用いてもよい。 In the present embodiment, there are 33 types of characteristic parameters (characteristic parameters PRM1 to PRM33, that is, PRMi (i: an integer of 1 to 33)). Hereinafter, each feature parameter will be described in detail. Note that the number of types of feature parameters and the definition of each feature parameter are examples, and different types of features and definitions of feature parameters may be used.
先ず、図7を参照し、第1の特徴パラメータPRM1から第12の特徴パラメータPRM12までの12種類の特徴パラメータについて説明する。 First, with reference to FIG. 7, 12 types of feature parameters from the first feature parameter PRM1 to the twelfth feature parameter PRM12 will be described.
第1の特徴パラメータPRM1(Z軸加速度第1極大点タイミング比)について:
特徴パラメータPRM1は、上下軸加速度について基準期間における最初のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき最初に現れる極大点の出現タイミングに関するパラメータである。例えば、PRM1は、当該タイミングTZAP1を基準期間StepTで規格化することで得られる無次元量TZAP1/StepTである。
Regarding the first characteristic parameter PRM1 (Z-axis acceleration first maximum point timing ratio):
The characteristic parameter PRM1 is a parameter related to the appearance timing of the maximum point that appears first when the vertical peak acceleration is based on the first peak in the reference period, that is, the zero cross point that changes from negative to positive. For example, PRM1 is a dimensionless quantity T ZAP1 / StepT obtained by normalizing the timing T ZAP1 with a reference period StepT.
第2の特徴パラメータPRM2(Z軸加速度第2極大点タイミング比)について:
特徴パラメータPRM2は、上下軸加速度について基準期間における第2のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき2番目に現れる極大点の出現タイミングに関するパラメータである。例えば、PRM2は、当該タイミングTZAP2を基準期間StepTで規格化することで得られる無次元量TZAP2/StepTである。
Regarding the second characteristic parameter PRM2 (Z-axis acceleration second maximum point timing ratio):
The characteristic parameter PRM2 is a parameter related to the appearance timing of the second maximum peak that appears when the vertical peak acceleration is based on the second peak in the reference period, that is, the zero cross point that changes from negative to positive. For example, PRM2 is a dimensionless quantity T ZAP2 / StepT obtained by normalizing the timing T ZAP2 with a reference period Step T.
第3の特徴パラメータPRM3(Z軸加速度第3極大点タイミング比)について:
特徴パラメータPRM3は、上下軸加速度について基準期間における第3のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき3番目に現れる極大点の出現タイミングに関するパラメータである。例えば、PRM3は、当該タイミングTZAP3を基準期間StepTで規格化することで得られる無次元量TZAP3/StepTである。
Regarding the third characteristic parameter PRM3 (Z-axis acceleration third maximum point timing ratio):
The characteristic parameter PRM3 is a parameter relating to the appearance timing of the third maximum point that appears when the vertical peak acceleration is based on the third peak in the reference period, that is, the zero cross point that changes from negative to positive. For example, PRM3 is a dimensionless amount T ZAP3 / StepT obtained by normalizing the timing T ZAP3 with a reference period StepT.
第4の特徴パラメータPRM4(Z軸加速度第2極小点タイミング比)について:
特徴パラメータPRM4は、上下軸加速度について基準期間における第2の谷、即ち、負から正に転じるゼロクロス点を基準としたとき2番目に現れる極大点と3番目に現れる極大点との間に存在する極小点の出現タイミングに関するパラメータである。例えば、PRM4は、当該タイミングTZAV2を基準期間StepTで規格化することで得られる無次元量TZAV2/StepTである。
Regarding the fourth characteristic parameter PRM4 (Z-axis acceleration second minimum point timing ratio):
The characteristic parameter PRM4 exists between the second maximum point that appears second and the maximum point that appears third when the vertical valley acceleration is based on the second valley in the reference period, that is, the zero cross point that changes from negative to positive. This is a parameter related to the appearance timing of the minimum point. For example, PRM4 is a dimensionless amount T ZAV2 / StepT obtained by normalizing the timing T ZAV2 with the reference period Step T.
第5の特徴パラメータPRM5(Z軸加速度第1極大点振幅比)について:
特徴パラメータPRM5は、上下軸加速度について基準期間における最初のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき最初に現れる極大点の大きさに関するパラメータである。例えば、PRM5は、当該点の値ZAP1を2番目に現れる極大点の値ZAP2で規格化することで得られる無次元量ZAP1/ZAP2である。
About the fifth characteristic parameter PRM5 (Z-axis acceleration first maximum point amplitude ratio):
The characteristic parameter PRM5 is a parameter related to the size of the maximum point that appears first when the vertical peak acceleration is based on the first peak in the reference period, that is, the zero cross point that changes from negative to positive. For example, PRM5 is a dimensionless quantity ZAP1 / ZAP2 obtained by normalizing the value ZAP1 of the point with the maximum value ZAP2 that appears second.
第6の特徴パラメータPRM6(Z軸加速度第3極大点振幅比)について:
特徴パラメータPRM6は、上下軸加速度について基準期間における第3のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき3番目に現れる極大点の大きさに関するパラメータである。例えば、PRM6は、当該点の値ZAP3を2番目に現れる極大点の値ZAP2で規格化することで得られる無次元量ZAP3/ZAP2である。
About the sixth characteristic parameter PRM6 (Z-axis acceleration third maximum point amplitude ratio):
The characteristic parameter PRM6 is a parameter relating to the size of the maximum point that appears third when the vertical axis acceleration is based on the third peak in the reference period, that is, the zero cross point that changes from negative to positive. For example, PRM6 is a dimensionless amount ZAP3 / ZAP2 obtained by normalizing the value ZAP3 of the point with the maximum value ZAP2 that appears second.
第7の特徴パラメータPRM7(Z軸加速度第2極小点振幅比)について:
特徴パラメータPRM7は、上下軸加速度について基準期間における第2の谷、即ち、負から正に転じるゼロクロス点を基準としたときに2番目に現れる極大点と3番目に現れる極大点との間に存在する極小点の大きさに関するパラメータである。例えば、PRM7は、当該点の値ZAV2を2番目に現れる極大点の値ZAP2で規格化することで得られる無次元量ZAV2/ZAP2である。
Regarding the seventh characteristic parameter PRM7 (Z-axis acceleration second minimum point amplitude ratio):
The characteristic parameter PRM7 exists between the second maximum and the third maximum appearing when the vertical axis acceleration is based on the second valley in the reference period, that is, the zero cross point that changes from negative to positive. This is a parameter related to the size of the local minimum point. For example, PRM7 is a dimensionless quantity ZAV2 / ZAP2 obtained by normalizing the value ZAV2 of the point with the value ZAP2 of the maximum point that appears second.
第8の特徴パラメータPRM8(Z軸加速度最小点タイミング比)について:
特徴パラメータPRM8は、上下軸加速度について基準期間における最小点の出現タイミングに関するパラメータである。例えば、PRM8は、当該タイミングTZAMNを基準期間StepTで規格化することで得られる無次元量TZAMN/StepTである。
Regarding the eighth feature parameter PRM8 (Z-axis acceleration minimum point timing ratio):
The characteristic parameter PRM8 is a parameter related to the appearance timing of the minimum point in the reference period for the vertical axis acceleration. For example, PRM8 is a dimensionless amount T ZAMN / StepT obtained by normalizing the timing T ZAMN with a reference period Step T.
第9の特徴パラメータPRM9(Z軸加速度最小点値)について:
特徴パラメータPRM9は、上下軸加速度について基準期間における最小点の大きさに関するパラメータである。例えば、PRM8は、当該点の値ZAMN(単位は[m/sec2])である。
Regarding the ninth feature parameter PRM9 (Z-axis acceleration minimum point value):
The characteristic parameter PRM9 is a parameter related to the size of the minimum point in the reference period for the vertical axis acceleration. For example, PRM8 is the value ZAMN (unit: [m / sec 2 ]) of the point.
第10の特徴パラメータPRM10(Z軸加速度第1極大点最小点差)について:
特徴パラメータPRM10は、上下軸加速度について基準期間における最初のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき最初に現れる極大点の大きさと同一基準期間における最小点の大きさの差に関するパラメータである。例えば、PRM10は、ZAP1-ZAMN(単位は[m/sec2])である。
About the tenth characteristic parameter PRM10 (Z-axis acceleration first maximum point minimum point difference):
The characteristic parameter PRM10 relates to the difference between the first peak in the reference period with respect to the vertical axis acceleration, that is, the size of the maximum point that appears first when the zero cross point that changes from negative to positive is used as a reference, and the size of the minimum point in the same reference period. It is a parameter. For example, PRM10 is ZAP1-ZAMN (unit: [m / sec 2 ]).
第11の特徴パラメータPRM11(Z軸加速度第2極大点最小点差)について:
特徴パラメータPRM11は、上下軸加速度について基準期間における第2のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき2番目に現れる極大点の大きさと同一基準期間における最小点の大きさの差に関するパラメータである。例えば、PRM11は、ZAP2-ZAMN(単位は[m/sec2])である。
Regarding the eleventh characteristic parameter PRM11 (Z-axis acceleration second maximum point minimum point difference):
The characteristic parameter PRM11 is the second peak in the reference period with respect to the vertical axis acceleration, that is, the size of the maximum point that appears second when the zero cross point that changes from negative to positive is used as the reference, and the size of the minimum point in the same reference period. It is a parameter related to the difference. For example, PRM11 is ZAP2-ZAMN (unit: [m / sec 2 ]).
第12の特徴パラメータPRM12(Z軸加速度第3極大点最小点差)について:
特徴パラメータPRM12は、上下軸加速度について基準期間における第3のピーク、即ち、負から正に転じるゼロクロス点を基準としたとき3番目に現れる極大点の大きさと同一基準期間における最小点の大きさの差に関するパラメータである。例えば、PRM12は、ZAP3-ZAMN(単位は[m/sec2])である。
Regarding the twelfth characteristic parameter PRM12 (Z-axis acceleration third maximum point minimum point difference):
The characteristic parameter PRM12 is the third peak in the reference period with respect to the vertical axis acceleration, that is, the size of the maximum point that appears third when the zero cross point that changes from negative to positive is used as the reference, and the size of the minimum point in the same reference period. It is a parameter related to the difference. For example, the PRM 12 is ZAP3-ZAMN (unit: [m / sec 2 ]).
次に、図8を参照し、第13の特徴パラメータPRM13から第16の特徴パラメータPRM16までの4種類の特徴パラメータについて説明する。 Next, the four types of feature parameters from the thirteenth feature parameter PRM13 to the sixteenth feature parameter PRM16 will be described with reference to FIG.
第13の特徴パラメータPRM13(Z軸速度最大点タイミング比)について:
特徴パラメータPRM13は、上下軸加速度の一階積分(上下軸速度)についての基準期間における最大点の出現タイミングに関するパラメータである。例えば、PRM13は、当該タイミングTZVMXを基準期間StepTで規格化することで得られる無次元量TZVMX/StepTである。
About the thirteenth characteristic parameter PRM13 (Z-axis speed maximum point timing ratio):
The characteristic parameter PRM13 is a parameter relating to the appearance timing of the maximum point in the reference period for the first-order integral (vertical axis velocity) of the vertical axis acceleration. For example, the PRM 13 is a dimensionless amount T ZVMX / StepT obtained by normalizing the timing T ZVMX with the reference period StepT.
第14の特徴パラメータPRM14(Z軸速度最大点前90%タイミング比)について:
特徴パラメータPRM14は、上下軸加速度の一階積分(上下軸速度)について基準期間において初めて最大点の90%に到達したタイミングに関するパラメータである。例えば、PRM14は、当該タイミングTZVMX90を基準期間StepTで規格化することで得られる無次元量TZVMX90/StepTである。
About the 14th characteristic parameter PRM14 (90% timing ratio before the maximum Z-axis speed):
The characteristic parameter PRM14 is a parameter related to the timing at which 90% of the maximum point is reached for the first time in the reference period with respect to the first-order integral (vertical axis velocity) of the vertical axis acceleration. For example, the PRM 14 is a dimensionless quantity T ZVMX90 / StepT obtained by normalizing the timing T ZVMX90 with the reference period StepT.
第15の特徴パラメータPRM15(Z軸速度最大点前80%タイミング比)について:
特徴パラメータPRM15は、上下軸加速度の一階積分(上下軸速度)について基準期間において初めて最大点の80%に到達したタイミングに関するパラメータである。例えば、PRM15は、当該タイミングTZVMX80を基準期間StepTで規格化することで得られる無次元量TZVMX80/StepTである。
Regarding the fifteenth characteristic parameter PRM15 (80% timing ratio before the maximum Z-axis speed):
The characteristic parameter PRM15 is a parameter related to the timing at which 80% of the maximum point is reached for the first time in the reference period for the first-order integral (vertical axis velocity) of the vertical axis acceleration. For example, the PRM 15 is a dimensionless amount T ZVMX80 / StepT obtained by normalizing the timing T ZVMX80 with the reference period StepT.
第16の特徴パラメータPRM16(Z軸速度最大点前50%タイミング比)について:
特徴パラメータPRM16は、上下軸加速度の一階積分(上下軸速度)について基準期間において初めて最大点の50%に到達したタイミングに関するパラメータである。例えば、PRM16は、当該タイミングTZVMX50を基準期間StepTで規格化することで得られる無次元量TZVMX50/StepTである。
Regarding the sixteenth characteristic parameter PRM16 (50% timing ratio before the maximum Z-axis speed):
The characteristic parameter PRM16 is a parameter related to the timing at which 50% of the maximum point is reached for the first time in the reference period with respect to the first-order integral (vertical axis velocity) of the vertical axis acceleration. For example, the PRM 16 is a dimensionless amount T ZVMX50 / StepT obtained by normalizing the timing TZVMX50 with the reference period StepT.
次に、図9を参照し、第17の特徴パラメータPRM17から第22の特徴パラメータPRM22までの6種類の特徴パラメータについて説明する。 Next, with reference to FIG. 9, six types of feature parameters from the 17th feature parameter PRM17 to the 22nd feature parameter PRM22 will be described.
第17の特徴パラメータPRM17(Z軸軌道最大点タイミング比)について:
特徴パラメータPRM17は、上下軸加速度の二階積分(上下軸軌道)についての基準期間における最大点の出現タイミングに関するパラメータである。例えば、PRM17は、当該タイミングTZTMXを基準期間StepTで規格化することで得られる無次元量TZTMX/StepTである。
Regarding the seventeenth characteristic parameter PRM17 (Z-axis trajectory maximum point timing ratio):
The characteristic parameter PRM17 is a parameter related to the appearance timing of the maximum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration. For example, the PRM 17 is a dimensionless amount TZTMX / StepT obtained by normalizing the timing TZTMX with the reference period StepT.
第18の特徴パラメータPRM18(Z軸軌道最小点タイミング比)について:
特徴パラメータPRM18は、上下軸加速度の二階積分(上下軸軌道)についての基準期間における最小点の出現タイミングに関するパラメータである。例えば、PRM18は、当該タイミングTZTMNを基準期間StepTで規格化することで得られる無次元量TZTMN/StepTである。
About the 18th characteristic parameter PRM18 (Z-axis trajectory minimum point timing ratio):
The characteristic parameter PRM18 is a parameter related to the appearance timing of the minimum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration. For example, the PRM 18 is a dimensionless amount T ZTMN / StepT obtained by normalizing the timing T ZTMN with a reference period Step T.
第19の特徴パラメータPRM19(Z軸軌道中立点タイミング比)について:
特徴パラメータPRM19は、上下軸加速度の二階積分(上下軸軌道)について(期間境界を含まない)基準期間におけるゼロ点(軌道中立点)の出現タイミングに関するパラメータである。例えば、PRM19は、当該タイミングTZTNTLを基準期間StepTで規格化することで得られる無次元量TZTNTL/StepTである。
About the 19th characteristic parameter PRM19 (Z-axis trajectory neutral point timing ratio):
The characteristic parameter PRM19 is a parameter related to the appearance timing of the zero point (orbit neutral point) in the reference period (not including the period boundary) for the second-order integral (vertical axis trajectory) of the vertical axis acceleration. For example, the PRM 19 is a dimensionless amount T ZTNTL / Step T obtained by normalizing the timing T ZTNTL with a reference period Step T.
第20の特徴パラメータPRM20(Z軸軌道最大点値)について:
特徴パラメータPRM20は、上下軸加速度の二階積分(上下軸軌道)についての基準期間における最大点の大きさに関するパラメータである。例えば、PRM20は、当該点の値ZTMX(単位は[m])である。
Regarding the twentieth feature parameter PRM20 (Z-axis trajectory maximum point value):
The characteristic parameter PRM20 is a parameter related to the size of the maximum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration. For example, the PRM 20 is the value ZTMX (unit: [m]) of the point.
第21の特徴パラメータPRM21(Z軸軌道最小点値)について:
特徴パラメータPRM21は、上下軸加速度の二階積分(上下軸軌道)についての基準期間における最小点の大きさに関するパラメータである。例えば、PRM21は、当該点の値ZTMN(単位は[m])である。
Regarding the twenty-first feature parameter PRM21 (Z-axis trajectory minimum point value):
The characteristic parameter PRM21 is a parameter related to the size of the minimum point in the reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration. For example, PRM21 is the value ZTMN (unit: [m]) of the point.
第22の特徴パラメータPRM22(Z軸軌道最大点最小点差)について:
特徴パラメータPRM22は、上下軸加速度の二階積分(上下軸軌道)について基準期間における最大点の大きさと同一基準期間における最小点の大きさの差に関するパラメータである。例えば、PRM22は、ZTMX-ZTMN(単位は[m])である。
Regarding the twenty-second feature parameter PRM22 (Z-axis trajectory maximum point minimum point difference):
The characteristic parameter PRM22 is a parameter related to the difference between the maximum point size in the reference period and the minimum point size in the same reference period for the second-order integral (vertical axis trajectory) of the vertical axis acceleration. For example, the PRM 22 is ZTMX-ZTMN (unit: [m]).
以上、第1の特徴パラメータPRM1から第22の特徴パラメータPRM22は、上下軸加速度の時間変化の特徴を示す特徴パラメータである。以下では、前後軸加速度の時間変化の特徴を示す特徴パラメータについて説明する。 As described above, the first feature parameter PRM1 to the twenty-second feature parameter PRM22 are feature parameters indicating the temporal change feature of the vertical axis acceleration. Below, the characteristic parameter which shows the characteristic of the time change of a longitudinal-axis acceleration is demonstrated.
図10を参照し、第23の特徴パラメータPRM23から第28の特徴パラメータPRM28までの6種類の特徴パラメータについて説明する。 Referring to FIG. 10, the six types of feature parameters from the 23rd feature parameter PRM23 to the 28th feature parameter PRM28 will be described.
第23の特徴パラメータPRM23(X軸加速度第1極大点タイミング比)について:
特徴パラメータPRM23は、前後軸加速度について基準期間において最小点の後で現れる最初のピーク、即ち、最小点負から正に転じるゼロクロス点を基準としたとき最初に現れる極大点の出現タイミングに関するパラメータである。例えば、PRM23は、当該タイミングTXAP1を基準期間StepTで規格化することで得られる無次元量TXAP1/StepTである。
About the 23rd characteristic parameter PRM23 (X-axis acceleration first maximum point timing ratio):
The characteristic parameter PRM23 is a parameter related to the appearance timing of the maximum point that appears first when the first peak that appears after the minimum point in the reference period with respect to the longitudinal acceleration, that is, the zero cross point that turns from the minimum point negative to positive, is used as a reference. . For example, the PRM 23 is a dimensionless amount T XAP1 / StepT obtained by normalizing the timing T XAP1 with the reference period StepT.
第24の特徴パラメータPRM24(X軸加速度第2極大点タイミング比)について:
特徴パラメータPRM24は、前後軸加速度について基準期間において最小点の後で2番目に現れるピーク、即ち、最小点負から正に転じるゼロクロス点を基準としたとき2番目に現れる極大点の出現タイミングに関するパラメータである。例えば、PRM24は、当該タイミングTXAP2を基準期間StepTで規格化することで得られる無次元量TXAP2/StepTである。
Regarding the 24th characteristic parameter PRM24 (X-axis acceleration second maximum point timing ratio):
The characteristic parameter PRM24 is a parameter relating to the appearance timing of the maximum point appearing second when the peak appears second after the minimum point in the reference period with respect to the longitudinal acceleration, that is, when the zero cross point that turns from the minimum point negative to positive is used as a reference. It is. For example, the PRM 24 is a dimensionless amount T XAP2 / StepT obtained by normalizing the timing T XAP2 with the reference period StepT.
第25の特徴パラメータPRM25(X軸加速度最小点タイミング比)について:
特徴パラメータPRM25は、前後軸加速度について基準期間において現れる最小点の出現タイミングに関するパラメータである。例えば、PRM25は、当該タイミングTXAMNを基準期間StepTで規格化することで得られる無次元量TXAMN/StepTである。
Regarding the 25th characteristic parameter PRM25 (X-axis acceleration minimum point timing ratio):
The characteristic parameter PRM25 is a parameter related to the appearance timing of the minimum point that appears in the reference period with respect to the longitudinal acceleration. For example, the PRM 25 is a dimensionless amount T XAMN / StepT obtained by normalizing the timing T XAMN with the reference period StepT.
第26の特徴パラメータPRM26(X軸加速度第1極大点値)について:
特徴パラメータPRM26は、前後軸加速度について基準期間において最小点の後で現れる最初のピーク、即ち、最小点負から正に転じるゼロクロス点を基準としたとき最初に現れる極大点の大きさに関するパラメータである。例えば、PRM26は、当該点の値XAP1(単位は[m/sec2])である。
About the 26th characteristic parameter PRM26 (X-axis acceleration first maximum point value):
The characteristic parameter PRM26 is a parameter relating to the size of the maximum point that appears first when the first peak that appears after the minimum point in the reference period with respect to the longitudinal acceleration, that is, the zero cross point that changes from the minimum point negative to positive, is used as a reference. . For example, the PRM 26 is the value XAP1 (unit: [m / sec 2 ]) of the point.
第27の特徴パラメータPRM27(X軸加速度第2極大点振幅比)について:
特徴パラメータPRM27は、前後軸加速度について基準期間において最小点の後で2番目に現れるピーク、即ち、最小点負から正に転じるゼロクロス点を基準としたとき2番目に現れる極大点の大きさに関するパラメータである。例えば、PRM27は、当該点の値XAP2を1番目の極大点の値で規格化することで得られる無次元量XAP2/XAP1である。
Regarding the 27th characteristic parameter PRM27 (X-axis acceleration second maximum point amplitude ratio):
The characteristic parameter PRM27 is a parameter relating to the size of the maximum point that appears second when the peak appears second after the minimum point in the reference period with respect to the longitudinal acceleration, that is, the zero cross point that turns from the minimum point negative to positive. It is. For example, the PRM 27 is a dimensionless amount XAP2 / XAP1 obtained by normalizing the value XAP2 of the point with the value of the first maximum point.
第28の特徴パラメータPRM28(X軸加速度最小点値)について:
特徴パラメータPRM28は、前後軸加速度について基準期間における最小点の大きさに関するパラメータである。例えば、PRM28は、当該点の値XAMN(単位は[m/sec2])である。
Regarding the twenty-eighth characteristic parameter PRM28 (X-axis acceleration minimum point value):
The characteristic parameter PRM28 is a parameter relating to the size of the minimum point in the reference period for the longitudinal acceleration. For example, the PRM 28 is the value XAMN (unit: [m / sec 2 ]) of the point.
最後に、図11を参照し、第29の特徴パラメータPRM29から第33の特徴パラメータPRM33までの5種類の特徴パラメータについて説明する。 Finally, the five types of feature parameters from the 29th feature parameter PRM29 to the 33rd feature parameter PRM33 will be described with reference to FIG.
第29の特徴パラメータPRM29(X軸速度最大点タイミング比)について:
特徴パラメータPRM29は、前後軸加速度の一階積分(前後軸速度)についての基準期間における最大点の出現タイミングに関するパラメータである。例えば、PRM29は、当該タイミングTXVMXを基準期間StepTで規格化することで得られる無次元量TZVMX/StepTである。
Regarding the 29th characteristic parameter PRM29 (X-axis speed maximum point timing ratio):
The characteristic parameter PRM29 is a parameter related to the appearance timing of the maximum point in the reference period for the first-order integration (front-rear axis velocity) of the front-rear axis acceleration. For example, the PRM 29 is a dimensionless amount T ZVMX / StepT obtained by normalizing the timing T XVMX with the reference period StepT.
第30の特徴パラメータPRM30(X軸速度最大点前90%タイミング比)について:
特徴パラメータPRM30は、前後軸加速度の一階積分(前後軸速度)について基準期間において初めて最大点の90%に到達したタイミングに関するパラメータである。例えば、PRM30は、当該タイミングTXVMX90を基準期間StepTで規格化することで得られる無次元量TXVMX90/StepTである。
About the 30th characteristic parameter PRM30 (90% timing ratio before the maximum X-axis speed):
The characteristic parameter PRM30 is a parameter related to the timing at which 90% of the maximum point is reached for the first time in the reference period with respect to the first order integration (front-rear axis velocity) of the front-rear axis acceleration. For example, the PRM 30 is a dimensionless amount T XVMX90 / StepT obtained by normalizing the timing T XVMX90 with the reference period StepT.
第31の特徴パラメータPRM31(X軸速度最大点前80%タイミング比)について:
特徴パラメータPRM31は、前後軸加速度の一階積分(前後軸速度)について基準期間において初めて最大点の80%に到達したタイミングに関するパラメータである。例えば、PRM31は、当該タイミングTXVMX80を基準期間StepTで規格化することで得られる無次元量TXVMX80/StepTである。
About the 31st characteristic parameter PRM31 (80% timing ratio before the maximum X-axis speed):
The characteristic parameter PRM31 is a parameter related to the timing at which 80% of the maximum point is reached for the first time in the reference period with respect to the first order integral (front-rear axis velocity) of the front-rear axis acceleration. For example, the PRM 31 is a dimensionless quantity T XVMX80 / StepT obtained by normalizing the timing T XVMX80 with the reference period StepT.
第32の特徴パラメータPRM32(X軸速度最小点タイミング比)について:
特徴パラメータPRM32は、前後軸加速度の一階積分(前後軸速度)についての基準期間における最小点の出現タイミングに関するパラメータである。例えば、PRM32は、当該タイミングTXVMNを基準期間StepTで規格化することで得られる無次元量TXVMN/StepTである。
Regarding the thirty-second feature parameter PRM32 (X-axis speed minimum point timing ratio):
The characteristic parameter PRM32 is a parameter related to the appearance timing of the minimum point in the reference period for the first order integral (front-rear axis velocity) of the front-rear axis acceleration. For example, the PRM 32 is a dimensionless amount T XVMN / StepT obtained by normalizing the timing T XVMN with the reference period StepT.
第33の特徴パラメータPRM33(X軸速度最小点前50%タイミング比)について:
特徴パラメータPRM33は、前後軸加速度の一階積分(前後軸速度)について基準期間において初めて最小点の50%に到達したタイミングに関するパラメータである。例えば、PRM33は、当該タイミングTXVMN50を基準期間StepTで規格化することで得られる無次元量TXVMN50/StepTである。
About the thirty-third feature parameter PRM33 (50% timing ratio before the X-axis speed minimum point):
The characteristic parameter PRM33 is a parameter related to the timing at which 50% of the minimum point is reached for the first time in the reference period with respect to the first order integration (front-rear axis velocity) of the front-rear axis acceleration. For example, the PRM 33 is a dimensionless amount T XVMN50 / StepT obtained by normalizing the timing T XVMN50 with the reference period StepT.
以上が本実施の形態において制御部110が導出する33種類の特徴パラメータである。これらの特徴パラメータPRM1~PRM33は、後述するようにステップS46(図13)において、重み付け加算される。本実施の形態においては、次式で表される特徴パラメータPRM1~PRM33についての重み付け加算、
図13に戻り、特徴パラメータ導出後の処理について説明する。ステップS45aにおいて、制御部110は、上下軸軌道時系列データ(図9)を用いて、被測定者の歩幅の広狭を評価する。制御部110は、例えば、一基準期間における上下軸軌道の最大値と最小値との差(図9におけるZTPP)の大きさに基づいて、被測定者の歩幅を推定する。当該推定は、予め実験的に求めた、上下軸軌道の最大値と最小値との差ZTPPと、歩幅との変換テーブル(変換式)等を用いて行えばよい。
Referring back to FIG. 13, the processing after deriving the characteristic parameters will be described. In step S45a, the
次に、ステップS45bにおいて、制御部110は、上下軸加速度時系列データ(図7)を用いて、被測定者の歩行周期のうち1歩分に相当する基準期間を求める。制御部110は、例えば、負から正へ転じるゼロクロス点と、負から正へ転じるゼロクロス点と、の間の時間間隔、即ち、基準期間StepTを求めればよい。
Next, in step S45b, the
次に、ステップS45cにおいて、制御部110は、上記式1による平均骨盤傾斜角θの算出に用いる重み係数を決定する。図14は、重み係数セットの決定に用いる参照テーブルである。制御部110は、先ず、ステップS45aで評価した被測定者の歩幅の広狭にもとづき、被測定者の歩幅が「小さい」か、「中程度」か、「大きい」かを判定する。次に、制御部110は、上記した歩幅の判定を踏まえ、さらに、基準期間が「短い」か、「中程度」か、「長い」か、を判定することにより、使用する重み係数セットKmn(m:1,2,3、n:1,2,3)を決定する。重み係数セットKmnは、それぞれ、式1で用いられる重み係数(k1,k2,...,k32,k33)および定数項k0を備えたセットである。例えば、重み係数セットK11は、歩幅が比較的小さくかつ基準期間が比較的短いという共通の特徴を備えた数多くの被測定者からなる群について式1により平均骨盤傾斜角θが精度良く算出できるように、実験的に求めた重み係数および定数項の値の集合である。同様、重み係数セットK33は、式1により平均骨盤傾斜角θが精度良く算出できるように、歩幅が比較的大きくかつ基準期間が比較的長いという共通の特徴を備えた数多くの被測定者からなる群について実験的に求めた重み係数および定数項の値の集合である。他の重み係数セットも、それぞれ、同様にして実験的に得た重み係数の値および定数項の値の集合である。
Next, in step S45c, the
なお、ステップS45a~S45cまでの処理は、被測定者の歩幅および/または基準期間に応じて複数の重み係数セットを備える場合に式1の評価に用いる重み係数を選択するために行う処理である。したがって、制御部110が用いる重み係数セットが1つである場合には、ステップS45a~S45cまでの処理はスキップされてよい。
Note that the processing from step S45a to S45c is processing performed to select a weighting factor used in the evaluation of
ステップS46において、制御部110は、ステップS44で算出した特徴パラメータPRM1~PRM33と、ステップS45cにおいて決定した重み係数セットKmnとを用いて、次式
図12に戻り、ステップS5において、制御部110は、ステップS4で求めた平均骨盤傾斜角θにもとづき、被測定者の歩行中の姿勢が、「前傾」であるか、「真中」であるか、「後傾」であるか、判定する(図5(A)参照。)。
Returning to FIG. 12, in step S5, the
最後に、制御部110は、ステップS5で得た評価結果をスマートフォン200へ出力(送信)する(ステップS6)。
Finally, the
スマートフォン200は、活動量計100からの情報を受信すると、合計点数とともに評価の結果を表示部240に表示する。スマートフォン200の表示部240には、例えば「あなたの歩行姿勢は、前傾傾向があります。」というようにメッセージとして表示される。なお、表示部240には、イラストやアニメーション等を用いて前傾/後傾の傾向が直感的にわかるような表示を行ってもよい。また、表示部240には、ステップS46で得た平均骨盤傾斜角θそのもの、あるいは、θの値を適切に処理して得られる評価数値を表示してもよい。
When the
この表示部240の表示内容を見て、ユーザは、歩行時の姿勢が前傾ぎみであるか、後傾ぎみであるか否かを知ることができる。
Referring to the display content of the
このように、発明者は、鋭意研究の結果、被験者の腰部に装着された加速度センサから出力される加速度の時間変化の波形の形状の特徴を捉えた1つまたは複数の特徴量(特徴パラメータ)の重み付け加算の結果が、被験者の歩行時における骨盤の傾斜角とよい相関性を有することを実験的に見いだした。この歩行姿勢計1によれば、被測定者の歩行中の姿勢について前傾/後傾であるか否かを適切に評価できる。また、この歩行姿勢計1では、加速度センサ112の出力に基づいて評価を行っているので、モーションキャプチャのような大がかりな設備によらず、手軽に評価できる。
Thus, as a result of earnest research, the inventor has obtained one or a plurality of feature quantities (feature parameters) that capture the characteristics of the waveform shape of the time-varying acceleration output from the acceleration sensor attached to the waist of the subject. It was experimentally found that the result of the weighted addition of has a good correlation with the inclination angle of the pelvis when the subject walks. According to this walking
上述の実施形態では、33種類の特徴パラメータPRM1~PRM33を算出したが、これに限られるものではない。例えば、33種類の特徴パラメータPRM1~PRM33全てを算出するのではなく、一部のみを算出し、その特徴パラメータのみを用いて、被測定者の歩行中の前傾/後傾の傾向を評価してもよい。 In the above-described embodiment, 33 types of characteristic parameters PRM1 to PRM33 are calculated, but the present invention is not limited to this. For example, instead of calculating all 33 types of feature parameters PRM1 to PRM33, only a part is calculated, and only the feature parameters are used to evaluate the tendency of the subject to lean forward / backward during walking. May be.
上述の実施形態では、被測定者の歩幅および/または基準期間に応じてそれぞれ重み係数セットKmnを設定したが、これに限られるものではない。基準期間に代えて歩行周期(左足基準期間と右足基準期間とを併せた2歩分の期間に相当する)を用い、被測定者の歩幅および/または歩行周期に応じてそれぞれ重み係数セットKmnを設定してもよい。 In the above-described embodiment, the weight coefficient set K mn is set in accordance with the step length and / or the reference period of the measurement subject, but the present invention is not limited to this. Instead of the reference period, a walking cycle (corresponding to a period of two steps including the left foot reference period and the right foot reference period) is used, and the weighting coefficient sets K mn are respectively set according to the step length and / or the walking cycle of the measurement subject. May be set.
上述の実施形態では、活動量計100とスマートフォン200とは、BLE通信によって互いに通信を行ったが、これに限られるものではない。例えば、活動量計100とスマートフォン200とは、NFC(Near Field Communication;近距離無線通信)によって、スマートフォン200と活動量計100とが互いに接近したときに通信を行うようにしてもよい。
In the above-described embodiment, the
また、上述の実施形態では、本発明の歩行姿勢計を、活動量計100とスマートフォン200とを含むシステムとして構成したが、これに限られるものではない。
In the above-described embodiment, the walking posture meter of the present invention is configured as a system including the
例えば、本発明の歩行姿勢計を、スマートフォン200のみで構成しても良い。その場合、スマートフォン200が加速度センサを含むものとする。また、スマートフォン200のメモリ220には、制御部210に、ヒトの歩行姿勢が正しい姿勢であるか否かを定量的に評価するプログラム、より詳しくは、歩行中の姿勢の前傾/後傾の程度を評価するプログラムをインストールする。これにより、本発明の歩行姿勢計を小型かつコンパクトに構成することができる。
For example, the walking posture meter of the present invention may be configured with only the
また、そのプログラムは、アプリケーションソフトウェアとして、CD、DVD、フラッシュメモリなどの記録媒体に記録することができる。この記録媒体に記録されたアプリケーションソフトウェアを、スマートフォン、パーソナルコンピュータ、PDA(パーソナル・デジタル・アシスタンツ)などの実質的なコンピュータ装置にインストールすることによって、それらのコンピュータ装置に、ヒトの歩行姿勢が正しい姿勢であるか否かを定量的に評価する方法を実行させることができる。 The program can be recorded as application software on a recording medium such as a CD, a DVD, or a flash memory. The application software recorded on this recording medium is installed in a substantial computer device such as a smartphone, personal computer, PDA (Personal Digital Assistance), etc., so that the human walking posture is correct in these computer devices. It is possible to execute a method for quantitatively evaluating whether or not.
1 歩行姿勢計
100 活動量計
112 加速度センサ
110 制御部CPU
120 メモリ
180 BLE通信部
200 スマートフォン
210 制御部
220 メモリ
290 ネットワーク通信部
230 操作部
240 表示部
280 BLE通信部
1 walking
DESCRIPTION OF
Claims (18)
被測定者の腰の正中線上に装着される加速度センサと、
前記加速度センサが出力する上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、前記被測定者の歩行中の姿勢の前後方向の傾きに対応する量を算出する演算部と、を備えた歩行姿勢計。 A walking posture meter that evaluates a walking posture of a measurement subject,
An acceleration sensor mounted on the midline of the waist of the person being measured;
Using the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output from the acceleration sensor, an amount corresponding to the forward / backward inclination of the posture of the person being measured is calculated. A walking posture meter, comprising:
前記演算部は、前記上下軸加速度の時間変化および前記前後軸加速度の時間変化の特徴を示す特徴パラメータを導出し、当該特徴パラメータに基づいて前記前後方向の傾きに対応する量を算出する、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 1,
The arithmetic unit derives a characteristic parameter indicating characteristics of the temporal change of the vertical axis acceleration and the temporal change of the longitudinal axis acceleration, and calculates an amount corresponding to the tilt in the front-rear direction based on the characteristic parameter; Walking posture meter characterized by
前記演算部は、複数種類の特徴パラメータを導出し、導出した複数種類の特徴パラメータを重み付け加算して得られる量に基づいて前記前後方向の傾きに対応する量を算出する、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 2,
The operation unit derives a plurality of types of feature parameters, and calculates an amount corresponding to the tilt in the front-rear direction based on an amount obtained by weighted addition of the derived types of feature parameters. Walking posture meter.
前記演算部は、前記上下軸加速度の時間的な変化波形に基づいて前記被測定者の歩幅を推定し、当該推定された歩幅に基づいて、前記重み付け加算に用いる重み係数を変化させる、ことを特徴とする歩行姿勢計。 In the walking posture meter according to claim 3,
The computing unit estimates a step length of the person to be measured based on a temporal change waveform of the vertical axis acceleration, and changes a weighting factor used for the weighted addition based on the estimated step length. A characteristic walking posture meter.
前記演算部は、前記上下軸加速度の時間的な変化波形および前記前後軸加速度の時間的な変化波形の少なくともいずれか一方に基づいて前記被測定者の歩行周期またはその歩行周期のうち1歩分に相当する基準期間を求め、当該求めた歩行周期または基準期間に基づいて、前記重み付け加算に用いる重み係数を変化させる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 4,
The calculation unit is configured to calculate one step of the walking cycle of the measurement subject or the walking cycle based on at least one of the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration. And a weighting factor used for the weighted addition is changed based on the obtained walking cycle or reference period.
前記前後方向の傾きに対応する量は、歩行中の前記被測定者の骨盤と水平とがなす角度である、ことを特徴とする歩行姿勢計。 In the walking posture meter according to any one of claims 1 to 5,
An amount corresponding to the tilt in the front-rear direction is an angle formed by the pelvis of the person being measured and a horizontal position during walking.
前記演算部は、上方を正とする前記上下軸加速度の時間的な変化波形において加速度値が負から正へ変化するゼロクロス点の出現タイミングから次の負から正へ変化するゼロクロス点の出現タイミングまでの期間を、前記被測定者の1歩分に相当する基準期間として求め、前記基準期間において前記上下軸加速度の時間的な変化波形に現れる3つの極大点の値および出現タイミング、前記3つの極大点の間の2つの極小点の値および出現タイミング、ならびに、最小点の値および出現タイミングのうちの少なくともいずれか1つを用いて前記特徴パラメータを導出する、ことを特徴とする歩行姿勢計。 In the walking posture meter according to any one of claims 2 to 5,
In the temporal change waveform of the vertical axis acceleration with the upper side being positive, the arithmetic unit is from the appearance timing of the zero cross point at which the acceleration value changes from negative to positive to the next occurrence timing of the zero cross point at which the negative value changes from positive to positive. Is determined as a reference period corresponding to one step of the person to be measured, and the values and timings of the three local maximum points appearing in the temporal change waveform of the vertical axis acceleration in the reference period, and the three local maximums A walking posture meter, wherein the feature parameter is derived using at least one of a value and an appearance timing of two minimum points between points, and a value and an appearance timing of a minimum point.
前記演算部は、前記基準期間の開始点から前記3つの極大点、前記2つの極小点、および、前記最小点のそれぞれの出現タイミングまでの時間と前記基準期間の長さとの比のうちの少なくともいずれか1つを前記特徴パラメータとして用いる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 7,
The calculation unit may include at least one of a ratio of the time from the start point of the reference period to the appearance timing of the three maximum points, the two minimum points, and the minimum point and the length of the reference period. Any one of them is used as the feature parameter.
前記演算部は、前記3つの極大点のそれぞれの値と前記最小点の値との差の絶対値のうちの少なくともいずれか1つを前記特徴パラメータとして用いる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 7,
The arithmetic unit uses at least one of absolute values of differences between the values of the three maximum points and the value of the minimum point as the feature parameter.
前記演算部は、前記上下軸加速度の時間的な変化波形から上下軸速度の時間的な変化波形を求め、
前記演算部は、上方を正とする前記上下軸加速度の時間的な変化波形において加速度値が負から正へ変化するゼロクロス点の出現タイミングから次の負から正へ変化するゼロクロス点の出現タイミングまでの期間を、前記被測定者の1歩分に相当する基準期間として求め、前記基準期間における前記上下軸速度の時間的な変化波形の最大点の出現タイミング、前記最大点の値の90%の値を有する点の出現タイミング、前記最大点の値の80%の値を有する点の出現タイミング、および、前記最大点の値の50%の値を有する点の出現タイミング、のうちの少なくともいずれか1つを用いて前記特徴パラメータを導出する、ことを特徴とする歩行姿勢計。 In the walking posture meter according to any one of claims 2 to 5,
The calculation unit obtains a temporal change waveform of the vertical axis velocity from a temporal change waveform of the vertical axis acceleration,
In the temporal change waveform of the vertical axis acceleration with the upper side being positive, the arithmetic unit is from the appearance timing of the zero cross point at which the acceleration value changes from negative to positive to the next occurrence timing of the zero cross point at which the negative value changes from positive to positive. Is determined as a reference period corresponding to one step of the person to be measured, and the appearance timing of the maximum point of the temporal change waveform of the vertical axis velocity in the reference period is 90% of the value of the maximum point. At least one of an appearance timing of a point having a value, an appearance timing of a point having a value of 80% of the value of the maximum point, and an appearance timing of a point having a value of 50% of the value of the maximum point A walking posture meter, wherein the feature parameter is derived using one.
前記演算部は、前記基準期間の開始点から前記最大点、前記90%の値を有する点、前記80%の値を有する点、および、前記50%の値を有する点のそれぞれの出現タイミングまでの時間と前記基準期間の長さとの比のうちの少なくともいずれか1つを前記特徴パラメータとして用いる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 10,
The calculation unit is configured to start each point of the reference period, the maximum point, the point having the value of 90%, the point having the value of 80%, and the point having the value of 50%. A walking pose meter using at least one of the ratio of the time of the reference period and the length of the reference period as the feature parameter.
前記演算部は、前記上下軸加速度の時間的な変化波形から上下軸軌道の時間的な変化波形を求め、
前記演算部は、上方を正とする前記上下軸加速度の時間的な変化波形において加速度値が負から正へ変化するゼロクロス点の出現タイミングから次の負から正へ変化するゼロクロス点の出現タイミングまでの期間を、前記被測定者の1歩分に相当する基準期間として求め、前記基準期間における前記上下軸軌道の時間的な変化波形の最大点の値および出現タイミング、最小点の値および出現タイミング、ならびに、前記最大点と前記最小点との差の絶対値、のうちの少なくともいずれか1つを用いて前記特徴パラメータを導出する、ことを特徴とする歩行姿勢計。 In the walking posture meter according to any one of claims 2 to 5,
The calculation unit obtains a temporal change waveform of the vertical axis trajectory from the temporal change waveform of the vertical axis acceleration,
In the temporal change waveform of the vertical axis acceleration with the upper side being positive, the arithmetic unit is from the appearance timing of the zero cross point at which the acceleration value changes from negative to positive to the next occurrence timing of the zero cross point at which the negative value changes from positive to positive. Is determined as a reference period corresponding to one step of the person to be measured, and the maximum point value and appearance timing, minimum point value and appearance timing of the temporal change waveform of the vertical axis trajectory in the reference period And the feature parameter is derived using at least one of an absolute value of a difference between the maximum point and the minimum point.
前記演算部は、前記基準期間の開始点から前記最大点、および、前記最小点のそれぞれの出現タイミングまでの時間と前記基準期間の長さとの比のうちの少なくともいずれか1つを前記特徴パラメータとして用いる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 12,
The computing unit calculates at least one of a ratio of a time from the start point of the reference period to the appearance timing of each of the maximum point and the minimum point and the length of the reference period as the characteristic parameter. It is used as a walking posture meter.
前記演算部は、上方を正とする前記上下軸加速度の時間的な変化波形において加速度値が負から正へ変化するゼロクロス点の出現タイミングから次の負から正へ変化するゼロクロス点の出現タイミングまでの期間を、前記被測定者の1歩分に相当する基準期間として求め、前記基準期間における前後軸加速度の時間的な変化波形の最小点の値および出現タイミング、ならびに、前記最小点のあとに現れる2つの極大点の値および出現タイミングのうちの少なくともいずれか1つを用いて前記特徴パラメータを導出する、ことを特徴とする歩行姿勢計。 In the walking posture meter according to any one of claims 2 to 5,
In the temporal change waveform of the vertical axis acceleration with the upper side being positive, the arithmetic unit is from the appearance timing of the zero cross point at which the acceleration value changes from negative to positive to the next occurrence timing of the zero cross point at which the negative value changes from positive to positive. Is determined as a reference period corresponding to one step of the person to be measured, and the minimum point value and appearance timing of the temporal change waveform of the longitudinal acceleration in the reference period, and after the minimum point A walking pose meter, wherein the feature parameter is derived using at least one of the value of the two maximum points that appear and the appearance timing.
前記演算部は、前記基準期間の開始点から前記最小点、および、前記2つの極大点のそれぞれの出現タイミングまでの時間と前記基準期間の長さとの比のうちの少なくともいずれか1つを前記特徴パラメータとして用いる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 14,
The calculation unit calculates at least one of a ratio of a time from the start point of the reference period to the appearance timing of each of the minimum point and the two maximum points and the length of the reference period. A walking pose meter that is used as a feature parameter.
前記演算部は、前記前後軸加速度の時間的な変化波形から前後軸速度の時間的な変化波形を求め、
前記演算部は、上方を正とする前記上下軸加速度の時間的な変化波形において加速度値が負から正へ変化するゼロクロス点の出現タイミングから次の負から正へ変化するゼロクロス点の出現タイミングまでの期間を、前記被測定者の1歩分に相当する基準期間として求め、前記基準期間における前後軸速度の時間的な変化波形の最大点の出現タイミング、前記最大点の値の90%の値を有する点の出現タイミング、前記最大点の値の80%の値を有する点の出現タイミング、最小点の出現タイミング、および、前記最小点の値の50%の値を有する点の出現タイミング、のうちの少なくともいずれか1つを用いて前記特徴パラメータを導出する、ことを特徴とする歩行姿勢計。 In the walking posture meter according to any one of claims 2 to 5,
The calculation unit obtains a temporal change waveform of the longitudinal speed from the temporal change waveform of the longitudinal acceleration,
In the temporal change waveform of the vertical axis acceleration with the upper side being positive, the arithmetic unit is from the appearance timing of the zero cross point at which the acceleration value changes from negative to positive to the next occurrence timing of the zero cross point at which the negative value changes from positive to positive. Is determined as a reference period corresponding to one step of the measured person, and the appearance timing of the maximum point of the temporal change waveform of the longitudinal speed in the reference period, a value of 90% of the value of the maximum point An appearance timing of a point having a value of 80% of the value of the maximum point, an appearance timing of a point having a value of 80% of the value of the maximum point, and an appearance timing of a point having a value of 50% of the value of the minimum point A walking posture meter, wherein the characteristic parameter is derived using at least one of them.
前記演算部は、前記基準期間の開始点から前記最大点、前記最大点の値の90%の値を有する点、前記最大点の値の80%の値を有する点、前記最小点、および、前記最小点の値の50%の値を有する点のそれぞれの出現タイミングまでの時間と前記基準期間の長さとの比のうちの少なくともいずれか1つを前記特徴パラメータとして用いる、ことを特徴とする歩行姿勢計。 The walking posture meter according to claim 16,
The calculation unit includes the maximum point from the start point of the reference period, a point having a value of 90% of the value of the maximum point, a point having a value of 80% of the value of the maximum point, the minimum point, and At least one of the ratios of the time until the appearance timing of each point having a value of 50% of the value of the minimum point and the length of the reference period is used as the feature parameter. Walking posture meter.
前記方法は、
被測定者の腰の正中線上に装着された加速度センサの出力を取得するステップと、
前記加速度センサが出力する上下軸加速度の時間的な変化波形と前後軸加速度の時間的な変化波形とを用いて、前記被測定者の歩行中の姿勢の前後方向の傾きに対応する量を算出するステップと、を備えたことを特徴とするプログラム。 A program for causing a computer to execute a method for evaluating a walking posture of a measurement subject,
The method
Obtaining an output of an acceleration sensor mounted on the midline of the subject's waist;
Using the temporal change waveform of the vertical axis acceleration and the temporal change waveform of the longitudinal axis acceleration output from the acceleration sensor, an amount corresponding to the forward / backward inclination of the posture of the person being measured is calculated. A program comprising the steps of:
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2013100620A JP2014217692A (en) | 2013-05-10 | 2013-05-10 | Walking posture meter and program |
| JP2013-100620 | 2013-05-10 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2014181602A1 true WO2014181602A1 (en) | 2014-11-13 |
Family
ID=51867088
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2014/058910 Ceased WO2014181602A1 (en) | 2013-05-10 | 2014-03-27 | Walking posture meter and program |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP2014217692A (en) |
| WO (1) | WO2014181602A1 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019111521A1 (en) * | 2017-12-06 | 2019-06-13 | 株式会社 資生堂 | Information processing device and program |
| CN113303789A (en) * | 2021-04-30 | 2021-08-27 | 武汉齐物科技有限公司 | Gait event detection method and device based on acceleration |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6750193B2 (en) * | 2015-06-15 | 2020-09-02 | 花王株式会社 | Walking cycle detection method and detection apparatus |
| JP2020113098A (en) * | 2019-01-15 | 2020-07-27 | 日本電信電話株式会社 | Pedometer measuring apparatus and method, and program |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009106390A (en) * | 2007-10-26 | 2009-05-21 | Panasonic Electric Works Co Ltd | Gait detection support system |
| JP2009106374A (en) * | 2007-10-26 | 2009-05-21 | Panasonic Electric Works Co Ltd | Gait information display system |
| JP2011078534A (en) * | 2009-10-06 | 2011-04-21 | Omron Healthcare Co Ltd | Walking posture determination device |
-
2013
- 2013-05-10 JP JP2013100620A patent/JP2014217692A/en active Pending
-
2014
- 2014-03-27 WO PCT/JP2014/058910 patent/WO2014181602A1/en not_active Ceased
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009106390A (en) * | 2007-10-26 | 2009-05-21 | Panasonic Electric Works Co Ltd | Gait detection support system |
| JP2009106374A (en) * | 2007-10-26 | 2009-05-21 | Panasonic Electric Works Co Ltd | Gait information display system |
| JP2011078534A (en) * | 2009-10-06 | 2011-04-21 | Omron Healthcare Co Ltd | Walking posture determination device |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019111521A1 (en) * | 2017-12-06 | 2019-06-13 | 株式会社 資生堂 | Information processing device and program |
| CN113303789A (en) * | 2021-04-30 | 2021-08-27 | 武汉齐物科技有限公司 | Gait event detection method and device based on acceleration |
| CN113303789B (en) * | 2021-04-30 | 2023-01-10 | 武汉齐物科技有限公司 | Gait event detection method and device based on acceleration |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2014217692A (en) | 2014-11-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6131706B2 (en) | Walking posture meter and program | |
| JP5724237B2 (en) | Walking change judgment device | |
| JP6075198B2 (en) | Walking posture meter and program | |
| CN105688396B (en) | Movable information display system and movable information display methods | |
| JP5504810B2 (en) | Walking posture determination device, control program, and control method | |
| JP6288706B2 (en) | Upper body motion measurement system and upper body motion measurement method | |
| US10898112B2 (en) | Gait posture meter and program | |
| JP2012024275A (en) | Walking posture determination device | |
| JP2010005033A (en) | Walking motion analyzer | |
| JP2011177349A (en) | Body motion detector, and display control method for body motion detector | |
| JP2014217691A5 (en) | ||
| JP2015062654A (en) | Gait estimation device and program, fall risk calculation device and program | |
| US20140142466A1 (en) | Physical motion detecting device and control method for physical motion detecting device | |
| WO2014181603A1 (en) | Walking posture meter and program | |
| WO2014181602A1 (en) | Walking posture meter and program | |
| WO2014181605A1 (en) | Walking posture meter and program | |
| WO2022219905A1 (en) | Measurement device, measurement system, measurement method, and recording medium | |
| JP6586173B2 (en) | Walking state detection device and walking state detection program | |
| US7370509B2 (en) | Endurance calculation device, endurance calculation method, and program | |
| JP2012191962A (en) | Position variation detection device, system including the same, and method for operating the position variation detection device | |
| JP5499500B2 (en) | Stride estimation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14794440 Country of ref document: EP Kind code of ref document: A1 |
|
| DPE1 | Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101) | ||
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 14794440 Country of ref document: EP Kind code of ref document: A1 |