Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7626221B2 - Interpolation device, gait measurement system, interpolation method, and program - Google Patents
[go: Go Back, main page]

JP7626221B2 - Interpolation device, gait measurement system, interpolation method, and program - Google Patents

Interpolation device, gait measurement system, interpolation method, and program Download PDF

Info

Publication number
JP7626221B2
JP7626221B2 JP2023529225A JP2023529225A JP7626221B2 JP 7626221 B2 JP7626221 B2 JP 7626221B2 JP 2023529225 A JP2023529225 A JP 2023529225A JP 2023529225 A JP2023529225 A JP 2023529225A JP 7626221 B2 JP7626221 B2 JP 7626221B2
Authority
JP
Japan
Prior art keywords
data
walking
interpolation
missing
waveform
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2023529225A
Other languages
Japanese (ja)
Other versions
JPWO2022269698A5 (en
JPWO2022269698A1 (en
Inventor
晨暉 黄
シンイ オウ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
NEC Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NEC Corp filed Critical NEC Corp
Publication of JPWO2022269698A1 publication Critical patent/JPWO2022269698A1/ja
Publication of JPWO2022269698A5 publication Critical patent/JPWO2022269698A5/en
Application granted granted Critical
Publication of JP7626221B2 publication Critical patent/JP7626221B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/1036Measuring load distribution, e.g. podologic studies
    • A61B5/1038Measuring plantar pressure during gait
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • A61B5/112Gait analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801Arrangements 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/6802Sensor mounted on worn items
    • A61B5/6804Garments; Clothes
    • A61B5/6807Footwear
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/68Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61HPHYSICAL THERAPY APPARATUS, e.g. DEVICES FOR LOCATING OR STIMULATING REFLEX POINTS IN THE BODY; ARTIFICIAL RESPIRATION; MASSAGE; BATHING DEVICES FOR SPECIAL THERAPEUTIC OR HYGIENIC PURPOSES OR SPECIFIC PARTS OF THE BODY
    • A61H2201/00Characteristics of apparatus not provided for in the preceding codes
    • A61H2201/16Physical interface with patient
    • A61H2201/1602Physical interface with patient kind of interface, e.g. head rest, knee support or lumbar support
    • A61H2201/165Wearable interfaces
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61HPHYSICAL THERAPY APPARATUS, e.g. DEVICES FOR LOCATING OR STIMULATING REFLEX POINTS IN THE BODY; ARTIFICIAL RESPIRATION; MASSAGE; BATHING DEVICES FOR SPECIAL THERAPEUTIC OR HYGIENIC PURPOSES OR SPECIFIC PARTS OF THE BODY
    • A61H2205/00Devices for specific parts of the body
    • A61H2205/12Feet
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V40/00Recognition of biometric, human-related or animal-related patterns in image or video data
    • G06V40/20Movements or behaviour, e.g. gesture recognition
    • G06V40/23Recognition of whole body movements, e.g. for sport training
    • G06V40/25Recognition of walking or running movements, e.g. gait recognition
    • GPHYSICS
    • G08SIGNALLING
    • G08BSIGNALLING SYSTEMS, e.g. PERSONAL CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02Alarms for ensuring the safety of persons
    • G08B21/04Alarms for ensuring the safety of persons responsive to non-activity, e.g. of elderly persons
    • G08B21/0438Sensor means for detecting
    • G08B21/0446Sensor means for detecting worn on the body to detect changes of posture, e.g. a fall, inclination, acceleration, gait
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H40/00ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices
    • G16H40/60ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices
    • G16H40/67ICT specially adapted for the management or administration of healthcare resources or facilities; ICT specially adapted for the management or operation of medical equipment or devices for the operation of medical equipment or devices for remote operation

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Description

本開示は、時系列データの欠損を補間する補間装置等に関する。 The present disclosure relates to an interpolation device that interpolates missing time series data.

体調管理を行うヘルスケアへの関心の高まりから、歩行パターンに含まれる特徴(歩容とも呼ぶ)を計測し、歩容に応じた情報をユーザに提供するサービスが注目されている。例えば、慣性センサを含む計測装置を靴等の履物に実装し、ユーザの歩容を解析する装置が開発されている。このような歩容解析においては、センサによって計測されたデータ(センサデータとも呼ぶ)を計測装置から無線送信する際に、通信障害などの要因によって、データの欠損が発生することがある。 With growing interest in health care for managing physical condition, attention is being paid to services that measure the characteristics contained in walking patterns (also called gait) and provide users with information corresponding to their gait. For example, devices have been developed that analyze a user's gait by mounting a measurement device including an inertial sensor on footwear such as shoes. In such gait analysis, when data measured by a sensor (also called sensor data) is wirelessly transmitted from the measurement device, data loss can occur due to factors such as communication failure.

特許文献1には、足が地面に接地している期間(立脚相)と、足が地面から離れている期間(遊脚相)とで異なる特性を示すドリフトを除去することを目的とした歩容計測システムについて開示されている。特許文献1のシステムは、慣性計測ユニットによって計測される加速度データから少なくとも一つの歩行フェーズを検出するとともに、加速度データを時間積分して速度データを計算する。特許文献1のシステムは、歩行フェーズと速度データに基づいて、それぞれの歩行フェーズに対応する補正量を計算する。特許文献1のシステムは、それぞれの歩行フェーズに対応する速度データから補正量を減じて補正速度データを計算し、算出した補正速度データを時間積分して軌跡データを計算する。 Patent Document 1 discloses a gait measurement system that aims to remove drift that shows different characteristics between the period when the foot is on the ground (stance phase) and the period when the foot is off the ground (swing phase). The system of Patent Document 1 detects at least one walking phase from acceleration data measured by an inertial measurement unit, and calculates velocity data by time-integrating the acceleration data. The system of Patent Document 1 calculates a correction amount corresponding to each walking phase based on the walking phase and the velocity data. The system of Patent Document 1 calculates corrected velocity data by subtracting the correction amount from the velocity data corresponding to each walking phase, and calculates trajectory data by time-integrating the calculated corrected velocity data.

特許文献2には、入力信号の欠落部分を補間する信号補間方式について開示されている。特許文献2の方式では、入力信号の欠落以前の信号の周期を検出し、検出された入力信号を蓄積する。特許文献2の方式では、入力信号の検出に応じて、蓄積された信号を読み出して、欠落部分を補間する。特許文献2の方式では、入力信号の欠落開始時点の極性およびレベルを識別し、欠落直前の入力信号との連続性を保つ極性およびレベルの蓄積波形信号を順次読み出す。特許文献2の方式では、蓄積波形信号のレベルを入力信号のレベルに次第に近似させることで、入力信号の欠落部分を連続性が保たれるように補間する。 Patent Document 2 discloses a signal interpolation method for interpolating missing portions of an input signal. In the method of Patent Document 2, the period of the signal before the input signal is missing is detected and the detected input signal is stored. In the method of Patent Document 2, the stored signal is read out in response to the detection of the input signal to interpolate the missing portion. In the method of Patent Document 2, the polarity and level at the start of the input signal missing are identified, and accumulated waveform signals of polarity and level that maintain continuity with the input signal immediately before the missing portion are sequentially read out. In the method of Patent Document 2, the level of the accumulated waveform signal is gradually made to approximate the level of the input signal, thereby interpolating the missing portion of the input signal so that continuity is maintained.

特許文献3には、ネットワークに接続されたセンサや機器によって計測されるデータの欠損を補間するデータ補間装置について開示されている。特許文献3の装置は、データの時系列モデルと、データの欠損を補間する補間方法とを格納する。補間方法は、時系列モデルにおいて、データの時間変化が共通する範囲として特定されるフェーズに対応する。特許文献3の装置は、外部装置から出力されるデータを外部装置ごとに蓄積し、蓄積されたデータの時間変化パターンがいずれの時系列モデルに該当するかを判定する。特許文献3の装置は、蓄積されたデータの一部に含まれる欠損が、判定された時系列モデルのいずれのフェーズに位置しているかを判定し、判定されたフェーズに応じた補間方法で欠損を補間する。 Patent Document 3 discloses a data interpolation device that interpolates missing data measured by sensors or devices connected to a network. The device of Patent Document 3 stores a time series model of data and an interpolation method for interpolating missing data. The interpolation method corresponds to a phase in the time series model that is specified as a range in which the time changes of data are common. The device of Patent Document 3 accumulates data output from an external device for each external device, and determines which time series model the time change pattern of the accumulated data corresponds to. The device of Patent Document 3 determines which phase of the determined time series model a missing part of the accumulated data is located in, and interpolates the missing part using an interpolation method according to the determined phase.

特許文献4には、生体情報に含まれる欠落部分を補完する生体情報測定装置について開示されている。特許文献4の装置は、脈波センサおよび心電センサによって計測された時系列データ(生体情報)に含まれる欠落部分を補完する。特許文献4の装置は、二つのセンサによって計測された時系列データのうち、一方のセンサによって計測されたデータの欠落部分を、他方のセンサによって計測されたデータに基づいて補完する。特許文献4の装置は、他方のセンサによって計測された時系列データの欠落部分の始点以前のデータのうち、欠落部分と同時刻の区間の時系列データとの相関を示す指標が所定の条件を満たす対応区間を特定する。特許文献4の装置は、一方のセンサによって計測されたデータのうち、対応区間と同時刻の区間の時系列データを用いて欠落部分を補完する。 Patent document 4 discloses a bioinformation measuring device that complements missing parts contained in bioinformation. The device of Patent document 4 complements missing parts contained in time series data (biological information) measured by a pulse wave sensor and an electrocardiogram sensor. The device of Patent document 4 complements missing parts of data measured by one of the two sensors based on data measured by the other sensor. The device of Patent document 4 identifies a corresponding section of data before the start point of the missing part of the time series data measured by the other sensor, where an index indicating a correlation between the missing part and the time series data of the same time section satisfies a predetermined condition. The device of Patent document 4 complements the missing part using time series data of the same time section as the corresponding section of data measured by one sensor.

国際公開第2020/105115号International Publication No. 2020/105115 特開昭61-227432号公報Japanese Patent Application Laid-Open No. 61-227432 特開2014-102779号公報JP 2014-102779 A 特開2010-264168号公報JP 2010-264168 A

特許文献1の手法によれば、立脚相や遊脚相のように、異なる歩行ピリオドにおいて明確に異なる特性を示すドリフトを除去することはできる。しかしながら、特許文献1の手法では、立脚相や遊脚相などの個々の歩行ピリオドにおいて局所的に発生した欠落部分を補間することができなかった。そのため、特許文献1の手法では、欠落部分に含まれる特徴までは復元することができなかった。 The method of Patent Document 1 can remove drifts that show clearly different characteristics in different gait periods, such as the stance phase and swing phase. However, the method of Patent Document 1 cannot fill in missing parts that occur locally in individual gait periods, such as the stance phase and swing phase. Therefore, the method of Patent Document 1 cannot restore the features contained in the missing parts.

特許文献2の手法では、入力信号に欠落があった場合、欠落開始時から欠落終了時にかけて連続性が得られるように、蓄積された信号を繰り返し読み出して、欠落部分を補間する。特許文献2の手法によれば、入力信号の周期性や形状が単純な場合には、欠落部分を補間できる。しかしながら、特許文献2の手法では、センサデータのように規則性や波形が複雑な場合には、欠落開始時から欠落終了時にかけて連続性が得られるように、欠落部分を補間することが難しかった。そのため、特許文献2の手法では、欠落部分に含まれる特徴が失われてしまう可能性があった。In the method of Patent Document 2, when there is a gap in the input signal, the accumulated signal is repeatedly read out to interpolate the missing portion so that continuity is obtained from the start of the gap to the end of the gap. According to the method of Patent Document 2, when the input signal has a simple periodicity or shape, the missing portion can be interpolated. However, with the method of Patent Document 2, when the input signal has a complex regularity or waveform, such as sensor data, it is difficult to interpolate the missing portion so that continuity is obtained from the start of the gap to the end of the gap. Therefore, with the method of Patent Document 2, there is a possibility that features contained in the missing portion may be lost.

特許文献3の手法では、時系列モデルと補間方法を予め格納しておく必要がある。そのため、特許文献3の手法では、時系列モデルや補間方法を特定できない欠損に関しては、データの欠損を補間することができなかった。In the method of Patent Document 3, it is necessary to store the time series model and the interpolation method in advance. Therefore, in the method of Patent Document 3, it is not possible to interpolate missing data for missing data for which the time series model or the interpolation method cannot be specified.

特許文献4の手法では、同時刻に計測された互いに関連し合う二つの時系列データを対応させて、一方の時系列データの欠落部分を他方のデータに基づいて補完する。そのため、特許文献4の手法は、同時刻に計測された時系列データがない限り、欠落部分を補完することができなかった。In the method of Patent Document 4, two related time series data measured at the same time are matched, and missing parts of one time series data are complemented based on the other data. Therefore, the method of Patent Document 4 cannot complement missing parts unless there is time series data measured at the same time.

本開示の目的は、センサデータの時系列データに含まれるデータの欠損区間の特徴を含めて、欠損区間のデータを補間できる補間装置等を提供することにある。 The objective of the present disclosure is to provide an interpolation device, etc., capable of interpolating data in missing sections, including the characteristics of the missing sections of data contained in time series data of sensor data.

本開示の一態様の補間装置は、足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成するとともに、時系列データにおけるデータの欠損区間を特定する歩容情報処理部と、特定された欠損区間の歩行フェーズを計算する欠損情報処理部と、欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、欠損区間の歩行フェーズのデータを用いて、欠損区間を補間する補間データを生成し、生成された補間データを欠損区間に補間する補間部と、を備える。 An interpolation device according to one aspect of the present disclosure includes a gait information processing unit that generates a walking waveform for each walking cycle using time-series data of sensor data relating to foot movement and identifies missing sections of data in the time-series data, a missing information processing unit that calculates the walking phase of the identified missing section, and an interpolation unit that generates interpolation data to interpolate the missing section using walking phase data of the missing section in a walking waveform of a walking cycle different from the walking cycle that includes the missing section, and interpolates the generated interpolation data into the missing section.

本開示の一態様の補間方法においては、コンピュータが、足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成し、時系列データにおけるデータの欠損区間を特定し、特定された欠損区間の歩行フェーズを計算し、欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、欠損区間の歩行フェーズのデータを用いて、欠損区間を補間する補間データを生成し、生成された補間データを欠損区間に補間する。In one embodiment of the interpolation method disclosed herein, a computer generates a gait waveform for each gait cycle using time series data of sensor data related to foot movement, identifies missing sections of data in the time series data, calculates the gait phase of the identified missing section, generates interpolation data that interpolates the missing section using data on the gait phase of the missing section in a gait waveform of a gait cycle different from the gait cycle that includes the missing section, and interpolates the generated interpolation data into the missing section.

本開示の一態様のプログラムは、足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成する処理と、時系列データにおけるデータの欠損区間を特定する処理と、特定された欠損区間の歩行フェーズを計算する処理と、欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、欠損区間の歩行フェーズのデータを用いて、欠損区間を補間する補間データを生成する処理と、生成された補間データを欠損区間に補間する処理と、をコンピュータに実行させる。 A program according to one embodiment of the present disclosure causes a computer to execute the following processes: generating a gait waveform for each gait cycle using time series data of sensor data relating to foot movement; identifying missing sections of data in the time series data; calculating the gait phase of the identified missing section; generating interpolation data that interpolates the missing section using data on the gait phase of the missing section in a gait waveform of a gait cycle different from the gait cycle that includes the missing section; and interpolating the generated interpolation data into the missing section.

本開示によれば、センサデータの時系列データに含まれるデータの欠損区間の特徴を含めて、欠損区間のデータを補間できる補間装置等を提供することが可能になる。 According to the present disclosure, it is possible to provide an interpolation device, etc., that can interpolate data in missing sections, including the characteristics of the missing sections of data contained in time series data of sensor data.

第1の実施形態に係る歩容計測システムの構成の一例を示すブロック図である。1 is a block diagram showing an example of a configuration of a gait measurement system according to a first embodiment. FIG. 第1の実施形態に係る歩容計測システムの計測装置の配置例を示す概念図である。FIG. 1 is a conceptual diagram showing an example of arrangement of measurement devices in a gait measurement system according to a first embodiment. 第1の実施形態に係る歩容計測システムの計測装置に設定される座標系の一例について説明するための概念図である。FIG. 2 is a conceptual diagram for explaining an example of a coordinate system set in a measurement device of the gait measurement system according to the first embodiment. 歩行周期について説明するための概念図である。FIG. 1 is a conceptual diagram for explaining a walking cycle. センサデータの時系列データに発生する欠損について説明するための概念図である。FIG. 11 is a conceptual diagram for explaining a defect that occurs in time-series data of sensor data. センサデータの時系列データに発生した欠損を解析的な手法で補間する一例について説明するための概念図である。FIG. 11 is a conceptual diagram for explaining an example of interpolating missing data that occurs in time-series data of sensor data using an analytical method. 第1の実施形態に係る歩容計測システムの補間装置による欠損区間の補間について説明するための概念図である。4 is a conceptual diagram for explaining interpolation of a missing section by the interpolation device of the gait measurement system according to the first embodiment. FIG. 第1の実施形態に係る歩容計測システムの計測装置の構成の一例を示すブロック図である。1 is a block diagram showing an example of the configuration of a measurement device of a gait measurement system according to a first embodiment. FIG. 第1の実施形態に係る歩容計測システムの補間装置の構成の一例を示すブロック図である。1 is a block diagram showing an example of a configuration of an interpolation device of a gait measurement system according to a first embodiment. FIG. 第1の実施形態に係る歩容計測システムの補間装置による欠損の補間の一例を示す概念図である。FIG. 2 is a conceptual diagram showing an example of interpolation of a defect by the interpolation device of the gait measurement system according to the first embodiment. 第1の実施形態に係る歩容計測システムの補間装置による欠損の補間について説明するための概念図である。4 is a conceptual diagram for explaining interpolation of a defect by the interpolation device of the gait measurement system according to the first embodiment. FIG. 第1の実施形態に係る歩容計測システムの歩容計測装置の構成の一例を示すブロック図である。1 is a block diagram showing an example of a configuration of a gait measurement device of a gait measurement system according to a first embodiment. FIG. 第1の実施形態に係る歩容計測システムの歩容計測装置の補間装置に含まれる歩容情報処理部の動作の一例について説明するためのフローチャートである。5 is a flowchart for illustrating an example of the operation of a gait information processing unit included in the interpolation device of the gait measurement system according to the first embodiment. 第1の実施形態に係る歩容計測システムの歩容計測装置の補間装置に含まれる欠損情報処理部の動作の一例について説明するためのフローチャートである。10 is a flowchart for illustrating an example of the operation of a missing information processing unit included in the interpolation device of the gait measurement system according to the first embodiment. 第1の実施形態に係る歩容計測システムの歩容計測装置の補間装置に含まれる補間部の動作の一例について説明するためのフローチャートである。5 is a flowchart for illustrating an example of the operation of an interpolation unit included in the interpolation device of the gait measurement system according to the first embodiment. 第1の実施形態に係る歩容計測システムの歩容計測装置の補間装置に含まれる補間部の動作の別の一例について説明するためのフローチャートである。10 is a flowchart for illustrating another example of the operation of the interpolation unit included in the interpolation device of the gait measurement system according to the first embodiment. 第2の実施形態に係る歩容計測システムの構成の一例を示すブロック図である。FIG. 11 is a block diagram showing an example of a configuration of a gait measurement system according to a second embodiment. 第2の実施形態に係る歩容計測システムの補間装置の構成の一例を示すブロック図である。FIG. 13 is a block diagram showing an example of the configuration of an interpolation device of a gait measurement system according to a second embodiment. 複数の歩行周期に亘るセンサデータの時系列データの一例である。1 is an example of time-series data of sensor data spanning a plurality of walking cycles. 複数の歩行周期に亘るセンサデータの時系列データから抽出された、一歩行周期分の時系列データを正規化した歩行波形の一例である。13 is an example of a gait waveform obtained by normalizing time series data for one walking cycle, extracted from time series data of sensor data spanning a plurality of walking cycles. 複数の歩行周期の歩行波形をサンプル群(全体)とする、歩行フェーズP1および歩行フェーズP2におけるデータ値の頻度分布である。13 shows a frequency distribution of data values in walking phase P1 and walking phase P2, with walking waveforms of multiple walking cycles forming a sample group (whole). 連続する歩行フェーズのデータNおよびデータN+1のデータ値の全体分布における偏差値の相関関係を示すグラフである。13 is a graph showing the correlation of standard deviations in the overall distribution of data values of data N and data N+1 of consecutive walking phases. 対象歩行フェーズからの距離に応じた偏差値の相関係数を示すグラフである。13 is a graph showing the correlation coefficient of the deviation value according to the distance from a target walking phase. 第2の実施形態に係る補間装置による欠損区間の補間について説明するためのグラフである。10 is a graph for explaining the interpolation of a missing section by an interpolation device according to a second embodiment; 第2の実施形態に係る補間装置の動作について説明するためのフローチャートである。10 is a flowchart illustrating an operation of the interpolation device according to the second embodiment. 第3の実施形態の補間装置の構成の一例を示すブロック図である。FIG. 13 is a block diagram showing an example of a configuration of an interpolation device according to a third embodiment. 各実施形態に係る制御や処理を実行するハードウェア構成の一例を示すブロック図である。FIG. 2 is a block diagram showing an example of a hardware configuration for executing control and processing according to each embodiment.

以下に、本発明を実施するための形態について図面を用いて説明する。ただし、以下に述べる実施形態には、本発明を実施するために技術的に好ましい限定がされているが、発明の範囲を以下に限定するものではない。なお、以下の実施形態の説明に用いる全図においては、特に理由がない限り、同様箇所には同一符号を付す。また、以下の実施形態において、同様の構成・動作に関しては繰り返しの説明を省略する場合がある。 Below, the form for implementing the present invention is explained using the drawings. However, the embodiment described below has technically preferable limitations for implementing the present invention, but does not limit the scope of the invention to the following. In addition, in all the drawings used to explain the following embodiments, the same symbols are used for similar parts unless there is a special reason. Also, in the following embodiments, repeated explanations of similar configurations and operations may be omitted.

(第1の実施形態)
まず、第1の実施形態に係る歩容計測システムについて図面を参照しながら説明する。本実施形態の歩容計測システムは、ユーザの履く履物に設置された計測装置によって、足の動きに関する物理量(センサデータ)を計測する。計測装置は、加速度センサや角速度センサを含む。例えば、足の動きに関する物理量は、加速度センサによって計測される3軸方向の加速度(空間加速度とも呼ぶ)や、角速度センサによって計測される3軸周りの角速度(空間角速度とも呼ぶ)を含む。本実施形態の歩容計測システムは、計測されたセンサデータの通信時などに発生したデータの欠損部分を補間する。
(First embodiment)
First, a gait measurement system according to a first embodiment will be described with reference to the drawings. The gait measurement system of this embodiment measures physical quantities (sensor data) related to foot movement by a measurement device installed in footwear worn by a user. The measurement device includes an acceleration sensor and an angular velocity sensor. For example, the physical quantities related to foot movement include acceleration in three axial directions (also called spatial acceleration) measured by an acceleration sensor and angular velocity around three axes (also called spatial angular velocity) measured by an angular velocity sensor. The gait measurement system of this embodiment interpolates missing parts of data that occur during communication of measured sensor data, etc.

(構成)
図1は、本実施形態の歩容計測システム1の構成を示すブロック図である。歩容計測システム1は、計測装置11、補間装置12、および歩容計測装置13を備える。補間装置12は、計測装置11および歩容計測装置13に有線で接続されてもよいし、無線で接続されてもよい。また、計測装置11、補間装置12、および歩容計測装置13は、単一の装置で構成されてもよい。また、歩容計測システム1は、計測装置11を除き、補間装置12および歩容計測装置13で構成されてもよい。
(composition)
1 is a block diagram showing the configuration of a gait measurement system 1 according to this embodiment. The gait measurement system 1 includes a measurement device 11, an interpolation device 12, and a gait measurement device 13. The interpolation device 12 may be connected to the measurement device 11 and the gait measurement device 13 by wire or wirelessly. The measurement device 11, the interpolation device 12, and the gait measurement device 13 may be configured as a single device. The gait measurement system 1 may be configured from the interpolation device 12 and the gait measurement device 13, excluding the measurement device 11.

計測装置11は、足部に設置される。例えば、計測装置11は、靴等の履物に設置される。例えば、計測装置11は、足弓の裏側の位置に配置される。計測装置11は、加速度センサおよび角速度センサを含む。計測装置11は、履物を履くユーザの足の動きに関する物理量として、加速度センサによって計測される加速度(空間加速度とも呼ぶ)や、角速度センサによって計測される角速度(空間角速度とも呼ぶ)を計測する。計測装置11が計測する足の動きに関する物理量には、加速度や角速度を積分することによって計算される速度や角度、位置(軌跡)も含まれる。計測装置11は、計測された物理量をデジタルデータ(センサデータとも呼ぶ)に変換する。計測装置11は、変換後のセンサデータを補間装置12に送信する。センサデータには、センサデータが取得された時刻に対応するタイムスタンプを含む。タイムスタンプは、センサデータに付与された時系列の番号である。例えば、計測装置11は、ユーザが携帯する携帯端末(図示しない)を介して、補間装置12に接続される。The measuring device 11 is installed on the foot. For example, the measuring device 11 is installed on footwear such as shoes. For example, the measuring device 11 is placed at the back of the arch of the foot. The measuring device 11 includes an acceleration sensor and an angular velocity sensor. The measuring device 11 measures the acceleration (also called spatial acceleration) measured by the acceleration sensor and the angular velocity (also called spatial angular velocity) measured by the angular velocity sensor as physical quantities related to the movement of the foot of the user wearing the footwear. The physical quantities related to the movement of the foot measured by the measuring device 11 also include the velocity, angle, and position (trajectory) calculated by integrating the acceleration and angular velocity. The measuring device 11 converts the measured physical quantities into digital data (also called sensor data). The measuring device 11 transmits the converted sensor data to the interpolation device 12. The sensor data includes a timestamp corresponding to the time when the sensor data was acquired. The timestamp is a time series number assigned to the sensor data. For example, the measuring device 11 is connected to the interpolation device 12 via a mobile terminal (not shown) carried by the user.

携帯端末(図示しない)は、ユーザによって携帯可能な通信機器である。例えば、携帯端末は、スマートフォンやスマートウォッチ、携帯電話等の通信機能を有する携帯型の通信機器である。携帯端末は、ユーザの足の動きに関するセンサデータを計測装置11から受信する。携帯端末は、受信されたセンサデータを、補間装置12が実装されたサーバやクラウド等に送信する。なお、補間装置12や歩容計測装置13の機能は、携帯端末にインストールされたアプリケーションソフトウェア等によって実現されていてもよい。その場合、携帯端末は、受信されたセンサデータを、自身にインストールされたアプリケーションソフトウェア等によって処理する。 The mobile terminal (not shown) is a communication device that can be carried by a user. For example, the mobile terminal is a portable communication device with communication functions such as a smartphone, a smart watch, or a mobile phone. The mobile terminal receives sensor data related to the movement of the user's feet from the measurement device 11. The mobile terminal transmits the received sensor data to a server or cloud in which the interpolation device 12 is implemented. Note that the functions of the interpolation device 12 and the gait measurement device 13 may be realized by application software or the like installed on the mobile terminal. In that case, the mobile terminal processes the received sensor data using application software or the like installed on the mobile terminal.

計測装置11は、例えば、加速度センサと角速度センサを含む慣性計測装置によって実現される。慣性計測装置の一例として、IMU(Inertial Measurement Unit)があげられる。IMUは、3軸方向の加速度を計測する加速度センサと、3軸周りの角速度を計測する角速度センサを含む。また、計測装置11は、VG(Vertical Gyro)やAHRS(Attitude Heading)などの慣性計測装置によって実現されてもよい。また、計測装置11は、GPS/INS(Global Positioning System/Inertial Navigation System)によって実現されてもよい。The measurement device 11 is realized, for example, by an inertial measurement device including an acceleration sensor and an angular velocity sensor. An example of an inertial measurement device is an IMU (Inertial Measurement Unit). The IMU includes an acceleration sensor that measures acceleration in three axial directions and an angular velocity sensor that measures angular velocity around three axes. The measurement device 11 may also be realized by an inertial measurement device such as a VG (Vertical Gyro) or an AHRS (Attitude Heading). The measurement device 11 may also be realized by a GPS/INS (Global Positioning System/Inertial Navigation System).

図2は、計測装置11を靴100の中に配置する一例を示す概念図である。図2の例では、計測装置11は、足弓の裏側に当たる位置に配置される。例えば、計測装置11は、靴100の中に挿入されるインソールに配置される。例えば、計測装置11は、靴100の底面に配置される。例えば、計測装置11は、靴100の本体に埋設されてもよい。計測装置11は、靴100から着脱できてもよいし、靴100から着脱できなくてもよい。なお、計測装置11は、足の動きに関するセンサデータを取得できさえすれば、足弓の裏側ではない位置に配置されてもよい。また、計測装置11は、ユーザが履く靴下や、ユーザが装着するアンクレット等の装飾品に設置されてもよい。また、計測装置11は、足に直に貼り付けられたり、足に埋め込まれたりしてもよい。図2においては、右足側の靴100に計測装置11が配置される例を示すが、両足分の靴100に計測装置11が配置されてもよい。両足分の靴100に計測装置11が配置されれば、両足分の足の動きに基づいて、歩容を計測できる。2 is a conceptual diagram showing an example of the arrangement of the measuring device 11 in the shoe 100. In the example of FIG. 2, the measuring device 11 is arranged at a position that hits the back side of the arch of the foot. For example, the measuring device 11 is arranged in an insole inserted into the shoe 100. For example, the measuring device 11 is arranged on the bottom surface of the shoe 100. For example, the measuring device 11 may be embedded in the body of the shoe 100. The measuring device 11 may be detachable from the shoe 100, or may not be detachable from the shoe 100. Note that the measuring device 11 may be arranged at a position other than the back side of the arch of the foot, as long as it can acquire sensor data related to the movement of the foot. The measuring device 11 may also be installed in socks worn by the user or in an accessory such as an anklet worn by the user. The measuring device 11 may also be attached directly to the foot or embedded in the foot. FIG. 2 shows an example in which the measuring device 11 is arranged in the shoe 100 on the right foot side, but the measuring device 11 may be arranged in the shoes 100 for both feet. If the measuring device 11 is placed in the shoes 100 for both feet, the gait can be measured based on the movements of both feet.

図3は、計測装置11を足弓の裏側に設置する場合に、計測装置11に設定されるローカル座標系(x軸、y軸、z軸)と、地面に対して設定される世界座標系(X軸、Y軸、Z軸)について説明するための概念図である。世界座標系(X軸、Y軸、Z軸)では、ユーザが直立した状態で、ユーザの横方向がX軸方向(右向きが正)、ユーザの正面の方向(進行方向)がY軸方向(前向きが正)、重力方向がZ軸方向(鉛直上向きが正)に設定される。本実施形態においては、計測装置11を基準とするx方向、y方向、およびz方向からなるローカル座標系を設定する。3 is a conceptual diagram for explaining the local coordinate system (x-axis, y-axis, z-axis) set in the measuring device 11 when the measuring device 11 is placed on the back side of the arch of the foot, and the world coordinate system (x-axis, y-axis, z-axis) set with respect to the ground. In the world coordinate system (x-axis, y-axis, z-axis), when the user is standing upright, the user's lateral direction is set to the x-axis direction (rightward is positive), the direction in front of the user (direction of travel) is set to the y-axis direction (forward is positive), and the direction of gravity is set to the z-axis direction (vertically upward is positive). In this embodiment, a local coordinate system consisting of the x-direction, y-direction, and z-direction is set with the measuring device 11 as the reference.

図4は、右足を基準とする一歩行周期について説明するための概念図である。図4は、右足の踵が地面に着地した時点を起点とし、次に右足の踵が地面に着地した時点を終点とする、右足の一歩行周期を表す。図4は、右足の一歩行周期を100%として、正規化された歩行周期である。片足の一歩行周期は、足の裏側の少なくとも一部が地面に接している立脚相と、足の裏側が地面から離れている遊脚相とに大別される。本実施形態においては、立脚相が60%を占め、遊脚相が40%を占めるように、歩行周期を正規化する。立脚相は、さらに、立脚初期T1、立脚中期T2、立脚終期T3、遊脚前期T4に細分される。遊脚相は、さらに、遊脚初期T5、遊脚中期T6、遊脚終期T7に細分される。なお、一歩行周期分の歩行波形は、踵が地面に着地した時点を起点としなくてもよい。例えば、一歩行周期分の歩行波形は、踵が持ち上がる時点を起点および終点としてもよい。 Figure 4 is a conceptual diagram for explaining a step cycle based on the right foot. Figure 4 shows one walking cycle of the right foot, starting from the time when the heel of the right foot lands on the ground and ending from the time when the heel of the right foot lands on the ground. Figure 4 shows a walking cycle normalized by setting one walking cycle of the right foot as 100%. One walking cycle of one foot is roughly divided into a stance phase in which at least a part of the sole of the foot is in contact with the ground, and a swing phase in which the sole of the foot is off the ground. In this embodiment, the walking cycle is normalized so that the stance phase accounts for 60% and the swing phase accounts for 40%. The stance phase is further divided into an early stance phase T1, a mid stance phase T2, a final stance phase T3, and a pre-swing phase T4. The swing phase is further divided into an early swing phase T5, a mid swing phase T6, and a final swing phase T7. The walking waveform for one step cycle does not have to start from the time when the heel hits the ground. For example, the walking waveform for one step cycle may start and end from the time when the heel lifts off the ground.

図4(a)は、右足の踵が接地する事象(踵接地)を表す(HS:Heel Strike)。図4(b)は、右足の足裏の接地面が接地した状態で、左足の爪先が地面から離れる事象(反対足爪先離地)を表す(OTO:Opposite Toe Off)。図4(c)は、右足の足裏の接地面が接地した状態で、右足の踵が持ち上がる事象(踵持ち上がり)を表す(HR:Heel Rise)。図4(d)は、左足の踵が接地した事象(反対足踵接地)である(OHS:Opposite Heel Strike)。図4(e)は、左足の足裏の接地面が接地した状態で、右足の爪先が地面から離れる事象(爪先離地)を表す(TO:Toe Off)。図4(f)は、左足の足裏の接地面が接地した状態で、左足と右足が交差する事象(足交差)を表す(FA:Foot Adjacent)。図4(g)は、左足の足裏が接地した状態で、右足の脛骨が地面に対してほぼ垂直になる事象(脛骨垂直)を表す(TV:Tibia Vertical)。図4(h)は、右足の踵が接地する事象(踵接地)を表す(HS:Heel Strike)。図4(h)は、図4(a)から始まる歩行周期の終点に相当するとともに、次の歩行周期の起点に相当する。 Figure 4(a) shows an event where the heel of the right foot touches the ground (heel strike) (HS: Heel Strike). Figure 4(b) shows an event where the toe of the left foot leaves the ground (opposite toe off) while the sole of the right foot is in contact with the ground (OTO: Opposite Toe Off). Figure 4(c) shows an event where the heel of the right foot rises (heel rise) while the sole of the right foot is in contact with the ground (HR: Heel Rise). Figure 4(d) shows an event where the heel of the left foot touches the ground (opposite heel strike) (OHS: Opposite Heel Strike). Figure 4(e) shows an event where the toe of the right foot leaves the ground (toe off) while the sole of the left foot is in contact with the ground (TO: Toe Off). Fig. 4(f) shows an event where the left and right feet cross (foot crossing) with the sole of the left foot touching the ground (FA: Foot Adjacent). Fig. 4(g) shows an event where the tibia of the right foot is almost perpendicular to the ground (TV: Tibia Vertical) with the sole of the left foot touching the ground. Fig. 4(h) shows an event where the heel of the right foot touches the ground (heel strike) (HS: Heel Strike). Fig. 4(h) corresponds to the end point of the walking cycle that starts from Fig. 4(a) and corresponds to the starting point of the next walking cycle.

立脚初期T1は、踵接地HSから反対足爪先離地OTOまでの期間である。立脚中期T2は、反対足爪先離地OTOから踵持ち上がりHRまでの期間である。立脚終期T3は、踵持ち上がりHRから反対足踵接地OHSまでの期間である。遊脚前期T4は、反対足踵接地OHSから爪先離地TOまでの期間である。遊脚初期T5は、爪先離地TOから足交差FAまでの期間である。遊脚中期T6は、足交差FAから脛骨垂直TVまでの期間である。遊脚終期T7は、脛骨垂直TVから踵接地HSまでの期間である。なお、歩行イベントが発現するタイミングは、人物や身体状態、歩行状態に応じて異なるため、想定される歩行周期と完全に一致するとは限らない。 The initial stance phase T1 is the period from heel-on-striking HS to toe-off OTO. The mid-stance phase T2 is the period from toe-off OTO to heel-lift HR. The final stance phase T3 is the period from heel-lift HR to heel-on-striking OHS. The early swing phase T4 is the period from heel-on-striking OHS to toe-off TO. The initial swing phase T5 is the period from toe-off TO to foot crossing FA. The mid-striking phase T6 is the period from foot crossing FA to tibia vertical TV. The final swing phase T7 is the period from tibia vertical TV to heel-on-striking HS. The timing at which walking events occur varies depending on the person, physical condition, and walking state, and does not necessarily match the expected walking cycle.

歩行動作においては、図4のように、周期的な動作が繰り返される。そのため、異なる歩行周期の同じタイミング(歩行フェーズとも呼ぶ)において使われる骨格と筋肉は、ほぼ同じである。そのため、歩行に伴って計測されるセンサデータの時系列データの波形は、各歩行フェーズにおいて相似する。本実施形態では、欠損の補間対象の歩行周期とは異なる歩行周期の時系列データに含まれる欠損のないデータを用いて、欠損の補間対象の歩行周期の時系列データの欠損区間を補間する。In walking, periodic movements are repeated as shown in FIG. 4. Therefore, the skeleton and muscles used at the same timing (also called walking phase) in different walking cycles are almost the same. Therefore, the waveforms of the time series data of the sensor data measured during walking are similar in each walking phase. In this embodiment, missing sections of the time series data of the walking cycle to be interpolated are interpolated using data without missing parts contained in the time series data of a walking cycle different from the walking cycle to be interpolated.

補間装置12は、計測装置11からセンサデータを受信する。補間装置12は、受信したセンサデータの時系列データを生成する。補間装置12は、センサデータに含まれるタイムスタンプに基づいて、時系列データを生成する。補間装置12は、生成された時系列データから一歩行周期分の波形を切り出す。例えば、補間装置12は、踵接地や踵持ち上がりのタイミングを起点とする一歩行周期分の波形を切り出す。補間装置12は、切り出された一歩行周期分の波形を正規化して、一歩行周期分の波形(歩行波形とも呼ぶ)を生成する。 The interpolation device 12 receives sensor data from the measurement device 11. The interpolation device 12 generates time series data of the received sensor data. The interpolation device 12 generates the time series data based on the timestamps included in the sensor data. The interpolation device 12 extracts a waveform for one step walking cycle from the generated time series data. For example, the interpolation device 12 extracts a waveform for one step walking cycle starting from the timing of heel strike or heel lift. The interpolation device 12 normalizes the extracted waveform for one step walking cycle to generate a waveform for one step walking cycle (also called a walking waveform).

補間装置12は、センサデータの時系列データにおいて、タイムスタンプの番号が抜けている部分を欠損区間として特定する。補間装置12は、欠損区間と同じ歩行フェーズにおいて欠損のない他の歩行周期における歩行波形を用いて、欠損区間のデータを補間するための補間データを生成する。補間装置12は、生成された補間データを用いて、欠損区間を補間する。The interpolation device 12 identifies a portion of the time series data of the sensor data where a time stamp number is missing as a missing section. The interpolation device 12 generates interpolation data for interpolating data of the missing section using a walking waveform in another walking cycle that is free of missing data in the same walking phase as the missing section. The interpolation device 12 uses the generated interpolation data to interpolate the missing section.

図5は、正常波形と、欠損を含む波形(欠損波形とも呼ぶ)とを比較するための概念図である。図5は、踵持ち上がりを起点とする一歩行周期の進行方向加速度(Y方向加速度)の歩行波形である。図5の欠損波形には、二か所の欠損区間が含まれる。図5のように、欠損区間のデータには、歩行の特徴が含まれる可能性がある。 Figure 5 is a conceptual diagram for comparing a normal waveform with a waveform that includes missing parts (also called a missing waveform). Figure 5 is a walking waveform of forward acceleration (Y-direction acceleration) in a step cycle starting from heel lift. The missing waveform in Figure 5 includes two missing sections. As shown in Figure 5, the data in the missing sections may include gait characteristics.

図6は、関連技術における、時系列データの欠損区間の補間の一例について説明するための概念図である。図6は、時系列データに含まれるセンサデータのタイムスタンプを参照して、解析的な手法で欠損区間のデータを補間する例を示す。図6には、欠損区間の始点と終点を線分(一次関数)で接続する線形補間(破線)と、欠損区間の始点と終点を曲線で接続する多項式補間(一点鎖線)とを示す。線形補間では、欠損区間の始点と終点が線分で結ばれる。線形補間では、時系列データに表れる振幅の凹凸を復元することができないため、凹凸を含む区間の特徴を復元することができない。多項式補間では、欠損区間の始点と終点が曲線で結ばれる。多項式補間では、時系列データの凹凸を大まかに復元することはできるものの、凹凸の部分に含まれる細かい特徴を復元することができない。 FIG. 6 is a conceptual diagram for explaining an example of interpolation of missing sections of time series data in the related art. FIG. 6 shows an example of interpolating data of missing sections by an analytical method with reference to the time stamps of sensor data included in the time series data. FIG. 6 shows linear interpolation (dashed line) that connects the start point and end point of the missing section with a line segment (linear function) and polynomial interpolation (dash-dotted line) that connects the start point and end point of the missing section with a curve. In linear interpolation, the start point and end point of the missing section are connected with a line segment. In linear interpolation, the unevenness of the amplitude appearing in the time series data cannot be restored, and therefore the characteristics of the section including the unevenness cannot be restored. In polynomial interpolation, the start point and end point of the missing section are connected with a curve. In polynomial interpolation, the unevenness of the time series data can be roughly restored, but the fine features included in the unevenness cannot be restored.

図7は、本実施形態における、時系列データの欠損区間の復元の一例について説明するための概念図である。図7の例では、歩行周期i~i+3の歩行波形の各々に含まれる欠損区間S1~S5のデータを相互に補間する(iは自然数)。歩行周期iの歩行波形は、欠損区間S1および欠損区間S5を含む。歩行周期i+1の歩行波形は、欠損区間S2および欠損区間S3を含む。歩行周期i+2の歩行波形は、欠損区間S4を含む。歩行周期i+3の歩行波形は、欠損区間を含まない。 Figure 7 is a conceptual diagram for explaining an example of restoration of missing sections of time series data in this embodiment. In the example of Figure 7, data of missing sections S1 to S5 contained in each of the walking waveforms of walking cycles i to i+3 are mutually interpolated (i is a natural number). The walking waveform of walking cycle i includes missing sections S1 and S5. The walking waveform of walking cycle i+1 includes missing sections S2 and S3. The walking waveform of walking cycle i+2 includes missing section S4. The walking waveform of walking cycle i+3 does not include any missing sections.

補間装置12は、歩行周期i+1~3の歩行波形に含まれる欠損区間S1に対応する歩行フェーズのデータを用いて、歩行周期iの歩行波形に含まれる欠損区間S1を補間する。補間装置12は、歩行周期i+1~3の歩行波形に含まれる欠損区間S5に対応する歩行フェーズのデータを用いて、歩行周期iの歩行波形に含まれる欠損区間S5を補間する。補間装置12は、歩行周期i+1の歩行波形に含まれる欠損区間S2を、歩行周期i、i+2~3の歩行波形に含まれる欠損区間S2に対応する歩行フェーズのデータで補間する。補間装置12は、歩行周期i+1の歩行波形に含まれる欠損区間S3を、歩行周期i、i+2~3の歩行波形に含まれる欠損区間S3に対応する歩行フェーズのデータで補間する。補間装置12は、歩行周期i+2の歩行波形に含まれる欠損区間S4を、歩行周期i~i+1、i+3の歩行波形に含まれる欠損区間S4に対応する歩行フェーズのデータで補間する。The interpolation device 12 uses data of the walking phase corresponding to the missing section S1 included in the walking waveform of the walking cycle i+1 to 3 to interpolate the missing section S1 included in the walking waveform of the walking cycle i. The interpolation device 12 uses data of the walking phase corresponding to the missing section S5 included in the walking waveform of the walking cycle i+1 to 3 to interpolate the missing section S5 included in the walking waveform of the walking cycle i. The interpolation device 12 interpolates the missing section S2 included in the walking waveform of the walking cycle i+1 with data of the walking phase corresponding to the missing section S2 included in the walking waveform of the walking cycle i, i+2 to 3. The interpolation device 12 interpolates the missing section S3 included in the walking waveform of the walking cycle i+1 with data of the walking phase corresponding to the missing section S3 included in the walking waveform of the walking cycle i, i+2 to 3. The interpolation device 12 interpolates the missing section S4 included in the walking waveform of the walking cycle i+2 with data of the walking phase corresponding to the missing section S4 included in the walking waveform of the walking cycle i to i+1, i+3.

例えば、補間装置12は、ある時系列データで発生した時系列データの欠損区間に関して、その欠損区間のデータが欠損していない少なくとも一つの歩行波形から、その欠損区間に対応する歩行フェーズのデータを少なくとも一つ選択する。例えば、補間装置12は、選択されたデータを用いて、欠損区間を補間するための補間データを生成する。For example, for a missing section of time series data that occurs in certain time series data, the interpolation device 12 selects at least one piece of walking phase data corresponding to the missing section from at least one walking waveform that does not have missing data for the missing section. For example, the interpolation device 12 uses the selected data to generate interpolation data for interpolating the missing section.

例えば、補間装置12は、ある時系列データで発生した時系列データの欠損区間に関して、その欠損区間のデータが欠損していない複数の歩行波形から、その欠損区間に対応する歩行フェーズのデータを複数選択する。補間装置12は、選択された複数のデータに対して加算平均や加重平均などの平均値を計算し、算出された平均値を用いて、欠損区間を補間するための補間データを生成する。例えば、補間装置12は、選択された複数のデータのうち、欠損区間の始点および終点の値が最も近いデータを用いて、補間データを生成する。For example, for a missing section of time series data occurring in certain time series data, the interpolation device 12 selects multiple pieces of walking phase data corresponding to the missing section from multiple walking waveforms in which data for the missing section is not missing. The interpolation device 12 calculates an average value, such as an arithmetic average or a weighted average, for the multiple selected pieces of data, and generates interpolation data for interpolating the missing section using the calculated average value. For example, the interpolation device 12 generates interpolation data using the data among the multiple selected pieces of data that has the closest values for the start point and end point of the missing section.

例えば、補間装置12は、欠損区間の歩行ピリオドに応じて、欠損区間のデータを補間する手法を選択してもよい。例えば、進行方向の加速度(Y方向加速度)の波形は、立脚相の期間においては変動が小さく、遊脚相の期間においては変動が大きい。そのため、補間装置12は、立脚相において発生した欠損区間に関しては解析的な手法で欠損区間のデータを補間し、遊脚相において発生した欠損区間に関しては他の歩行周期の歩行波形に基づいて欠損区間のデータを補間するように構成されてもよい。 For example, the interpolation device 12 may select a method for interpolating data of a missing section depending on the walking period of the missing section. For example, the waveform of the acceleration in the direction of travel (Y-direction acceleration) fluctuates little during the stance phase and fluctuates greatly during the swing phase. Therefore, the interpolation device 12 may be configured to interpolate data of a missing section occurring during the stance phase using an analytical method, and to interpolate data of a missing section occurring during the swing phase based on the walking waveform of another walking cycle.

例えば、補間装置12は、図示しないサーバ等に実装される。例えば、補間装置12は、アプリケーションサーバによって実現されてもよい。例えば、補間装置12は、携帯端末(図示しない)にインストールされたアプリケーションソフトウェア等によって実現されてもよい。For example, the interpolation device 12 is implemented in a server or the like (not shown). For example, the interpolation device 12 may be realized by an application server. For example, the interpolation device 12 may be realized by application software or the like installed on a mobile terminal (not shown).

歩容計測装置13は、補間装置12から、欠損区間を含まない歩行波形を取得する。歩容計測装置13は、取得した歩行波形を用いて、歩容に関する計測を実行する。例えば、歩容計測装置13は、歩行波形から歩行イベントを検出し、検出された歩行イベントに基づいて、ユーザの歩容を計測する。例えば、歩容計測装置13は、検出された歩行イベントに基づいて、ストライド長やステップ長、歩幅、歩隔、足角、歩行速度などの歩容を計測する。例えば、歩容計測装置13は、計測された歩容に基づいて、回内/回外の度合や、外反母趾の進行度、身体の対称性、身体の柔軟度などの身体的な状態を推定してもよい。例えば、歩容計測装置13は、計測された歩容に基づいて、上肢や下肢、上腕、前腕、上腿、下腿の長さなどの身体パラメータを計測してもよい。歩容計測装置13による計測対象項目には、特に限定を加えない。The gait measurement device 13 acquires a gait waveform that does not include a missing section from the interpolation device 12. The gait measurement device 13 uses the acquired gait waveform to perform measurements related to gait. For example, the gait measurement device 13 detects a gait event from the gait waveform and measures the user's gait based on the detected gait event. For example, the gait measurement device 13 measures gait such as stride length, step length, stride width, foot angle, and walking speed based on the detected gait event. For example, the gait measurement device 13 may estimate physical conditions such as the degree of pronation/supination, the progression of hallux valgus, body symmetry, and body flexibility based on the measured gait. For example, the gait measurement device 13 may measure physical parameters such as the length of the upper limbs, lower limbs, upper arms, forearms, upper legs, and lower legs based on the measured gait. There are no particular limitations on the items to be measured by the gait measurement device 13.

歩容計測装置13は、歩容に関する計測結果(歩容情報とも呼ぶ)を出力する。例えば、歩容計測装置13は、表示装置(図示しない)や携帯端末(図示しない)に歩容情報を出力する。表示装置に出力された歩容情報は、表示装置や携帯端末の画面に表示される。例えば、歩容計測装置13は、外部システム(図示しない)に歩容情報を出力する。歩容計測装置13から出力される歩容情報は、任意の用途に使用できる。歩容計測装置13が歩容情報を出力する通信機能については、特に限定を加えない。 The gait measurement device 13 outputs measurement results related to gait (also referred to as gait information). For example, the gait measurement device 13 outputs the gait information to a display device (not shown) or a mobile terminal (not shown). The gait information output to the display device is displayed on the screen of the display device or the mobile terminal. For example, the gait measurement device 13 outputs the gait information to an external system (not shown). The gait information output from the gait measurement device 13 can be used for any purpose. There are no particular limitations on the communication function through which the gait measurement device 13 outputs the gait information.

例えば、歩容計測装置13は、図示しないサーバ等に実装される。例えば、歩容計測装置13は、アプリケーションサーバによって実現されてもよい。例えば、歩容計測装置13は、携帯端末(図示しない)にインストールされたアプリケーションソフトウェア等によって実現されてもよい。For example, the gait measurement device 13 is implemented in a server or the like (not shown). For example, the gait measurement device 13 may be realized by an application server. For example, the gait measurement device 13 may be realized by application software or the like installed on a mobile terminal (not shown).

〔計測装置〕
次に、計測装置11の詳細構成について図面を参照しながら説明する。図8は、計測装置11の詳細構成の一例を示すブロック図である。計測装置11は、加速度センサ111、角速度センサ112、制御部113、および送信部115を有する。なお、計測装置11は、図示しない電源を含む。
[Measuring equipment]
Next, a detailed configuration of the measurement device 11 will be described with reference to the drawings. Fig. 8 is a block diagram showing an example of a detailed configuration of the measurement device 11. The measurement device 11 has an acceleration sensor 111, an angular velocity sensor 112, a control unit 113, and a transmission unit 115. The measurement device 11 also includes a power supply (not shown).

加速度センサ111は、3軸方向の加速度(空間加速度とも呼ぶ)を計測するセンサである。加速度センサ111は、計測した加速度を制御部113に出力する。例えば、加速度センサ111には、圧電型や、ピエゾ抵抗型、静電容量型等の方式のセンサを用いることができる。なお、加速度センサ111に用いられるセンサは、加速度を計測できれば、その計測方式に限定を加えない。The acceleration sensor 111 is a sensor that measures acceleration in three axial directions (also called spatial acceleration). The acceleration sensor 111 outputs the measured acceleration to the control unit 113. For example, a piezoelectric type, a piezo-resistive type, a capacitance type, or other type of sensor can be used for the acceleration sensor 111. Note that there are no limitations on the measurement method of the sensor used for the acceleration sensor 111 as long as it can measure acceleration.

角速度センサ112は、3軸方向の角速度(空間角速度とも呼ぶ)を計測するセンサである。角速度センサ112は、計測した角速度を制御部113に出力する。例えば、角速度センサ112には、振動型や静電容量型等の方式のセンサを用いることができる。なお、角速度センサ112に用いられるセンサは、角速度を計測できれば、その計測方式に限定を加えない。 Angular velocity sensor 112 is a sensor that measures angular velocity (also called spatial angular velocity) in three axial directions. Angular velocity sensor 112 outputs the measured angular velocity to control unit 113. For example, a vibration type, a capacitance type, or other type of sensor can be used for angular velocity sensor 112. Note that there are no limitations on the measurement method of the sensor used for angular velocity sensor 112 as long as it can measure angular velocity.

制御部113は、加速度センサ111および角速度センサ112の各々から、3軸方向の加速度と3軸周りの角速度を取得する。制御部113は、取得した加速度および角速度をデジタルデータに変換し、変換後のデジタルデータ(センサデータとも呼ぶ)を送信部115に出力する。センサデータには、デジタルデータに変換された加速度データと、デジタルデータに変換された角速度データとが少なくとも含まれる。加速度データは、3軸方向の加速度ベクトルを含む。角速度データは、3軸周りの角速度ベクトルを含む。なお、加速度データおよび角速度データには、それらのデータの取得時間が紐付けられる。また、制御部113は、取得した加速度データおよび角速度データに対して、実装誤差や温度補正、直線性補正などの補正を加えたセンサデータを出力するように構成してもよい。また、制御部113は、取得した加速度データおよび角速度データを用いて、3軸周りの角度データを生成してもよい。The control unit 113 acquires the acceleration in three axial directions and the angular velocity around the three axes from each of the acceleration sensor 111 and the angular velocity sensor 112. The control unit 113 converts the acquired acceleration and angular velocity into digital data and outputs the converted digital data (also called sensor data) to the transmission unit 115. The sensor data includes at least the acceleration data converted into digital data and the angular velocity data converted into digital data. The acceleration data includes an acceleration vector in the three axial directions. The angular velocity data includes an angular velocity vector around the three axes. The acceleration data and the angular velocity data are linked to the acquisition time of the data. The control unit 113 may also be configured to output sensor data that has been corrected for mounting errors, temperature correction, linearity correction, and the like, for the acquired acceleration data and angular velocity data. The control unit 113 may also generate angle data around the three axes using the acquired acceleration data and angular velocity data.

例えば、制御部113は、計測装置11の全体制御やデータ処理を行うマイクロコンピュータまたはマイクロコントローラである。例えば、制御部113は、CPU(Central Processing Unit)やRAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ等を有する。制御部113は、加速度センサ111および角速度センサ112を制御して角速度や加速度を計測する。例えば、制御部113は、計測された角速度および加速度等の物理量(アナログデータ)をAD変換(Analog-to-Digital Conversion)する。制御部113は、変換後のデジタルデータをフラッシュメモリ等の記憶部(図示しない)に記憶させる。なお、加速度センサ111および角速度センサ112によって計測された物理量(アナログデータ)は、加速度センサ111および角速度センサ112の各々においてデジタルデータに変換されてもよい。記憶部に記憶されたデジタルデータは、所定のタイミングで送信部115に出力される。For example, the control unit 113 is a microcomputer or microcontroller that performs overall control of the measuring device 11 and data processing. For example, the control unit 113 has a CPU (Central Processing Unit), a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, etc. The control unit 113 controls the acceleration sensor 111 and the angular velocity sensor 112 to measure angular velocity and acceleration. For example, the control unit 113 performs AD conversion (Analog-to-Digital Conversion) of physical quantities (analog data) such as the measured angular velocity and acceleration. The control unit 113 stores the converted digital data in a storage unit (not shown) such as a flash memory. Note that the physical quantities (analog data) measured by the acceleration sensor 111 and the angular velocity sensor 112 may be converted to digital data in each of the acceleration sensor 111 and the angular velocity sensor 112. The digital data stored in the storage unit is output to the transmission unit 115 at a predetermined timing.

送信部115は、制御部113からセンサデータを取得する。送信部115は、取得したセンサデータを補間装置12に送信する。例えば、送信部115は、ケーブルなどの有線を介して、センサデータを補間装置12に送信する。例えば、送信部115は、無線通信を介して、センサデータを補間装置12に送信する。例えば、送信部115は、Bluetooth(登録商標)やWiFi(登録商標)などの規格に則した無線通信機能(図示しない)を介して、センサデータを補間装置12に送信するように構成される。なお、送信部115の通信機能は、Bluetooth(登録商標)やWiFi(登録商標)以外の規格に則していてもよい。The transmitting unit 115 acquires sensor data from the control unit 113. The transmitting unit 115 transmits the acquired sensor data to the interpolation device 12. For example, the transmitting unit 115 transmits the sensor data to the interpolation device 12 via a wired connection such as a cable. For example, the transmitting unit 115 transmits the sensor data to the interpolation device 12 via wireless communication. For example, the transmitting unit 115 is configured to transmit the sensor data to the interpolation device 12 via a wireless communication function (not shown) conforming to a standard such as Bluetooth (registered trademark) or WiFi (registered trademark). Note that the communication function of the transmitting unit 115 may be conforming to a standard other than Bluetooth (registered trademark) or WiFi (registered trademark).

〔補間装置〕
次に、補間装置12の詳細構成について図面を参照しながら説明する。図9は、補間装置12の詳細構成について説明するためのブロック図である。補間装置12は、受信部121、歩容情報処理部122、欠損情報処理部123、記憶部125、補間部126、および送信部128を有する。
[Interpolation device]
Next, a detailed configuration of the interpolation device 12 will be described with reference to the drawings. Fig. 9 is a block diagram for explaining the detailed configuration of the interpolation device 12. The interpolation device 12 has a receiving unit 121, a gait information processing unit 122, a missing information processing unit 123, a storage unit 125, an interpolation unit 126, and a transmitting unit 128.

受信部121は、計測装置11からセンサデータを受信する。受信部121は、受信したセンサデータを歩容情報処理部122に出力する。例えば、受信部121は、ケーブルなどの有線を介して、センサデータを計測装置11から受信する。例えば、受信部121は、無線通信を介して、センサデータを計測装置11から受信する。例えば、受信部121は、Bluetooth(登録商標)やWiFi(登録商標)などの規格に則した無線通信機能(図示しない)を介して、センサデータを計測装置11から受信するように構成される。なお、受信部121の通信機能は、Bluetooth(登録商標)やWiFi(登録商標)以外の規格に則していてもよい。The receiving unit 121 receives sensor data from the measuring device 11. The receiving unit 121 outputs the received sensor data to the gait information processing unit 122. For example, the receiving unit 121 receives the sensor data from the measuring device 11 via a wired connection such as a cable. For example, the receiving unit 121 receives the sensor data from the measuring device 11 via wireless communication. For example, the receiving unit 121 is configured to receive the sensor data from the measuring device 11 via a wireless communication function (not shown) that complies with standards such as Bluetooth (registered trademark) or WiFi (registered trademark). Note that the communication function of the receiving unit 121 may be in accordance with standards other than Bluetooth (registered trademark) or WiFi (registered trademark).

歩容情報処理部122は、計測装置11が設置された靴100を履いたユーザの歩行に伴って計測されたセンサデータを、受信部121から取得する。歩容情報処理部122は、センサデータに含まれるタイムスタンプに基づいて、取得したセンサデータの時系列データを生成する。例えば、歩容情報処理部122は、取得されたセンサデータの座標系を、ローカル座標系から世界座標系に変換する。The gait information processing unit 122 acquires sensor data measured in association with the walking of a user wearing the shoe 100 in which the measuring device 11 is installed from the receiving unit 121. The gait information processing unit 122 generates time-series data of the acquired sensor data based on the timestamp included in the sensor data. For example, the gait information processing unit 122 converts the coordinate system of the acquired sensor data from the local coordinate system to the world coordinate system.

歩容情報処理部122は、センサデータを用いて、足の動きに関する物理量の時系列データを生成する。例えば、歩容情報処理部122は、空間加速度や空間角速度などの時系列データを生成する。また、歩容情報処理部122は、空間加速度や空間角速度を積分し、空間速度や空間角度(足底角)、空間軌跡などの時系列データを生成する。これらの時系列データが歩行波形に相当する。歩容情報処理部122は、一般的な歩行周期や、ユーザに固有の歩行周期に合わせて設定された所定のタイミングや時間間隔において、時系列データを生成する。歩容情報処理部122が時系列データを生成するタイミングは、任意に設定できる。例えば、歩容情報処理部122は、ユーザの歩行が継続されている期間、時系列データを生成し続けるように構成される。また、歩容情報処理部122は、特定のタイミングにおいて、時系列データを生成するように構成されてもよい。The gait information processing unit 122 uses the sensor data to generate time series data of physical quantities related to foot movement. For example, the gait information processing unit 122 generates time series data such as spatial acceleration and spatial angular velocity. The gait information processing unit 122 also integrates the spatial acceleration and spatial angular velocity to generate time series data such as spatial velocity, spatial angle (sole angle), and spatial trajectory. These time series data correspond to the walking waveform. The gait information processing unit 122 generates time series data at a predetermined timing or time interval set according to a general walking cycle or a walking cycle specific to the user. The timing at which the gait information processing unit 122 generates the time series data can be set arbitrarily. For example, the gait information processing unit 122 is configured to continue generating time series data during the period during which the user continues walking. The gait information processing unit 122 may also be configured to generate time series data at a specific timing.

歩容情報処理部122は、生成された時系列データから一歩行周期分の波形を切り出す。例えば、歩容情報処理部122は、踵接地を起点とする一歩行周期分の波形を切り出す。例えば、歩容情報処理部122は、一歩行周期分の進行方向加速度(Y方向加速度)の歩行波形において、踵接地のタイミングを検出する。一歩行周期分の進行方向加速度の歩行波形において、踵接地のタイミングは、最小ピークが検出されるタイミングと、最小ピークの次に現れる極大ピークが検出されるタイミングとの中点のタイミングである。例えば、歩容情報処理部122は、先行する踵接地のタイミングを起点とし、後続する踵接地のタイミングを終点とする時系列データを、踵接地を起点とする一歩行周期分の波形として切り出す。The gait information processing unit 122 extracts a waveform for one gait cycle from the generated time series data. For example, the gait information processing unit 122 extracts a waveform for one gait cycle starting from heel strike. For example, the gait information processing unit 122 detects the timing of heel strike in a gait waveform of forward acceleration (Y-direction acceleration) for one gait cycle. In a gait waveform of forward acceleration for one gait cycle, the timing of heel strike is the midpoint between the timing at which the minimum peak is detected and the timing at which the maximum peak that appears after the minimum peak is detected. For example, the gait information processing unit 122 extracts time series data starting from the timing of the preceding heel strike and ending at the timing of the succeeding heel strike as a waveform for one gait cycle starting from heel strike.

例えば、歩容情報処理部122は、一歩行周期分の進行方向加速度の歩行波形において、爪先離地のタイミングを検出する。一歩行周期分の進行方向加速度の歩行波形において、爪先離地のタイミングは、最大ピークに含まれる二つの山の間に表れる谷のタイミングである。例えば、歩容情報処理部122は、一歩行周期分の進行方向加速度の歩行波形において、最大ピークに含まれる二つの山の間に表れる谷のタイミングを、爪先離地のタイミングとして検出する。For example, the gait information processing unit 122 detects the timing of toe-off in the walking waveform of the forward acceleration for one walking cycle. In the walking waveform of the forward acceleration for one walking cycle, the timing of toe-off is the timing of a valley that appears between two peaks included in the maximum peak. For example, the gait information processing unit 122 detects the timing of toe-off in the walking waveform of the forward acceleration for one walking cycle as the timing of toe-off.

例えば、歩容情報処理部122は、踵持ち上がりを起点とする一歩行周期分の波形を抽出してもよい。足底角の歩行波形において、爪先が踵よりも上に位置する状態(背屈)を負と定義し、爪先が踵よりも下に位置する状態(底屈)を正と定義する。足底角の歩行波形が極小となる時刻は、立脚相開始のタイミングに相当する。足底角の歩行波形が極大となる時刻は、遊脚相開始のタイミングに相当する。立脚相開始の時刻と遊脚相開始の時刻との中点の時刻が、立脚相の中央のタイミングに相当する。立脚相の中央のタイミングは、踵持ち上がりのタイミングに相当する。例えば、歩容情報処理部122は、二歩行分の時系列データに含まれる二つの踵持ち上がりのタイミングのうち、先行するタイミングを起点とし、後続するタイミングを終点とする時系列データを、切り出す。切り出された時系列データは、踵持ち上がりを起点とする一歩行周期分の歩行波形である。For example, the gait information processing unit 122 may extract a waveform for one step cycle starting from a heel lift. In the gait waveform of the plantar angle, a state in which the toe is located above the heel (dorsiflexion) is defined as negative, and a state in which the toe is located below the heel (plantar flexion) is defined as positive. The time at which the gait waveform of the plantar angle becomes minimum corresponds to the timing of the start of the stance phase. The time at which the gait waveform of the plantar angle becomes maximum corresponds to the timing of the start of the swing phase. The time at the midpoint between the start of the stance phase and the start of the swing phase corresponds to the timing of the center of the stance phase. The timing of the center of the stance phase corresponds to the timing of the heel lift. For example, the gait information processing unit 122 extracts time series data of two heel lifts included in the time series data for two steps, starting from the earlier timing and ending from the later timing. The extracted time series data is a gait waveform for one stride cycle starting from heel lift.

歩容情報処理部122は、切り出された一歩行周期分の波形を正規化して、一歩行周期分の波形(歩行波形とも呼ぶ)を生成する。歩容情報処理部122は、一歩行周期分の波形の時間軸を、起点のタイミングが0であり、終点のタイミングが100である歩行周期に正規化する。正規化された波形の歩行周期は、100等分されてパーセンテージで表現される。歩行周期の各パーセンテージを、歩行フェーズと呼ぶ。例えば、歩容情報処理部122は、立脚相と遊脚相の比率が60:40になるように、歩行波形の歩行周期を正規化する。例えば、歩容情報処理部122は、歩行波形に表れる歩行イベントに基づいて、歩行波形を正規化してもよい。例えば、歩容情報処理部122は、踵持ち上がりのタイミングを起点とする歩行波形に関して、爪先離地のタイミングが30%になり、踵接地のタイミングが70%になるように、その歩行波形を正規化する。なお、歩行イベントが発現するタイミングは、人物や身体状態、歩行状態に応じて異なるため、想定される歩行フェーズと完全に一致するとは限らない。The gait information processing unit 122 normalizes the waveform of the extracted step cycle to generate a waveform of the step cycle (also called a walking waveform). The gait information processing unit 122 normalizes the time axis of the waveform of the step cycle to a walking cycle whose starting point timing is 0 and whose ending point timing is 100. The walking cycle of the normalized waveform is divided into 100 equal parts and expressed as a percentage. Each percentage of the walking cycle is called a walking phase. For example, the gait information processing unit 122 normalizes the walking cycle of the walking waveform so that the ratio of the stance phase to the swing phase is 60:40. For example, the gait information processing unit 122 may normalize the walking waveform based on walking events that appear in the walking waveform. For example, the gait information processing unit 122 normalizes the walking waveform with respect to the walking waveform starting from the timing of heel lift so that the timing of toe off is 30% and the timing of heel contact is 70%. Note that the timing at which a walking event occurs varies depending on the person, their physical condition, and their walking state, and therefore does not necessarily match perfectly with the expected walking phase.

また、歩容情報処理部122は、一歩行周期分の歩行波形の振幅(加速度や角速度の値)を正規化してもよい。例えば、歩容情報処理部122は、一歩行周期分の歩行波形に関して、振幅の変動幅が一定の範囲内になるように、歩行波形の振幅を正規化する。例えば、歩容情報処理部122は、一歩行周期分の歩行波形に関して、振幅の最大値が1、振幅の最大値が-1になるように、歩行波形の振幅を正規化する。歩行波形の振幅を正規化すれば、複数の歩行周期の歩行波形の振幅が揃うため、欠損区間にデータを補間しやすくなる。また、歩容情報処理部122は、一歩行周期分の歩行波形のベースラインを正規化してもよい。例えば、歩容情報処理部122は、ベースラインの傾きが0となるように、一歩行周期分の歩行波形を正規化する。歩行波形のベースラインを正規化すれば、複数の歩行周期の歩行波形のベースラインの傾きが揃う。そのため、歩行波形のベースラインを正規化すれば、欠損区間のデータを補間しやすくなる。 The gait information processing unit 122 may also normalize the amplitude (acceleration or angular velocity value) of the walking waveform for one step gait cycle. For example, the gait information processing unit 122 normalizes the amplitude of the walking waveform for one step gait cycle so that the amplitude fluctuation range is within a certain range. For example, the gait information processing unit 122 normalizes the amplitude of the walking waveform for one step gait cycle so that the maximum amplitude value is 1 and the maximum amplitude value is -1. If the amplitude of the walking waveform is normalized, the amplitudes of the walking waveform for multiple gait cycles are aligned, making it easier to interpolate data in missing sections. The gait information processing unit 122 may also normalize the baseline of the walking waveform for one step gait cycle. For example, the gait information processing unit 122 normalizes the walking waveform for one step gait cycle so that the slope of the baseline is 0. If the baseline of the walking waveform is normalized, the slope of the baseline of the walking waveform for multiple gait cycles is aligned. Therefore, if the baseline of the walking waveform is normalized, it becomes easier to interpolate data in missing sections.

歩容情報処理部122は、正規化された歩行波形を記憶部125に記憶させる。例えば、歩容情報処理部122は、欠損区間のない歩行波形のみを記憶部125に記憶させる。その場合、欠損区間のある歩行波形(欠損波形とも呼ぶ)は、欠損情報処理部123が欠損区間を特定してから記憶部125に記憶させればよい。例えば、欠損波形を含めた全ての歩行波形を記憶部125に記憶させてもよい。その場合、記憶部125に記憶された欠損波形を欠損情報処理部123が取得し、その欠損波形に含まれる欠損区間を特定すればよい。記憶部125に記憶された歩行波形の欠損の無い区間のデータ(正常データとも呼ぶ)は、欠損波形に含まれる欠損区間を補間するために用いられる。The gait information processing unit 122 stores the normalized walking waveform in the storage unit 125. For example, the gait information processing unit 122 stores only the walking waveform without a missing section in the storage unit 125. In that case, the missing information processing unit 123 may identify the missing section and then store the walking waveform with a missing section (also called a missing waveform) in the storage unit 125. For example, all walking waveforms including the missing waveform may be stored in the storage unit 125. In that case, the missing information processing unit 123 may acquire the missing waveform stored in the storage unit 125 and identify the missing section included in the missing waveform. The data of the missing section of the walking waveform stored in the storage unit 125 (also called normal data) is used to interpolate the missing section included in the missing waveform.

欠損情報処理部123は、欠損波形を取得する。欠損情報処理部123は、歩容情報処理部122から欠損波形を取得してもよいし、記憶部125から欠損波形を取得してもよい。欠損情報処理部123は、取得した欠損波形において、タイムスタンプの番号が抜けている部分を欠損区間として特定する。例えば、欠損情報処理部123は、一歩行周期の歩行波形において、欠損区間の起点から何番目のデータが欠損しているかを判別し、欠損している歩行フェーズを計算する。例えば、欠損情報処理部123は、タイムスタンプの番号に基づいて、欠損区間の歩行ピリオドを特定してもよい。The missing information processing unit 123 acquires the missing waveform. The missing information processing unit 123 may acquire the missing waveform from the gait information processing unit 122, or may acquire the missing waveform from the memory unit 125. The missing information processing unit 123 identifies a portion of the acquired missing waveform where a timestamp number is missing as a missing section. For example, the missing information processing unit 123 determines which data is missing from the start point of the missing section in the walking waveform of a step cycle, and calculates the missing walking phase. For example, the missing information processing unit 123 may identify the walking period of the missing section based on the timestamp number.

欠損情報処理部123は、欠損波形に含まれる欠損区間に関する情報(欠損情報とも呼ぶ)を、その欠損波形に対応付けて記憶部125に記憶させる。例えば、欠損情報処理部123は、欠損情報として、欠損区間の起点と終点の歩行フェーズを、その欠損波形に対応付けて記憶部125に記憶させる。例えば、欠損情報処理部123は、欠損情報として、欠損区間の起点と終点の歩行フェーズと、その歩行フェーズにおける振幅の値とを、その欠損波形に対応付けて記憶部125に記憶させてもよい。また、欠損情報処理部123は、補間部126に欠損情報を出力するように、構成されてもよい。The missing information processing unit 123 stores information about the missing section included in the missing waveform (also called missing information) in the storage unit 125 in association with the missing waveform. For example, the missing information processing unit 123 stores the walking phases at the start and end of the missing section in association with the missing waveform as missing information in the storage unit 125. For example, the missing information processing unit 123 may store the walking phases at the start and end of the missing section and the amplitude value at the walking phase in association with the missing waveform as missing information in the storage unit 125. The missing information processing unit 123 may also be configured to output the missing information to the interpolation unit 126.

記憶部125には、歩行波形が記憶される。記憶部125には、欠損区間を含まない歩行波形(正常波形とも呼ぶ)が記憶される。また、記憶部125には、欠損区間を含む歩行波形(欠損波形とも呼ぶ)や、欠損区間が修復された歩行波形(修復波形とも呼ぶ)が記憶される。欠損波形は、欠損区間の補間対象である。欠損波形には、その欠損波形の欠損情報が対応付けられる。正常波形と修復波形は、歩容計測装置13による歩容計測に用いられる。記憶部125に記憶された歩行波形は、所定のタイミングで送信部128によって取得され、送信部128から歩容計測装置13に送信される。 A walking waveform is stored in the memory unit 125. A walking waveform that does not include a missing section (also called a normal waveform) is stored in the memory unit 125. In addition, a walking waveform that includes a missing section (also called a missing waveform) and a walking waveform in which a missing section has been repaired (also called a repaired waveform) are stored in the memory unit 125. A missing waveform is the subject of interpolation of the missing section. Missing information of the missing waveform is associated with the missing waveform. The normal waveform and the repaired waveform are used for gait measurement by the gait measurement device 13. The walking waveform stored in the memory unit 125 is acquired by the transmission unit 128 at a predetermined timing and transmitted from the transmission unit 128 to the gait measurement device 13.

補間部126は、欠損波形を取得する。補間部126は、記憶部125または欠損情報処理部123から欠損波形を取得する。補間部126が取得する欠損波形には、欠損情報処理部123によって特定された欠損区間に関する情報(欠損情報)が含まれる。補間部126は、取得した欠損波形とは異なる歩行周期の歩行波形のデータを用いて、その欠損波形に含まれる欠損区間を補間する。補間部126は、欠損波形とは異なる歩行周期の歩行波形から、その欠損区間に対応する歩行フェーズにおける正常データを抽出する。補間部126は、抽出された補間データを用いて、欠損区間を補間するための補間データを生成する。The interpolation unit 126 acquires a missing waveform. The interpolation unit 126 acquires a missing waveform from the memory unit 125 or the missing information processing unit 123. The missing waveform acquired by the interpolation unit 126 includes information (missing information) on the missing section identified by the missing information processing unit 123. The interpolation unit 126 uses data of a walking waveform with a walking cycle different from the acquired missing waveform to interpolate the missing section included in the missing waveform. The interpolation unit 126 extracts normal data in the walking phase corresponding to the missing section from the walking waveform with a walking cycle different from the missing waveform. The interpolation unit 126 generates interpolation data for interpolating the missing section using the extracted interpolation data.

例えば、欠損波形は、歩行フェーズが10~20%の区間に欠損区間を含むものとする。その場合、補間部126は、記憶部125に記憶された歩行波形のうち、歩行フェーズが10~20%の区間に欠損を含まない歩行波形を少なくとも一つ選択する。補間部126は、選択された少なくとも一つの歩行波形から、歩行フェーズが10~20%の区間の正常データを抽出する。例えば、補間部126は、ある欠損波形に含まれる欠損区間に関して、その欠損区間のデータが欠損していない少なくとも一つの歩行波形から、その欠損区間に対応する歩行フェーズの正常データを少なくとも一つ選択する。補間部126は、選択された正常データを用いて、欠損区間を補間するための補間データを生成する。For example, the missing waveform includes a missing section in the section where the walking phase is 10-20%. In this case, the interpolation unit 126 selects at least one walking waveform that does not include a missing section in the section where the walking phase is 10-20% from among the walking waveforms stored in the memory unit 125. The interpolation unit 126 extracts normal data for the section where the walking phase is 10-20% from the at least one selected walking waveform. For example, for a missing section included in a certain missing waveform, the interpolation unit 126 selects at least one normal data for the walking phase corresponding to the missing section from at least one walking waveform that is not missing data for the missing section. The interpolation unit 126 generates interpolation data for interpolating the missing section using the selected normal data.

例えば、補間部126は、ある欠損波形に含まれる欠損区間に関して、その欠損区間のデータが欠損していない複数の歩行波形から、その欠損区間に対応する歩行フェーズのデータを複数選択する。例えば、補間部126は、欠損区間に含まれる単独の歩行フェーズごとに、補間データを計算する。補間部126は、選択された複数のデータに加算平均や加重平均などの平均値を計算し、算出された平均値を用いて、欠損区間を補間するための補間データを生成する。例えば、補間部126は、選択された複数のデータのうち、欠損区間の始点および終点の値が最も近いデータを用いて、補間データを生成する。For example, for a missing section included in a certain missing waveform, the interpolation unit 126 selects multiple pieces of data for the walking phase corresponding to the missing section from multiple walking waveforms in which data for the missing section is not missing. For example, the interpolation unit 126 calculates interpolation data for each individual walking phase included in the missing section. The interpolation unit 126 calculates an average value, such as an arithmetic average or a weighted average, for the multiple selected pieces of data, and generates interpolation data for interpolating the missing section using the calculated average value. For example, the interpolation unit 126 generates interpolation data using data, of the multiple selected pieces of data, that has values for the start point and end point of the missing section that are closest to each other.

例えば、補間部126は、欠損区間の歩行ピリオドに応じて、欠損区間のデータを補間する手法を選択してもよい。例えば、進行方向の加速度(Y方向加速度)の波形は、立脚相の期間においては変動が小さく、遊脚相の期間においては変動が大きい。すなわち、遊脚相の期間(第1歩行ピリオドとも呼ぶ)は、立脚相の期間(第2歩行ピリオドとも呼ぶ)と比べてデータの変動が激しいため、解析的にデータを補間することが難しい。言い換えると、立脚相の期間(第2歩行ピリオド)は、遊脚相の期間(第1歩行ピリオド)と比べてデータの変動が小さいため、解析的にデータを補間することができる。そのため、補間部126は、立脚相の期間(第2歩行ピリオド)に含まれる欠損区間に関しては解析的な手法で補間してもよい。そして、遊脚相の期間(第1歩行ピリオド)に含まれる欠損区間に関しては、他の歩行周期の歩行波形に基づいて補間すればよい。なお、補間部126は、変動が小さい区間であれば、遊脚相(第1歩行ピリオド)に含まれる欠損区間に対しても、解析的な手法で補間してもよい。

For example, the interpolation unit 126 may select a method for interpolating data of a missing section according to the walking period of the missing section. For example, the waveform of the acceleration in the direction of travel (Y-direction acceleration) fluctuates less during the stance phase and fluctuates more during the swing phase. That is, the swing phase (also called the first walking period) has more drastic data fluctuations than the stance phase (also called the second walking period), making it difficult to analytically interpolate data. In other words, the stance phase (second walking period) has less data fluctuations than the swing phase (first walking period), making it possible to analytically interpolate data. Therefore, the interpolation unit 126 may use an analytical method to interpolate the missing section included in the stance phase (second walking period). Then, the missing section included in the swing phase (first walking period) may be interpolated based on the walking waveform of another walking cycle. Note that the interpolation unit 126 may also perform analytical interpolation for a missing section included in the swing phase (first walking period) if the section has small fluctuations.

図10は、解析的に補間する区間における欠損区間の補間について説明するための概念図である。図10の例では、歩行フェーズが0~20%の区間や、歩行フェーズが80%~100%の区間においては、振幅の変動が小さい。それに対し、歩行フェーズが20~80%の区間においては、振幅の変動が大きい。このような場合、補間部126は、歩行フェーズが0~20%の区間や、歩行フェーズが80%~100%の区間については、解析的な補間を行ってもよい。歩行フェーズが20~80%の区間については、振幅の変動が大きいため、補間部126は、他の歩行波形のデータを用いて補間データを生成した方がよい。 Figure 10 is a conceptual diagram for explaining the interpolation of missing sections in analytically interpolated sections. In the example of Figure 10, the amplitude fluctuation is small in the section where the walking phase is 0-20% and the section where the walking phase is 80%-100%. In contrast, the amplitude fluctuation is large in the section where the walking phase is 20-80%. In such a case, the interpolation unit 126 may perform analytical interpolation for the section where the walking phase is 0-20% and the section where the walking phase is 80%-100%. Since the amplitude fluctuation is large in the section where the walking phase is 20-80%, it is better for the interpolation unit 126 to generate interpolated data using data of other walking waveforms.

補間部126は、生成された補間データを用いて、欠損波形に含まれる欠損区間を補間する。例えば、補間部126は、補間データを欠損区間に挿入することで、欠損波形を補間する。例えば、補間部126は、欠損区間の始点から終点に向けて、補間データに含まれる正常データで欠損区間を順番に補間する。例えば、補間部126は、欠損区間の終点から始点に向けて、補間データに含まれる正常データで欠損区間を順番に補間してもよい。例えば、補間部126は、欠損区間の終点および始点から、補間データに含まれる正常データを順番に補間してもよい。The interpolation unit 126 uses the generated interpolation data to interpolate the missing sections included in the missing waveform. For example, the interpolation unit 126 interpolates the missing waveform by inserting the interpolation data into the missing sections. For example, the interpolation unit 126 sequentially interpolates the missing sections with normal data included in the interpolation data from the start point of the missing section to the end point of the missing section. For example, the interpolation unit 126 may sequentially interpolate the missing sections with normal data included in the interpolation data from the end point of the missing section to the start point of the missing section. For example, the interpolation unit 126 may sequentially interpolate the normal data included in the interpolation data from the end point and the start point of the missing section.

例えば、補間部126は、補間データによって構成される波形を欠損区間に挿入する。例えば、補間部126は、補間データの始点および終点の各々が欠損区間の始点および終点の各々と一致するように、補間データによって構成される波形を挿入する。例えば、補間データの始点および終点の各々が、欠損区間の始点および終点の各々と一致しない場合が想定される。このような場合、補間部126は、補間データのベースラインを回転させ、補間データの始点および終点の各々と、欠損区間の始点および終点の各々とが一致するようにして、補間データを欠損区間に挿入する。例えば、補間部126は、欠損区間の起点および終点の振幅を、他の歩行周期の正常波形に基づいて計算する。例えば、補間部126は、算出された欠損区間の起点および終点の振幅と、その欠損区間の実際の起点および終点の振幅との差に応じて、ベースラインを調整する。For example, the interpolation unit 126 inserts a waveform formed by the interpolation data into the missing section. For example, the interpolation unit 126 inserts a waveform formed by the interpolation data so that the start point and the end point of the interpolation data match the start point and the end point of the missing section. For example, it is assumed that the start point and the end point of the interpolation data do not match the start point and the end point of the missing section. In such a case, the interpolation unit 126 rotates the baseline of the interpolation data so that the start point and the end point of the interpolation data match the start point and the end point of the missing section, and inserts the interpolation data into the missing section. For example, the interpolation unit 126 calculates the amplitude of the start point and the end point of the missing section based on the normal waveform of another walking cycle. For example, the interpolation unit 126 adjusts the baseline according to the difference between the calculated amplitude of the start point and the end point of the missing section and the actual amplitude of the start point and the end point of the missing section.

図11は、補間データによって構成される波形を、欠損区間に挿入する一例を示す概念図である。(1)補間部126は、欠損波形に含まれる欠損区間の補間データを生成する。(2)補間部126は、生成された補間データを、欠損波形の欠損区間に挿入する。(3)補間部126は、欠損波形の欠損区間に挿入された補間データにベースラインを設定する。(4)補間部126は、欠損波形の欠損区間に挿入された補間データのベースラインを回転させて、欠損区間の起点および終点の各々に、補間データの起点および終点の各々が一致するように起点終点補正を行う。 Figure 11 is a conceptual diagram showing an example of inserting a waveform composed of interpolated data into a missing section. (1) The interpolation unit 126 generates interpolated data for the missing section included in the missing waveform. (2) The interpolation unit 126 inserts the generated interpolated data into the missing section of the missing waveform. (3) The interpolation unit 126 sets a baseline for the interpolated data inserted into the missing section of the missing waveform. (4) The interpolation unit 126 rotates the baseline of the interpolated data inserted into the missing section of the missing waveform, and performs start-end point correction so that the start and end points of the missing section match the start and end points of the interpolated data, respectively.

補間部126は、欠損区間が補間された歩行波形(修復波形)を記憶部125に記憶させる。例えば、補間部126は、修復波形において、修復された欠損区間を特定できるように、欠損区間に含まれる歩行フェーズにラベルを付してもよい。例えば、補間部126は、欠損区間に仮の補間データが挿入された修復波形(仮修復波形とも呼ぶ)を記憶部125に記憶させてもよい。例えば、歩容計測装置13が、仮修復波形に含まれる欠損区間に挿入された補間データを補正してもよい。記憶部125に格納された修復波形は、欠損を含まない他の歩行波形と同様に、歩容計測装置13による歩容計測に用いられる。The interpolation unit 126 stores the walking waveform (repaired waveform) in which the missing section has been interpolated in the storage unit 125. For example, the interpolation unit 126 may label the walking phase included in the missing section in the repaired waveform so that the repaired missing section can be identified. For example, the interpolation unit 126 may store the repaired waveform (also called the provisional repaired waveform) in which provisional interpolation data has been inserted in the missing section in the storage unit 125. For example, the gait measurement device 13 may correct the interpolation data inserted in the missing section included in the provisional repaired waveform. The repaired waveform stored in the storage unit 125 is used for gait measurement by the gait measurement device 13 in the same way as other walking waveforms that do not include a loss.

送信部128は、記憶部125から歩行波形を取得する。送信部128は、正常波形および修復波形を含む歩行波形を記憶部125から取得する。送信部128は、取得した歩行波形を歩容計測装置13に送信する。例えば、送信部128は、ケーブルなどの有線を介して、歩容計測装置13に歩行波形を送信する。例えば、送信部128は、無線通信を介して、歩容計測装置13に歩行波形を送信する。送信部128と歩容計測装置13の間の通信方式に関しては、特に限定しない。The transmission unit 128 acquires a gait waveform from the memory unit 125. The transmission unit 128 acquires a gait waveform including a normal waveform and a restored waveform from the memory unit 125. The transmission unit 128 transmits the acquired gait waveform to the gait measurement device 13. For example, the transmission unit 128 transmits the gait waveform to the gait measurement device 13 via a wire such as a cable. For example, the transmission unit 128 transmits the gait waveform to the gait measurement device 13 via wireless communication. There are no particular limitations on the communication method between the transmission unit 128 and the gait measurement device 13.

〔歩容計測装置〕
次に、歩容計測システム1が備える歩容計測装置13の詳細構成について図面を参照しながら説明する。図12は、歩容計測装置13の詳細構成の一例を示すブロック図である。歩容計測装置13は、取得部131、検出部132、および歩容計測部133を有する。実際には、歩容計測部133による推測結果を出力する出力部(通信インターフェース)が設けられるが、図12においては出力部の機能が歩容計測部133に含まれるものとする。
[Gait measurement device]
Next, a detailed configuration of gait measurement device 13 provided in gait measurement system 1 will be described with reference to the drawings. Fig. 12 is a block diagram showing an example of a detailed configuration of gait measurement device 13. Gait measurement device 13 has an acquisition unit 131, a detection unit 132, and a gait measurement unit 133. In reality, an output unit (communication interface) that outputs the estimation result by gait measurement unit 133 is provided, but in Fig. 12, the function of the output unit is included in gait measurement unit 133.

取得部131は、補間装置12から歩行波形を取得する。取得部131は、正常波形および修復波形を含む歩行波形を補間装置12から取得する。取得部131は、取得した歩行波形を検出部132に出力する。例えば、取得部131は、ケーブルなどの有線を介して、補間装置12から歩行波形を受信する。例えば、取得部131は、無線通信を介して、補間装置12から歩行波形を受信する。補間装置12と取得部131の間の通信方式に関しては、特に限定しない。The acquisition unit 131 acquires a walking waveform from the interpolation device 12. The acquisition unit 131 acquires a walking waveform including a normal waveform and a restored waveform from the interpolation device 12. The acquisition unit 131 outputs the acquired walking waveform to the detection unit 132. For example, the acquisition unit 131 receives the walking waveform from the interpolation device 12 via a wire such as a cable. For example, the acquisition unit 131 receives the walking waveform from the interpolation device 12 via wireless communication. There are no particular limitations on the communication method between the interpolation device 12 and the acquisition unit 131.

検出部132は、取得部131から歩行波形を取得する。検出部132は、取得した歩行波形から、歩行ピリオドや歩行イベントを検出する。例えば、検出部132は、立脚相や遊脚相などの歩行ピリオドを歩行波形から検出する。例えば、検出部132は、踵接地や、反対足爪先離地、踵持ち上がり、反対足踵接地、爪先離地、足交差、脛骨垂直などの歩行イベントを、歩行波形から検出する。例えば、検出部132は、検出された歩行イベントに基づいて、立脚初期や、立脚中期、立脚終期、遊脚前期、遊脚初期や、遊脚中期、遊脚終期などの歩行ピリオドを、歩行波形から検出する。The detection unit 132 acquires a walking waveform from the acquisition unit 131. The detection unit 132 detects walking periods and walking events from the acquired walking waveform. For example, the detection unit 132 detects walking periods such as a stance phase and a swing phase from the walking waveform. For example, the detection unit 132 detects walking events such as heel strike, opposite foot toe off, heel lift, opposite foot heel strike, toe off, foot crossing, and tibia vertical from the walking waveform. For example, the detection unit 132 detects walking periods such as early stance, mid stance, end stance, early swing, early swing, mid swing, and end swing from the walking waveform based on the detected walking events.

例えば、検出部132は、歩行フェーズのパーセンテージに基づいて、歩行ピリオドや歩行イベントを歩行波形から検出する。同じ人物が同じ条件で歩行していれば、歩行ピリオドや歩行イベントが検出されるタイミングは、ほぼ一定である。人物が異なっていても、歩行ピリオドや歩行イベントが検出されるタイミングは、同じ傾向を示す。そのため、歩行フェーズのパーセンテージに基づいて、歩行ピリオドや歩行イベントを特定できる。For example, the detection unit 132 detects walking periods and walking events from the walking waveform based on the percentage of the walking phase. If the same person is walking under the same conditions, the timing at which walking periods and walking events are detected is almost constant. Even if the person is different, the timing at which walking periods and walking events are detected shows the same tendency. Therefore, walking periods and walking events can be identified based on the percentage of the walking phase.

例えば、検出部132は、歩行波形から抽出される特徴に基づいて、歩行ピリオドや歩行イベントを検出してもよい。例えば、検出部132は、足底角の歩行波形が極小となるタイミングを、立脚相開始のタイミングとして検出する。例えば、検出部132は、足底角の歩行波形が極大となるタイミングを、遊脚相開始のタイミングとして検出する。例えば、検出部132は、立脚相開始と遊脚相開始の中点のタイミングを、踵持ち上がりのタイミングとして検出する。For example, the detection unit 132 may detect walking periods and walking events based on features extracted from the walking waveform. For example, the detection unit 132 detects the timing at which the walking waveform of the plantar angle becomes minimum as the timing of the start of the stance phase. For example, the detection unit 132 detects the timing at which the walking waveform of the plantar angle becomes maximum as the timing of the start of the swing phase. For example, the detection unit 132 detects the timing of the midpoint between the start of the stance phase and the start of the swing phase as the timing of heel lift.

例えば、検出部132は、踵持ち上がりのタイミングを起点とする歩行波形に基づいて、歩行イベントを検出する。検出部132は、進行方向加速度(Y方向加速度)の歩行波形において、最大ピークに含まれる二つの山の間に谷が表れるタイミングを、爪先離地のタイミングとして検出する。例えば、検出部132は、進行方向加速度(Y方向加速度)の歩行波形において、歩行周期が60%を超えたあたりの最小ピークを、遊脚終期における足の急減速のタイミングとして検出する。例えば、検出部132は、進行方向加速度(Y方向加速度)の歩行波形において、歩行周期が70%のあたりの極大ピークを、ヒールロッカーのタイミングとして検出する。例えば、検出部132は、進行方向加速度(Y方向加速度)の歩行波形において、極小ピークと極大ピークの中点のタイミングを、踵接地のタイミングとして検出する。例えば、検出部132は、垂直方向加速度(Z方向加速度)の歩行波形において、爪先離地と踵接地の間の最大ピークのタイミングを、脛骨垂直のタイミングとして検出する。例えば、検出部132は、進行方向加速度(Y方向加速度)の歩行波形において、爪先離地と脛骨垂直のタイミングの間のなだらかな極大ピークのタイミングを、足交差のタイミングとして検出する。例えば、検出部132は、ロール角速度の歩行波形から、反対足踵接地および反対足爪先離地のタイミングを検出する。例えば、検出部132は、ロール角速度の歩行波形において、起点から爪先離地のタイミングに亘る曲線における加速変曲点のタイミングを、反対足踵接地のタイミングとして検出する。例えば、検出部132は、ロール角速度の歩行波形において、踵接地のタイミングから終点に亘る曲線における減速変曲点のタイミングを、反対足爪先離地のタイミングとして検出する。For example, the detection unit 132 detects a walking event based on a walking waveform starting from the timing of heel lift. The detection unit 132 detects the timing at which a valley appears between two peaks included in the maximum peak in the walking waveform of the forward acceleration (Y-direction acceleration) as the timing of toe-off. For example, the detection unit 132 detects the minimum peak at about 60% of the walking cycle in the walking waveform of the forward acceleration (Y-direction acceleration) as the timing of sudden deceleration of the foot at the end of the swing phase. For example, the detection unit 132 detects the maximum peak at about 70% of the walking cycle in the walking waveform of the forward acceleration (Y-direction acceleration) as the timing of heel rocker. For example, the detection unit 132 detects the midpoint between the minimum peak and the maximum peak in the walking waveform of the forward acceleration (Y-direction acceleration) as the timing of heel contact. For example, the detection unit 132 detects the timing of the maximum peak between the toe-off and the heel-on-ground in the walking waveform of the vertical acceleration (Z-direction acceleration) as the timing of the tibia vertical. For example, the detection unit 132 detects the timing of a gentle maximum peak between the toe-off and the tibia vertical in the walking waveform of the forward acceleration (Y-direction acceleration) as the timing of the foot crossing. For example, the detection unit 132 detects the timing of the opposite foot heel-on-ground and the opposite foot toe-off from the walking waveform of the roll angular velocity. For example, the detection unit 132 detects the timing of the acceleration inflection point in the curve from the starting point to the toe-off timing in the walking waveform of the roll angular velocity as the timing of the opposite foot heel-on-ground. For example, the detection unit 132 detects the timing of the deceleration inflection point in the curve from the heel-on-ground timing to the end point in the walking waveform of the roll angular velocity as the timing of the opposite foot toe-off.

歩容計測部133は、検出部132によって検出された歩行ピリオドや歩行イベントに関する情報を取得する。歩容計測部133は、検出部132から取得される情報を用いて、歩容に関する計測を行う。例えば、歩容計測部133は、歩行波形から検出された歩行イベントに基づいて、ユーザの歩容を計測する。例えば、歩容計測部133は、検出された歩行イベントに基づいて、ストライド長やステップ長、歩幅、歩隔、足角、歩行速度などの歩容を計測する。例えば、歩容計測部133は、計測された歩容に基づいて、回内/回外の度合や、外反母趾の進行度、身体の対称性、身体の柔軟度などの身体的な状態を推定してもよい。例えば、歩容計測部133は、計測された歩容に基づいて、筋力低下状況や、骨密度、基礎代謝などの身体的な状態を推定してもよい。例えば、歩容計測部133は、計測された歩容に基づいて、上肢や下肢、上腕、前腕、上腿、下腿の長さなどの身体パラメータを計測してもよい。歩容計測部133による計測対象の項目には、特に限定を加えない。The gait measurement unit 133 acquires information on the walking period and walking events detected by the detection unit 132. The gait measurement unit 133 performs measurements related to the gait using the information acquired from the detection unit 132. For example, the gait measurement unit 133 measures the user's gait based on the walking events detected from the walking waveform. For example, the gait measurement unit 133 measures gait such as stride length, step length, stride width, foot angle, and walking speed based on the detected walking events. For example, the gait measurement unit 133 may estimate physical conditions such as the degree of pronation/supination, the progress of hallux valgus, body symmetry, and body flexibility based on the measured gait. For example, the gait measurement unit 133 may estimate physical conditions such as muscle weakness, bone density, and basal metabolism based on the measured gait. For example, the gait measurement unit 133 may measure body parameters such as the lengths of the upper limbs, lower limbs, upper arms, forearms, upper legs, and lower legs based on the measured gait. There are no particular limitations on the items to be measured by the gait measurement unit 133.

歩容計測装置13は、歩容計測部133によって計測された結果(歩容情報とも呼ぶ)を出力する。例えば、歩容計測装置13は、表示装置(図示しない)や携帯端末(図示しない)に歩容情報を出力する。表示装置に出力された歩容情報は、表示装置や携帯端末の画面に表示される。例えば、歩容計測装置13は、外部システム(図示しない)に歩容情報を出力する。歩容計測装置13から出力される歩容情報は、任意の用途に使用できる。歩容計測装置13が歩容情報を出力する通信機能については、特に限定を加えない。 The gait measurement device 13 outputs the results (also referred to as gait information) measured by the gait measurement unit 133. For example, the gait measurement device 13 outputs the gait information to a display device (not shown) or a mobile terminal (not shown). The gait information output to the display device is displayed on the screen of the display device or the mobile terminal. For example, the gait measurement device 13 outputs the gait information to an external system (not shown). The gait information output from the gait measurement device 13 can be used for any purpose. There are no particular limitations on the communication function through which the gait measurement device 13 outputs the gait information.

(動作)
次に、歩容計測システム1の補間装置12の動作について、図面を参照しながら説明する。計測装置11および歩容計測装置13の動作については、説明を省略する。以下においては、補間装置12に含まれる、歩容情報処理部122、欠損情報処理部123、および補間部126の行う処理について個別に説明する。以下の補間装置12の動作には、上述した補間装置12とは異なる処理や順序が含まれる場合がある。
(Operation)
Next, the operation of the interpolation device 12 of the gait measurement system 1 will be described with reference to the drawings. Descriptions of the operations of the measurement device 11 and the gait measurement device 13 will be omitted. Below, the processes performed by the gait information processing unit 122, the missing information processing unit 123, and the interpolation unit 126 included in the interpolation device 12 will be described individually. The following operation of the interpolation device 12 may include processes and sequences different from those of the interpolation device 12 described above.

〔歩容情報処理〕
まず、歩容情報処理部122による歩容情報処理について、フローチャートを参照しながら説明する。図13は、歩容情報処理部122による歩容情報処理について説明するためのフローチャートである。図13のフローチャートに沿った処理の説明においては、歩容情報処理部122を動作主体として説明する。
[Gait Information Processing]
First, the gait information processing by the gait information processing unit 122 will be described with reference to a flowchart. Fig. 13 is a flowchart for explaining the gait information processing by the gait information processing unit 122. In explaining the processing according to the flowchart of Fig. 13, the gait information processing unit 122 will be described as the subject of action.

図13において、まず、歩容情報処理部122は、センサデータの時系列データを取得する(ステップS111)。In FIG. 13, first, the gait information processing unit 122 acquires time series data of sensor data (step S111).

時系列データに欠損がある場合(ステップS112でYes)、歩容情報処理部122は、欠損を含む時系列データを欠損情報処理部123に出力する(ステップS113)。歩容情報処理部122は、欠損を含む時系列データを記憶部125に記憶させてもよい。ステップS113の次は、ステップS118に進む。If there is a missing part in the time series data (Yes in step S112), the gait information processing unit 122 outputs the time series data including the missing part to the missing part information processing unit 123 (step S113). The gait information processing unit 122 may store the time series data including the missing part in the storage unit 125. After step S113, the process proceeds to step S118.

一方、時系列データに欠損がない場合(ステップS112でNo)、歩容情報処理部122は、取得した時系列データから歩行イベントを検出する(ステップS114)。例えば、歩容情報処理部122は、踵接地や爪先離地、踵持ち上がりなどの歩行イベントを歩行波形から検出する。On the other hand, if there is no missing data in the time series data (No in step S112), the gait information processing unit 122 detects walking events from the acquired time series data (step S114). For example, the gait information processing unit 122 detects walking events such as heel contact, toe lift, and heel lift from the walking waveform.

ステップS114の次に、歩容情報処理部122は、検出された歩行イベントに基づいて、一歩行周期分の歩行波形を切り出す(ステップS115)。例えば、歩容情報処理部122は、踵接地を起点とする一歩行周期分の歩行波形を切り出す。例えば、歩容情報処理部122は、爪先離地を起点とする一歩行周期分の歩行波形を切り出す。 After step S114, the gait information processing unit 122 extracts a gait waveform for one step cycle based on the detected walking event (step S115). For example, the gait information processing unit 122 extracts a gait waveform for one step cycle starting from heel contact. For example, the gait information processing unit 122 extracts a gait waveform for one step cycle starting from toe off.

次に、歩容情報処理部122は、切り出された一周期分の歩行波形を正規化する(ステップS116)。例えば、歩容情報処理部122は、歩行波形の歩行フェーズや振幅を正規化する。Next, the gait information processing unit 122 normalizes the extracted walking waveform for one cycle (step S116). For example, the gait information processing unit 122 normalizes the walking phase and amplitude of the walking waveform.

次に、歩容情報処理部122は、正規化された歩行波形を記憶部125に記憶させる(ステップS117)。記憶部125に記憶された歩行波形は、欠損波形に含まれる欠損区間の補間や、歩容計測装置13による歩容計測に用いられる。Next, the gait information processing unit 122 stores the normalized gait waveform in the memory unit 125 (step S117). The gait waveform stored in the memory unit 125 is used to interpolate missing sections included in the missing waveform and for gait measurement by the gait measurement device 13.

処理を継続する場合(ステップS118でYes)、ステップS111に戻る。処理を継続しない場合(ステップS118でNo)、図13のフローチャートに沿った処理は終了である。ステップS118の処理の継続を決定する条件に関しては、予め設定しておけばよい。If the process is to be continued (Yes in step S118), the process returns to step S111. If the process is not to be continued (No in step S118), the process according to the flowchart in FIG. 13 ends. The conditions for determining whether to continue the process in step S118 may be set in advance.

〔欠損情報処理〕
次に、欠損情報処理部123による欠損情報処理について、フローチャートを参照しながら説明する。図14は、欠損情報処理部123による欠損情報処理について説明するためのフローチャートである。図14のフローチャートに沿った処理の説明においては、欠損情報処理部123を動作主体として説明する。
[Missing information processing]
Next, the missing information processing by the missing information processing unit 123 will be described with reference to a flowchart. Fig. 14 is a flowchart for explaining the missing information processing by the missing information processing unit 123. In the explanation of the processing according to the flowchart of Fig. 14, the missing information processing unit 123 will be described as the subject of operation.

図14において、まず、欠損情報処理部123は、欠損を含む時系列データ(欠損波形)を歩容情報処理部122から取得する(ステップS121)。欠損情報処理部123は、記憶部125から欠損波形を取得してもよい。14, first, the missing information processing unit 123 acquires time series data including missing parts (missing waveform) from the gait information processing unit 122 (step S121). The missing information processing unit 123 may acquire the missing waveform from the memory unit 125.

次に、欠損情報処理部123は、タイムスタンプに基づいて、取得した欠損波形に含まれる欠損区間を特定する(ステップS122)。Next, the missing information processing unit 123 identifies missing sections contained in the acquired missing waveform based on the timestamp (step S122).

次に、欠損情報処理部123は、特定された欠損区間の歩行フェーズを計算する(ステップS123)。Next, the missing information processing unit 123 calculates the walking phase of the identified missing section (step S123).

次に、欠損情報処理部123は、欠損波形に含まれる欠損区間の歩行フェーズに関する情報(欠損情報)を記憶部125に記憶させる(ステップS124)。記憶部125に記憶された欠損情報は、補間部126による欠損区間の補間に用いられる。Next, the missing information processing unit 123 stores information (missing information) about the walking phase of the missing section included in the missing waveform in the storage unit 125 (step S124). The missing information stored in the storage unit 125 is used by the interpolation unit 126 to interpolate the missing section.

〔補間処理〕
次に、補間部126による補間処理について、フローチャートを参照しながら説明する。ここでは、全ての区間において補間データを生成する例(図15)と、一部の区間において解析的な手法を用いて補間する例(図16)とについて説明する。
[Interpolation Processing]
Next, the interpolation process by the interpolation unit 126 will be described with reference to a flowchart. Here, an example in which interpolation data is generated in all sections (FIG. 15) and an example in which interpolation is performed in some sections using an analytical method (FIG. 16) will be described.

図15は、全ての区間において補間データを生成する場合における、補間部126の補間処理について説明するためのフローチャートである。図15のフローチャートに沿った説明においては、補間部126を動作主体として説明する。 Figure 15 is a flowchart for explaining the interpolation process of the interpolation unit 126 when generating interpolation data for all sections. In the explanation following the flowchart of Figure 15, the interpolation unit 126 will be described as the subject of operation.

図15において、まず、補間部126は、欠損区間を含む時系列データ(欠損波形)を記憶部125から取得する(ステップS131)。補間部126は、歩容情報処理部122から欠損波形を取得してもよい。欠損波形には、欠損区間の歩行フェーズを示す欠損情報が対応付けられる。15, first, the interpolation unit 126 acquires time series data including a missing section (missing waveform) from the storage unit 125 (step S131). The interpolation unit 126 may acquire the missing waveform from the gait information processing unit 122. The missing waveform is associated with missing information indicating the walking phase of the missing section.

次に、補間部126は、欠損区間の歩行フェーズに欠損がない、他の歩行周期の歩行波形を記憶部125から取得する(ステップS132)。Next, the interpolation unit 126 acquires from the memory unit 125 a walking waveform of another walking cycle in which there is no missing walking phase in the missing section (step S132).

次に、補間部126は、取得された歩行波形のデータを用いて、欠損区間の補間データを生成する(ステップS133)。Next, the interpolation unit 126 uses the acquired walking waveform data to generate interpolation data for the missing section (step S133).

次に、補間部126は、補間データの起点と終点を補正して、補間データを欠損区間に挿入する(ステップS134)。Next, the interpolation unit 126 corrects the start and end points of the interpolated data and inserts the interpolated data into the missing section (step S134).

次に、補間部126は、補正された波形をリサンプリングして、時系列データに戻す(ステップS135)。ステップS135は、歩行速度や歩幅などを計算する際に、元の時系列データが必要となるために行われる。元の時系列データが必要ない場合、ステップS135は省略してもよい。Next, the interpolation unit 126 resamples the corrected waveform to return it to time series data (step S135). Step S135 is performed because the original time series data is required when calculating walking speed, stride length, etc. If the original time series data is not required, step S135 may be omitted.

図16は、一部の区間において解析的な手法を用いて補間する場合における、補間部126の補間処理について説明するためのフローチャートである。図16のフローチャートに沿った説明においては、補間部126を動作主体として説明する。 Figure 16 is a flowchart for explaining the interpolation process of the interpolation unit 126 when an analytical method is used to interpolate in some sections. In the explanation following the flowchart of Figure 16, the interpolation unit 126 will be described as the subject of operation.

図16において、まず、補間部126は、欠損区間を含む時系列データ(欠損波形)を記憶部125から取得する(ステップS141)。補間部126は、歩容情報処理部122から欠損波形を取得してもよい。欠損波形には、欠損区間の歩行フェーズを示す欠損情報が対応付けられる。16, first, the interpolation unit 126 acquires time series data including a missing section (missing waveform) from the storage unit 125 (step S141). The interpolation unit 126 may acquire the missing waveform from the gait information processing unit 122. The missing waveform is associated with missing information indicating the walking phase of the missing section.

解析的な手法を使用可能な欠損区間である場合(ステップS142でYes)、補間部126は、解析的な手法で欠損区間のデータを補間する(ステップS143)。ステップS143の次は、ステップS147に進む。If the missing section is one for which an analytical method can be used (Yes in step S142), the interpolation unit 126 interpolates the data of the missing section using an analytical method (step S143). After step S143, the process proceeds to step S147.

一方、解析的な手法を使用可能な欠損区間ではない場合(ステップS142でNo)、補間部126は、欠損区間の歩行フェーズに欠損がない、他の歩行周期の歩行波形を記憶部125から取得する(ステップS144)。On the other hand, if the missing section is not one for which analytical methods can be used (No in step S142), the interpolation unit 126 acquires from the memory unit 125 a walking waveform of another walking cycle in which there is no missing part in the walking phase of the missing section (step S144).

次に、補間部126は、取得された歩行波形のデータを用いて、欠損区間の補間データを生成する(ステップS145)。Next, the interpolation unit 126 uses the acquired walking waveform data to generate interpolation data for the missing section (step S145).

次に、補間部126は、補間データの起点と終点を補正して、補間データを欠損区間に挿入する(ステップS146)。Next, the interpolation unit 126 corrects the start and end points of the interpolated data and inserts the interpolated data into the missing section (step S146).

次に、補間部126は、補正された波形をリサンプリングして、時系列データに戻す(ステップS146)。ステップS146は、歩行速度や歩幅などを計算する際に、元の時系列データが必要となるために行われる。元の時系列データが必要ない場合、ステップS146は省略してもよい。Next, the interpolation unit 126 resamples the corrected waveform to return it to time series data (step S146). Step S146 is performed because the original time series data is required when calculating walking speed, stride length, etc. If the original time series data is not required, step S146 may be omitted.

以上のように、本実施形態の歩容計測システムは、計測装置、補間装置、および歩容計測装置を備える。計測装置は、ユーザの履物に配置される。計測装置は、ユーザの歩行に応じて空間加速度および空間角速度を計測する。計測装置は、計測された空間加速度および空間角速度に基づくセンサデータを生成する。計測装置は、生成されたセンサデータを補間装置に出力する。補間装置は、受信部、歩容情報処理部、欠損情報処理部、記憶部、補間部、および送信部を有する。受信部は、計測装置から送信されたセンサデータを受信する。歩容情報処理部は、センサデータの時系列データを生成する。歩容情報処理部は、足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成するとともに、時系列データにおけるデータの欠損区間を特定する。歩容情報処理部は、生成された歩行波形や、特定された欠損区間に関する情報を記憶部に記憶させる。記憶部には、歩行波形や、欠損区間に関する情報、欠損区間が補間された修復波形が記憶される。欠損情報処理部は、特定された欠損区間の歩行フェーズを計算する。補間部は、欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、欠損区間の歩行フェーズのデータを用いて、欠損区間を補間する補間データを生成する。補間部は、生成された補間データを欠損区間に補間する。
歩容計測装置は、センサデータの時系列データを用いて補間装置が生成した欠損区間を含まない歩行波形と、補間装置が欠損区間のデータを補間した歩行波形とを補間装置から取得する。歩容計測装置は、取得した歩行波形から検出される歩行イベントに基づいてユーザの歩容を計測する。歩容計測装置は、計測されたユーザの歩容に関する情報を出力する。
As described above, the gait measurement system of this embodiment includes a measurement device, an interpolation device, and a gait measurement device. The measurement device is placed on a user's footwear. The measurement device measures spatial acceleration and spatial angular velocity according to the user's gait. The measurement device generates sensor data based on the measured spatial acceleration and spatial angular velocity. The measurement device outputs the generated sensor data to the interpolation device. The interpolation device has a receiving unit, a gait information processing unit, a missing information processing unit, a storage unit, an interpolation unit, and a transmission unit. The receiving unit receives the sensor data transmitted from the measurement device. The gait information processing unit generates time-series data of the sensor data. The gait information processing unit generates a gait waveform for each gait cycle using the time-series data of the sensor data related to the movement of the feet, and identifies a missing section of data in the time-series data. The gait information processing unit stores the generated gait waveform and information on the identified missing section in the storage unit. The storage unit stores the gait waveform, information on the missing section, and a repair waveform in which the missing section is interpolated. The missing information processing unit calculates the gait phase of the identified missing section. The interpolation unit generates interpolation data for interpolating the missing section by using data of a walking phase of the missing section in a walking waveform having a walking cycle different from the walking cycle including the missing section. The interpolation unit interpolates the missing section with the generated interpolation data.
The gait measurement device acquires from the interpolation device a gait waveform that does not include a missing section that is generated by the interpolation device using time-series data of the sensor data, and a gait waveform in which the interpolation device interpolates data of the missing section. The gait measurement device measures the user's gait based on walking events detected from the acquired gait waveform. The gait measurement device outputs information related to the measured gait of the user.

本実施形態の手法では、歩行の周期性に着目して、ユーザの歩行に応じて計測されるセンサデータの時系列データに発生したデータの欠損区間を補間するための補間データを生成する。本実施形態の手法では、異なる歩行周期の歩行波形の同じ歩行フェーズのデータを用いて、補間データを生成する。そのため、本実施形態の手法によれば、歩行の周期性に着目して生成された補間データを用いるため、センサデータの時系列データに含まれるデータの欠損区間の特徴を含めて、欠損区間のデータを補間できる。 In the method of this embodiment, by focusing on the periodicity of walking, interpolation data is generated to fill in missing sections of data that occur in the time series data of sensor data measured in response to the user's walking. In the method of this embodiment, interpolation data is generated using data from the same walking phase of walking waveforms with different walking cycles. Therefore, according to the method of this embodiment, since interpolation data generated by focusing on the periodicity of walking is used, data of missing sections can be interpolated, including the characteristics of the missing sections of data contained in the time series data of sensor data.

無線通信機能を有するウエラブルセンサによって計測されたセンサデータを通信する場合、通信障害によりデータの欠損が起こりうる。通信障害で起こりうるデータの欠損は、離散的ではなく、連続的であることが多い。例えば、センサデータの時系列データから、特徴抽出などに用いられる重要なデータが丸ごと欠損することもありうる。本実施形態の手法によれば、歩行動作の相似性と周期性を利用してデータの欠損区間を補間するため、数%から十数%の歩行フェーズに亘る長い欠損区間のデータを補間できる。また、本実施形態の手法は、回路的な要因で計測された外れ値の除去や、原因不明の要因で計測された異常な値の平準化にも用いることができる。When communicating sensor data measured by a wearable sensor with wireless communication capabilities, data loss may occur due to communication failure. Data loss due to communication failure is often continuous, not discrete. For example, important data used for feature extraction may be lost entirely from time-series sensor data. According to the method of this embodiment, the similarity and periodicity of walking movements are used to interpolate data loss intervals, so that data loss intervals spanning several percent to several tens of percent of the walking phase can be interpolated. In addition, the method of this embodiment can also be used to remove outliers measured due to circuitry-related factors and to smooth out abnormal values measured due to unknown factors.

本実施形態の一態様において、歩容情報処理部は、時系列データから一歩行周期分の波形を切り出す。歩容情報処理部は、切り出された波形に含まれる歩行フェーズを正規化して歩行波形を生成する。本態様によれば、センサデータの時系列データから生成される複数の歩行波形の歩行フェーズが統一されるので、複数の歩行波形の間で相互にデータを補間しやすくなる。In one aspect of this embodiment, the gait information processing unit extracts a waveform for one step walking cycle from the time series data. The gait information processing unit normalizes the walking phase included in the extracted waveform to generate a walking waveform. According to this aspect, the walking phases of multiple walking waveforms generated from the time series data of the sensor data are unified, making it easier to mutually interpolate data between the multiple walking waveforms.

本実施形態の一態様において、補間部は、時系列データに含まれる欠損区間に関して、欠損区間に含まれる歩行フェーズのデータが欠損していない少なくとも一つの歩行波形から、欠損区間に含まれる前記歩行フェーズのデータを少なくとも一つ選択する。補間部は、選択された歩行フェーズのデータを用いて、補間データを生成する。本態様によれば、他の歩行周期のデータを用いて、欠損区間を補間できる。In one aspect of this embodiment, the interpolation unit selects, for a missing section included in the time series data, at least one piece of data for the walking phase included in the missing section from at least one walking waveform in which data for the walking phase included in the missing section is not missing. The interpolation unit generates interpolated data using the data for the selected walking phase. According to this aspect, the missing section can be interpolated using data for other walking cycles.

本実施形態の一態様において、補間部は、複数の歩行周期の歩行波形から選択された歩行フェーズのデータの平均値を計算する。補間部は、算出された歩行フェーズのデータの平均値を用いて、補間データを生成する。本態様によれば、複数の歩行周期のデータの平均値に基づいて欠損区間のデータを補間するため、歩行の特徴がより反映された補間データを生成できる。In one aspect of this embodiment, the interpolation unit calculates the average value of data for a walking phase selected from walking waveforms of multiple walking cycles. The interpolation unit generates interpolated data using the average value of the calculated walking phase data. According to this aspect, data for missing sections is interpolated based on the average value of data for multiple walking cycles, so that interpolated data that better reflects the characteristics of walking can be generated.

本実施形態の一態様において、補間部は、データの変動が小さい歩行ピリオドに関しては、解析的な手法で補間データを生成する。補間部は、データの変動が大きい歩行ピリオドに関しては、欠損区間を含む歩行周期とは異なる歩行周期のデータを用いて補間データを生成する。本態様によれば、データの変動が小さい歩行ピリオドに関しては、解析的な手法で補間データを生成するため、計算時間や計算量を低減できる。 In one aspect of this embodiment, the interpolation unit generates interpolation data using an analytical method for walking periods with small data fluctuations. For walking periods with large data fluctuations, the interpolation unit generates interpolation data using data from a walking cycle that is different from the walking cycle that includes the missing section. According to this aspect, for walking periods with small data fluctuations, the interpolation data is generated using an analytical method, which reduces the calculation time and amount of calculations.

(第2の実施形態)
次に、第2の実施形態に係る歩容計測システムについて図面を参照しながら説明する。本実施形態の歩容計測システムは、複数の歩行周期の歩行波形に関する歩行フェーズごとのデータ値(振幅の値)の全体分布における偏差値に基づいて、欠損区間の補間データを生成する点において、第1の実施形態とは異なる。
Second Embodiment
Next, a gait measurement system according to a second embodiment will be described with reference to the drawings. The gait measurement system of this embodiment differs from the first embodiment in that it generates interpolation data for missing sections based on the deviation value in the overall distribution of data values (amplitude values) for each walking phase related to the walking waveform of a plurality of walking cycles.

(構成)
図17は、本実施形態の歩容計測システム2の構成を示すブロック図である。歩容計測システム2は、計測装置21、補間装置22、および歩容計測装置23を備える。補間装置22は、計測装置21および歩容計測装置23に有線で接続されてもよいし、無線で接続されてもよい。また、計測装置21、補間装置22、および歩容計測装置23は、単一の装置で構成されてもよい。また、歩容計測システム2は、計測装置21を除き、補間装置22および歩容計測装置23で構成されてもよい。
(composition)
17 is a block diagram showing the configuration of a gait measurement system 2 of this embodiment. The gait measurement system 2 includes a measurement device 21, an interpolation device 22, and a gait measurement device 23. The interpolation device 22 may be connected to the measurement device 21 and the gait measurement device 23 by wire or wirelessly. The measurement device 21, the interpolation device 22, and the gait measurement device 23 may be configured as a single device. The gait measurement system 2 may be configured from the interpolation device 22 and the gait measurement device 23, excluding the measurement device 21.

計測装置21は、第1の実施形態の計測装置11と同様の構成である。計測装置21は、足部に設置される。計測装置21は、履物を履くユーザの足の動きに関する物理量として、加速度センサによって計測される加速度(空間加速度とも呼ぶ)や、角速度センサによって計測される角速度(空間角速度とも呼ぶ)を計測する。計測装置21が計測する足の動きに関する物理量には、加速度や角速度を積分することによって計算される速度や角度、位置(軌跡)も含まれる。計測装置21は、計測された物理量をデジタルデータ(センサデータとも呼ぶ)に変換する。計測装置21は、変換後のセンサデータを補間装置22に送信する。The measuring device 21 has the same configuration as the measuring device 11 of the first embodiment. The measuring device 21 is installed on the foot. The measuring device 21 measures acceleration (also called spatial acceleration) measured by an acceleration sensor and angular velocity (also called spatial angular velocity) measured by an angular velocity sensor as physical quantities related to the movement of the foot of a user wearing the footwear. The physical quantities related to the movement of the foot measured by the measuring device 21 also include velocity, angle, and position (trajectory) calculated by integrating the acceleration and angular velocity. The measuring device 21 converts the measured physical quantities into digital data (also called sensor data). The measuring device 21 transmits the converted sensor data to the interpolation device 22.

補間装置22は、計測装置21からセンサデータを受信する。補間装置22は、受信したセンサデータの時系列データを生成する。補間装置22は、センサデータに含まれるタイムスタンプに基づいて、時系列データを生成する。補間装置22は、生成された時系列データから一歩行周期分の波形を切り出す。例えば、補間装置22は、踵接地や踵持ち上がりのタイミングを起点とする一歩行周期分の波形を抽出する。補間装置22は、切り出された一歩行周期分の波形を正規化して、一歩行周期分の波形(歩行波形とも呼ぶ)を生成する。 The interpolation device 22 receives sensor data from the measurement device 21. The interpolation device 22 generates time series data of the received sensor data. The interpolation device 22 generates the time series data based on the timestamp included in the sensor data. The interpolation device 22 extracts a waveform for one step walking cycle from the generated time series data. For example, the interpolation device 22 extracts a waveform for one step walking cycle starting from the timing of heel strike or heel lift. The interpolation device 22 normalizes the extracted waveform for one step walking cycle to generate a waveform for one step walking cycle (also called a walking waveform).

補間装置22は、センサデータの時系列データにおいて、タイムスタンプの番号が抜けている部分を欠損区間として特定する。補間装置22は、欠損区間に含まれる歩行フェーズのデータ値の全体分布に基づいて、欠損区間のデータを補間するための補間データを生成する。全体分布とは、欠損区間の補間に用いられる複数の歩行周期の歩行波形を全体とする、欠損区間に含まれる歩行フェーズごとのデータ値の分布である。言い換えると、補間装置22は、欠損区間の歩行フェーズに欠損がない複数の歩行周期の歩行波形のデータ値(振幅の値)を用いて、欠損区間に対応する歩行フェーズのデータ値の全体分布における統計量を計算する。補間装置22は、欠損区間に含まれる歩行フェーズごとのデータ値の全体分布における統計量に基づいて、欠損区間の補間データを生成する。特に、補間装置22は、欠損区間に含まれる歩行フェーズごとのデータ値の全体分布における代表値に基づいて、欠損区間の補間データを生成する。例えば、補間装置22は、欠損区間に含まれる歩行フェーズごとのデータ値(振幅の値)の全体分布における偏差値に基づいて、欠損区間の補間データを生成する。例えば、補間装置22は、偏差値の代わりに、標準偏差や分散、偏差などに基づいて、欠損区間の補間データを生成してもよい。補間装置22は、生成された補間データを用いて、欠損区間を補間する。The interpolation device 22 identifies a portion of the time series data of the sensor data where the number of the time stamp is missing as a missing section. The interpolation device 22 generates interpolation data for interpolating the data of the missing section based on the overall distribution of the data values of the walking phase included in the missing section. The overall distribution is the distribution of the data values for each walking phase included in the missing section, with the walking waveforms of the multiple walking cycles used to interpolate the missing section as the whole. In other words, the interpolation device 22 calculates the statistics in the overall distribution of the data values of the walking phase corresponding to the missing section using the data values (amplitude values) of the walking waveforms of the multiple walking cycles in which the walking phase of the missing section has no missing data. The interpolation device 22 generates the interpolation data for the missing section based on the statistics in the overall distribution of the data values for each walking phase included in the missing section. In particular, the interpolation device 22 generates the interpolation data for the missing section based on a representative value in the overall distribution of the data values for each walking phase included in the missing section. For example, the interpolation device 22 generates the interpolation data for the missing section based on the deviation value in the overall distribution of the data values (amplitude values) for each walking phase included in the missing section. For example, the interpolation device 22 may generate the interpolation data for the missing section based on the standard deviation, variance, deviation, etc. instead of the deviation value. The interpolation device 22 uses the generated interpolation data to interpolate the missing section.

歩容計測装置23は、第1の実施形態の歩容計測装置13と同様の構成である。歩容計測装置23は、補間装置22から、欠損区間を含まない歩行波形を取得する。歩容計測装置23は、取得した歩行波形を用いて、歩容に関する計測を実行する。歩容計測装置23は、歩容に関する計測結果(歩容情報とも呼ぶ)を出力する。 The gait measurement device 23 has the same configuration as the gait measurement device 13 of the first embodiment. The gait measurement device 23 acquires a gait waveform that does not include a missing section from the interpolation device 22. The gait measurement device 23 performs measurements related to gait using the acquired gait waveform. The gait measurement device 23 outputs measurement results related to gait (also called gait information).

例えば、補間装置22および歩容計測装置23は、図示しないサーバ等に実装される。例えば、補間装置22および歩容計測装置23は、アプリケーションサーバによって実現されてもよい。例えば、補間装置22および歩容計測装置23は、携帯端末(図示しない)にインストールされたアプリケーションソフトウェア等によって実現されてもよい。For example, the interpolation device 22 and the gait measurement device 23 are implemented in a server or the like (not shown). For example, the interpolation device 22 and the gait measurement device 23 may be realized by an application server. For example, the interpolation device 22 and the gait measurement device 23 may be realized by application software or the like installed on a mobile terminal (not shown).

〔補間装置〕
次に、補間装置22の詳細構成について図面を参照しながら説明する。図18は、補間装置22の詳細構成の一例を示すブロックである。補間装置22は、受信部221、歩容情報処理部222、欠損情報処理部223、記憶部225、補間部226、および送信部228を備える。
[Interpolation device]
Next, a detailed configuration of the interpolation device 22 will be described with reference to the drawings. Fig. 18 is a block diagram showing an example of the detailed configuration of the interpolation device 22. The interpolation device 22 includes a receiving unit 221, a gait information processing unit 222, a missing information processing unit 223, a storage unit 225, an interpolation unit 226, and a transmitting unit 228.

受信部221は、第1の実施形態の受信部121と同様の構成である。受信部221は、計測装置(図示しない)からセンサデータを受信する。受信部221は、受信したセンサデータを歩容情報処理部222に出力する。The receiving unit 221 has a configuration similar to that of the receiving unit 121 of the first embodiment. The receiving unit 221 receives sensor data from a measuring device (not shown). The receiving unit 221 outputs the received sensor data to the gait information processing unit 222.

歩容情報処理部222は、第1の実施形態の歩容情報処理部122と同様の構成である。歩容情報処理部222は、受信部221からセンサデータを取得する。歩容情報処理部222は、取得したセンサデータの時系列データを生成する。歩容情報処理部222は、センサデータに含まれるタイムスタンプに基づいて、時系列データを生成する。歩容情報処理部222は、生成された時系列データから一歩行周期分の波形を切り出す。歩容情報処理部222は、切り出された一歩行周期分の波形を正規化して、一歩行周期分の歩行波形を生成する。歩容情報処理部222は、正規化された歩行波形を記憶部225に記憶させる。 The gait information processing unit 222 has a configuration similar to that of the gait information processing unit 122 of the first embodiment. The gait information processing unit 222 acquires sensor data from the receiving unit 221. The gait information processing unit 222 generates time series data of the acquired sensor data. The gait information processing unit 222 generates the time series data based on a timestamp included in the sensor data. The gait information processing unit 222 extracts a waveform for one step gait cycle from the generated time series data. The gait information processing unit 222 normalizes the extracted waveform for one step gait cycle to generate a gait waveform for one step gait cycle. The gait information processing unit 222 stores the normalized gait waveform in the memory unit 225.

欠損情報処理部223は、第1の実施形態の欠損情報処理部123と同様の構成である。欠損情報処理部223は、欠損波形を取得する。欠損情報処理部223は、欠損波形に含まれる欠損区間に関する情報(欠損情報とも呼ぶ)を、その欠損波形に対応付けて記憶部225に記憶させる。The missing information processing unit 223 has a configuration similar to that of the missing information processing unit 123 of the first embodiment. The missing information processing unit 223 acquires a missing waveform. The missing information processing unit 223 stores information about a missing section included in the missing waveform (also called missing information) in the memory unit 225 in association with the missing waveform.

記憶部225は、第1の実施形態の記憶部125と同様の構成である。記憶部225には、歩行波形が記憶される。記憶部225に記憶される歩行波形は、欠損区間を含まない歩行波形(正常波形とも呼ぶ)や、欠損区間を含む歩行波形(欠損波形とも呼ぶ)、欠損区間が修復された歩行波形(修復波形とも呼ぶ)が含まれる。欠損波形には、その欠損波形の欠損情報が対応付けられる。記憶部225に記憶された歩行波形は、所定のタイミングで送信部228によって取得され、送信部228から歩容計測装置(図示しない)に送信される。The memory unit 225 has the same configuration as the memory unit 125 of the first embodiment. A gait waveform is stored in the memory unit 225. The gait waveforms stored in the memory unit 225 include a gait waveform that does not include a missing section (also called a normal waveform), a gait waveform that includes a missing section (also called a missing waveform), and a gait waveform in which the missing section has been repaired (also called a repaired waveform). The missing waveform is associated with the missing information of the missing waveform. The gait waveform stored in the memory unit 225 is acquired by the transmission unit 228 at a predetermined timing and transmitted from the transmission unit 228 to a gait measurement device (not shown).

補間部226は、補間対象の欠損波形に含まれる欠損区間の歩行フェーズに欠損がない複数の歩行波形を、記憶部225から取得する。補間部226は、取得された複数の歩行波形のデータ値(振幅の値)を用いて、複数の歩行波形のデータ群における、欠損区間に対応する歩行フェーズのデータ値の全体分布における統計量を計算する。特に、補間部226は、統計量として、偏差値や、標準偏差、分散、偏差などの代表値を計算する。例えば、補間部226は、複数の歩行波形のデータ群の全体分布における、欠損区間に含まれる歩行フェーズのデータ値の偏差値を計算する。補間部226は、算出された統計量に基づいて、補間データを生成する。例えば、補間部226は、偏差値に基づいて、補間データを生成する。例えば、補間部226は、偏差値の代わりに、標準偏差や分散、偏差などに基づいて、欠損区間の補間データを生成してもよい。補間部226は、生成された補間データを用いて、欠損区間を補間する。The interpolation unit 226 acquires from the storage unit 225 a plurality of walking waveforms having no missing parts in the walking phase of the missing section included in the missing waveform to be interpolated. The interpolation unit 226 uses the data values (amplitude values) of the acquired plurality of walking waveforms to calculate statistics in the overall distribution of the data values of the walking phase corresponding to the missing section in the data group of the plurality of walking waveforms. In particular, the interpolation unit 226 calculates representative values such as deviation values, standard deviations, variances, and deviations as statistics. For example, the interpolation unit 226 calculates the deviation values of the data values of the walking phase included in the missing section in the overall distribution of the data group of the plurality of walking waveforms. The interpolation unit 226 generates the interpolation data based on the calculated statistics. For example, the interpolation unit 226 generates the interpolation data based on the deviation values. For example, the interpolation unit 226 may generate the interpolation data of the missing section based on the standard deviation, variance, deviation, etc. instead of the deviation values. The interpolation unit 226 uses the generated interpolation data to interpolate the missing section.

ここで、統計量(代表値)の一つである偏差値に基づいて補間データを生成する方法について、いくつかの図面(図19~図23)を用いて説明する。図19は、ユーザの歩行に伴って計測された進行方向加速度(Y方向加速度)の時系列データである。図19の時系列データは、20秒間の計測期間における複数の歩行周期の波形を含む。各歩行周期の波形は、データ値(振幅の値)に違いがあるものの、同様のパターンを示す。Here, a method for generating interpolated data based on a deviation value, which is one of the statistics (representative values), will be explained using several drawings (Figs. 19 to 23). Fig. 19 shows time series data of forward acceleration (Y-direction acceleration) measured as the user walks. The time series data in Fig. 19 includes waveforms of multiple walking cycles over a measurement period of 20 seconds. The waveforms of each walking cycle show a similar pattern, although there are differences in the data values (amplitude values).

図20は、図19の時系列データに含まれる複数の歩行周期の波形のうち一つを切り出し、切り出された波形を正規化した波形(歩行波形)である。図20の歩行波形には、歩行の特徴を含むパターンがみられる。図20には、後続する説明で用いられる歩行フェーズ(P1、P2)を示す。以下において、歩行フェーズP1および歩行フェーズP2は、図20の歩行周期のみならず、複数の歩行周期に関しても適用される。 Figure 20 shows a waveform (walking waveform) obtained by extracting one of the multiple walking cycle waveforms contained in the time series data of Figure 19 and normalizing the extracted waveform. The walking waveform of Figure 20 shows a pattern that includes walking characteristics. Figure 20 shows walking phases (P1, P2) that will be used in the following explanation. In the following, walking phase P1 and walking phase P2 apply not only to the walking cycle of Figure 20 but also to multiple walking cycles.

図21は、図19の時系列データから切り出された複数の歩行周期の歩行波形をサンプル群(全体)とする、歩行フェーズP1および歩行フェーズP2におけるデータ値(振幅の値)の頻度分布である。歩行フェーズP1に関しては、データ値の頻度分布が正規分布に近い分布を示す。歩行フェーズP2に関しても、データ値の頻度分布が正規分布に近い分布を示す。その他の全ての歩行フェーズに関しても、1%ごとのデータ値の頻度分布をコルモゴロフ-スミルノフ検定で正規化判別したところ、陽性結果が得られた。すなわち、複数の歩行周期の歩行波形の各歩行フェーズにおけるデータ値の頻度分布は、歩行周期に関わらず、正規分布に近い分布を示す。これは、複数の歩行周期の歩行波形に関して、歩行フェーズごとのデータ値の全体分布に基づいて、補間データを推定可能であることを示す。言い換えると、複数の歩行周期の歩行波形に関して、歩行フェーズごとのデータ値の全体分布における統計量に基づいて、補間データを推定可能である。 Figure 21 shows the frequency distribution of data values (amplitude values) in walking phases P1 and P2, with the walking waveforms of multiple walking cycles extracted from the time series data of Figure 19 as the sample group (whole). For walking phase P1, the frequency distribution of data values shows a distribution close to a normal distribution. For walking phase P2, the frequency distribution of data values shows a distribution close to a normal distribution. For all other walking phases, the frequency distribution of data values for each 1% was normalized and judged using the Kolmogorov-Smirnov test, and a positive result was obtained. In other words, the frequency distribution of data values in each walking phase of the walking waveforms of multiple walking cycles shows a distribution close to a normal distribution regardless of the walking cycle. This shows that it is possible to estimate interpolated data for walking waveforms of multiple walking cycles based on the overall distribution of data values for each walking phase. In other words, it is possible to estimate interpolated data for walking waveforms of multiple walking cycles based on statistics in the overall distribution of data values for each walking phase.

図22は、図19のセンサデータの時系列データにおいて、連続する歩行フェーズのデータNおよびデータN+1のデータ値の全体分布における偏差値の相関関係を示すグラフである(Nは自然数)。図22のように、連続するデータNとデータN+1の偏差値は、ほぼ比例関係を示す。すなわち、歩行フェーズP2におけるデータ値の全体分布における偏差値は、連続するデータNとデータN+1に関して同様の値を示す。 Figure 22 is a graph showing the correlation between the standard deviations in the overall distribution of data values of data N and data N+1 in consecutive walking phases in the time series data of the sensor data in Figure 19 (N is a natural number). As shown in Figure 22, the standard deviations of consecutive data N and data N+1 show an almost proportional relationship. In other words, the standard deviations in the overall distribution of data values in walking phase P2 show similar values for consecutive data N and data N+1.

図23は、対象歩行フェーズ(歩行フェーズP2)からの距離に応じた偏差値の相関係数を示すグラフである。図23の実線は、図19の時系列データから切り出された複数の歩行周期の歩行波形のサンプル群に関する、歩行フェーズP2のデータ値の偏差値と、歩行フェーズP2とは異なる歩行フェーズのデータ値との相関係数を示す。図23の破線は、相関係数(実線)を多項式にフィッティングさせた曲線である。図23のように、対象歩行フェーズからの距離が近い歩行フェーズほど、偏差値の相関係数が大きく、偏差値が同様の値を示すことが分かる。言い換えると、対象歩行フェーズからの距離が近い歩行フェーズほど、相関係数が大きく、偏差値の値が近くなる。すなわち、対象歩行フェーズの偏差値に近い偏差値の歩行フェーズのデータ値を用いて、その対象歩行フェーズの補間データを生成できる。23 is a graph showing the correlation coefficient of the deviation value according to the distance from the target walking phase (walking phase P2). The solid line in FIG. 23 shows the correlation coefficient between the deviation value of the data value of walking phase P2 and the data value of a walking phase different from walking phase P2 for the sample group of walking waveforms of multiple walking cycles cut out from the time series data in FIG. 19. The dashed line in FIG. 23 is a curve obtained by fitting the correlation coefficient (solid line) to a polynomial. As shown in FIG. 23, it can be seen that the closer the walking phase is from the target walking phase, the larger the correlation coefficient of the deviation value and the more similar the deviation value is. In other words, the closer the walking phase is from the target walking phase, the larger the correlation coefficient and the closer the deviation value is. In other words, the data value of a walking phase with a deviation value close to the deviation value of the target walking phase can be used to generate interpolated data for the target walking phase.

例えば、補間部226は、算出された歩行フェーズの偏差値から逆算して補間データを計算し、欠損区間を補間するための補間データを生成する。例えば、補間部226は、偏差値の代わりに、標準偏差や分散、偏差などの代表値に基づいて、補間データを生成してもよい。例えば、補間部226は、以下の式1を逆算して、歩行フェーズjにおける補間データの値djを計算する。以下の式1は、欠損区間に含まれる歩行フェーズjのデータ値の偏差値Tjの計算式である。

Figure 0007626221000001
For example, the interpolation unit 226 calculates the interpolation data by performing a reverse calculation from the deviation value of the calculated walking phase, and generates the interpolation data for interpolating the missing section. For example, the interpolation unit 226 may generate the interpolation data based on a representative value such as a standard deviation, a variance, or a deviation, instead of the deviation value. For example, the interpolation unit 226 calculates the value d j of the interpolation data in the walking phase j by performing a reverse calculation of the following formula 1. The following formula 1 is a calculation formula for the deviation value T j of the data value of the walking phase j included in the missing section.
Figure 0007626221000001

上記の式1において、μは、対象のデータ群に含まれる複数の歩行波形における歩行フェーズjのデータ値の平均値である。σは、対象のデータ群に含まれる複数の歩行波形における歩行フェーズjのデータ値の標準偏差である。偏差値の代わりに、標準偏差や分散、偏差などの代表値に基づいて補間データを生成する場合、標準偏差や分散、偏差などの計算式を逆算して、補間データの値を計算すればよい。 In the above formula 1, μ is the average value of the data values of walking phase j in the multiple walking waveforms included in the target data group. σ is the standard deviation of the data values of walking phase j in the multiple walking waveforms included in the target data group. When generating interpolated data based on representative values such as standard deviation, variance, and deviation instead of deviation values, the value of the interpolated data can be calculated by inversely calculating the formula for standard deviation, variance, and deviation.

例えば、補間部226は、補間対象の欠損波形に含まれる欠損区間の起点および終点のデータ値の偏差値に基づいて、欠損区間の補間データを推定してもよい。図24は、補間対象の欠損波形に含まれる欠損区間の起点(歩行フェーズM1)および終点(歩行フェーズM2)の偏差値に基づく欠損区間の補間の一例について説明するためのグラフである。例えば、補間部226は、欠損区間に関して、歩行フェーズM1のデータ値の全体分布における偏差値と、歩行フェーズM2のデータ値の全体分布における偏差値との平均値を計算する。補間部226は、歩行フェーズM1と歩行フェーズM2のデータ値の偏差値の平均値を用いて、データ値の全体分布から補間データを推定する。図24の手法では、欠損区間の歩行フェーズM1と歩行フェーズM2の距離が近いほど、補間データの推定精度が高い。For example, the interpolation unit 226 may estimate the interpolation data of the missing section based on the deviation value of the data value of the start point and the end point of the missing section included in the missing waveform to be interpolated. FIG. 24 is a graph for explaining an example of the interpolation of the missing section based on the deviation value of the start point (walking phase M1) and the end point (walking phase M2) of the missing section included in the missing waveform to be interpolated. For example, the interpolation unit 226 calculates the average value of the deviation value in the entire distribution of the data value of the walking phase M1 and the deviation value in the entire distribution of the data value of the walking phase M2 for the missing section. The interpolation unit 226 estimates the interpolation data from the entire distribution of the data value using the average value of the deviation values of the data values of the walking phase M1 and the walking phase M2. In the method of FIG. 24, the closer the distance between the walking phase M1 and the walking phase M2 of the missing section, the higher the estimation accuracy of the interpolation data.

例えば、本実施形態の補間部226による補間データの生成方法は、第1の実施形態の補間部126による補間データの生成方法と組み合わせてもよい。例えば、第1の実施形態の手法で補間データを生成し、生成された補間データの起点および終点を本実施形態の手法で補正してもよい。For example, the method of generating interpolated data by the interpolation unit 226 of this embodiment may be combined with the method of generating interpolated data by the interpolation unit 126 of the first embodiment. For example, interpolated data may be generated by the method of the first embodiment, and the start point and end point of the generated interpolated data may be corrected by the method of this embodiment.

補間部226は、欠損区間が補間された歩行波形(修復波形)を記憶部225に記憶させる。例えば、補間部226は、修復波形において、修復された欠損区間を特定できるように、欠損区間にラベルを付してもよい。例えば、補間部226は、欠損区間に仮の補間データを挿入した修復波形(仮修復波形とも呼ぶ)を記憶部225に記憶させてもよい。記憶部225に格納された修復波形は、他の歩行波形と同様に、歩容計測装置(図示しない)による歩容計測に用いられる。The interpolation unit 226 stores the walking waveform (repaired waveform) in which the missing sections have been interpolated in the memory unit 225. For example, the interpolation unit 226 may label the missing sections in the repaired waveform so that the repaired missing sections can be identified. For example, the interpolation unit 226 may store in the memory unit 225 a repaired waveform (also called a provisional repaired waveform) in which provisional interpolation data has been inserted into the missing sections. The repaired waveform stored in the memory unit 225 is used for gait measurement by a gait measurement device (not shown) in the same way as other walking waveforms.

送信部228は、記憶部225から歩行波形を取得する。送信部228は、正常波形および修復波形を含む歩行波形を記憶部225から取得する。送信部228は、取得した歩行波形を歩容計測装置(図示しない)に送信する。送信部228は、ケーブルなどの有線を介して歩容計測装置に歩行波形を送信してもよいし、無線通信を介して歩容計測装置に歩行波形を送信してもよい。送信部228と歩容計測装置の間の通信方式に関しては、特に限定しない。The transmission unit 228 acquires a gait waveform from the memory unit 225. The transmission unit 228 acquires a gait waveform including a normal waveform and a restored waveform from the memory unit 225. The transmission unit 228 transmits the acquired gait waveform to a gait measurement device (not shown). The transmission unit 228 may transmit the gait waveform to the gait measurement device via a wired connection such as a cable, or may transmit the gait waveform to the gait measurement device via wireless communication. There are no particular limitations on the communication method between the transmission unit 228 and the gait measurement device.

(動作)
次に、補間装置22の動作について、図面を参照しながら説明する。以下においては、補間装置22に含まれる補間部226による、欠損区間の起点および終点の偏差値に基づいて欠損区間の補間を行う例について説明する。以下の補間装置22の動作は、上述した補間部226とは異なる処理や順序を含む場合もある。歩容情報処理部222および欠損情報処理部223の処理は、第1の実施形態と同様なので省略する。
(Operation)
Next, the operation of the interpolation device 22 will be described with reference to the drawings. In the following, an example will be described in which the interpolation unit 226 included in the interpolation device 22 performs interpolation of a missing section based on the deviation values of the start point and end point of the missing section. The following operation of the interpolation device 22 may include processing and an order different from those of the above-mentioned interpolation unit 226. The processing of the gait information processing unit 222 and the missing information processing unit 223 is omitted since it is the same as in the first embodiment.

〔補間処理〕
補間部226による補間処理について、フローチャートを参照しながら説明する。図25は、補間部226の補間処理について説明するためのフローチャートである。図25のフローチャートに沿った説明においては、補間部226を動作主体として説明する。
[Interpolation Processing]
The interpolation process by the interpolation unit 226 will be described with reference to a flowchart. Fig. 25 is a flowchart for explaining the interpolation process by the interpolation unit 226. In the explanation based on the flowchart of Fig. 25, the interpolation unit 226 will be described as the subject of operation.

図25において、まず、補間部226は、欠損区間を含む時系列データ(欠損波形)を記憶部225から取得する(ステップS231)。補間部226は、歩容情報処理部222から欠損波形を取得してもよい。欠損波形には、欠損区間の歩行フェーズを示す欠損情報が対応付けられる。25, first, the interpolation unit 226 acquires time series data including a missing section (missing waveform) from the storage unit 225 (step S231). The interpolation unit 226 may acquire the missing waveform from the gait information processing unit 222. The missing waveform is associated with missing information indicating the walking phase of the missing section.

次に、補間部226は、欠損波形に含まれる欠損区間の歩行フェーズに欠損がない複数の歩行周期の歩行波形を記憶部225から取得する(ステップS232)。Next, the interpolation unit 226 acquires from the memory unit 225 walking waveforms of multiple walking cycles in which the walking phase of the missing section included in the missing waveform is not missing (step S232).

次に、補間部226は、取得された複数の歩行波形のデータを用いて、欠損区間に対応する歩行フェーズのデータの全体分布を導出する(ステップS233)。Next, the interpolation unit 226 uses the acquired data of multiple walking waveforms to derive the overall distribution of data for the walking phase corresponding to the missing section (step S233).

次に、補間部226は、欠損区間に含まれる歩行フェーズのデータ値の全体分布における、欠損区間の起点および終点の歩行フェーズの偏差値を計算する(ステップS234)。Next, the interpolation unit 226 calculates the standard deviation of the walking phase at the start and end of the missing section in the overall distribution of data values of the walking phase included in the missing section (step S234).

次に、補間部226は、算出された偏差値を用いて、欠損区間の補間データを生成する(ステップS235)。Next, the interpolation unit 226 uses the calculated deviation value to generate interpolation data for the missing section (step S235).

次に、補間部226は、生成された補間データを欠損区間に挿入する(ステップS236)。例えば、補間部226は、補間データの起点と終点を補正して、補間データを欠損区間に挿入する。Next, the interpolation unit 226 inserts the generated interpolation data into the missing section (step S236). For example, the interpolation unit 226 corrects the start point and end point of the interpolation data and inserts the interpolation data into the missing section.

次に、補間部226は、補正された波形をリサンプリングして、時系列データに戻す(ステップS237)。ステップS237は、歩行速度や歩幅などを計算する際に、元の時系列データが必要となるために行われる。元の時系列データが必要ない場合、ステップS237は省略してもよい。Next, the interpolation unit 226 resamples the corrected waveform to return it to time series data (step S237). Step S237 is performed because the original time series data is required when calculating walking speed, stride length, etc. If the original time series data is not required, step S237 may be omitted.

以上のように、本実施形態の歩容計測システムは、計測装置、補間装置、および歩容計測装置を備える。計測装置は、ユーザの履物に配置される。計測装置は、ユーザの歩行に応じて空間加速度および空間角速度を計測する。計測装置は、計測された空間加速度および空間角速度に基づくセンサデータを生成する。計測装置は、生成されたセンサデータを補間装置に出力する。補間装置は、受信部、歩容情報処理部、欠損情報処理部、記憶部、補間部、および送信部を有する。受信部は、計測装置から送信されたセンサデータを受信する。歩容情報処理部は、センサデータの時系列データを生成する。歩容情報処理部は、足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成するとともに、時系列データにおけるデータの欠損区間を特定する。歩容情報処理部は、生成された歩行波形や、特定された欠損区間に関する情報を記憶部に記憶させる。記憶部には、歩行波形や、欠損区間に関する情報、欠損区間が補間された修復波形が記憶される。欠損情報処理部は、特定された欠損区間の歩行フェーズを計算する。補間部は、欠損区間の補間に用いられる複数の歩行周期の歩行波形のデータ群における、欠損区間に含まれる歩行フェーズごとのデータの代表値に基づいて、欠損区間の補間データを生成する。例えば、補間部は、欠損区間に含まれる歩行フェーズごとのデータの偏差値を用いて、欠損区間の補間データを生成する。補間部は、生成された補間データを欠損区間に補間する。歩容計測装置は、センサデータの時系列データを用いて補間装置が生成した欠損区間を含まない歩行波形と、補間装置が欠損区間のデータを補間した歩行波形とを補間装置から取得する。歩容計測装置は、取得した歩行波形から検出される歩行イベントに基づいてユーザの歩容を計測する。歩容計測装置は、計測されたユーザの歩容に関する情報を出力する。As described above, the gait measurement system of this embodiment includes a measurement device, an interpolation device, and a gait measurement device. The measurement device is placed on the user's footwear. The measurement device measures spatial acceleration and spatial angular velocity according to the user's gait. The measurement device generates sensor data based on the measured spatial acceleration and spatial angular velocity. The measurement device outputs the generated sensor data to the interpolation device. The interpolation device has a receiving unit, a gait information processing unit, a missing information processing unit, a storage unit, an interpolation unit, and a transmission unit. The receiving unit receives the sensor data transmitted from the measurement device. The gait information processing unit generates time-series data of the sensor data. The gait information processing unit generates a gait waveform for each gait cycle using the time-series data of the sensor data related to the movement of the foot, and identifies a missing section of data in the time-series data. The gait information processing unit stores information about the generated gait waveform and the identified missing section in the storage unit. The storage unit stores the gait waveform, information about the missing section, and a repair waveform in which the missing section is interpolated. The missing information processing unit calculates the gait phase of the identified missing section. The interpolation unit generates interpolation data for the missing section based on a representative value of data for each walking phase included in the missing section in a data group of walking waveforms of multiple walking cycles used for interpolating the missing section. For example, the interpolation unit generates the interpolation data for the missing section using a deviation value of data for each walking phase included in the missing section. The interpolation unit interpolates the generated interpolation data into the missing section. The gait measurement device acquires from the interpolation device a gait waveform that does not include the missing section generated by the interpolation device using time-series data of the sensor data, and a gait waveform in which the interpolation device has interpolated data for the missing section. The gait measurement device measures the user's gait based on walking events detected from the acquired gait waveform. The gait measurement device outputs information related to the measured gait of the user.

本実施形態の手法では、歩行の周期性に着目し、ユーザの歩行に応じて計測されるセンサデータの時系列データに発生したデータの欠損区間を補間する補間データを生成する。本実施形態の手法では、欠損区間に含まれる歩行フェーズごとのデータの代表値に基づいて、補間データを生成する。本実施形態の手法では、時系列データから生成された複数の歩行周期の歩行波形のデータの分布に基づいて、補間データを生成する。本実施形態の手法によれば、全体的な分布に基づいて補間データを生成することで、突発的な異常を含む歩行波形のデータが平滑化されるため、欠損区間のデータをより高精度に復元できる。 In the method of this embodiment, attention is focused on the periodicity of walking, and interpolation data is generated to fill in missing sections of data that occur in time-series data of sensor data measured in response to the user's walking. In the method of this embodiment, interpolation data is generated based on a representative value of data for each walking phase included in the missing section. In the method of this embodiment, interpolation data is generated based on the distribution of walking waveform data for multiple walking cycles generated from the time-series data. According to the method of this embodiment, by generating interpolation data based on the overall distribution, walking waveform data including sudden abnormalities is smoothed, so that data in the missing section can be restored with higher accuracy.

(第3の実施形態)
次に、第3の実施形態に係る補間装置について図面を参照しながら説明する。本実施形態の補間装置は、第1~第2の実施形態の補間装置を簡略化した構成である。図26は、本実施形態の補間装置32の構成の一例を示すブロック図である。補間装置32は、歩容情報処理部322、欠損情報処理部323、および補間部326を備える。
Third Embodiment
Next, an interpolation device according to a third embodiment will be described with reference to the drawings. The interpolation device of this embodiment has a simplified configuration of the interpolation devices of the first and second embodiments. Fig. 26 is a block diagram showing an example of the configuration of an interpolation device 32 of this embodiment. The interpolation device 32 includes a gait information processing unit 322, a missing information processing unit 323, and an interpolation unit 326.

歩容情報処理部322は、足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成するとともに、時系列データにおけるデータの欠損区間を特定する。欠損情報処理部323は、特定された欠損区間の歩行フェーズを計算する。補間部326は、欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、欠損区間の歩行フェーズのデータを用いて、欠損区間を補間する補間データを生成する。補間部326は、生成された補間データを欠損区間に補間する。The gait information processing unit 322 generates a walking waveform for each walking cycle using time series data of sensor data related to foot movement, and identifies missing sections of data in the time series data. The missing information processing unit 323 calculates the walking phase of the identified missing section. The interpolation unit 326 generates interpolation data that interpolates the missing section using data on the walking phase of the missing section in a walking waveform of a walking cycle different from the walking cycle that includes the missing section. The interpolation unit 326 interpolates the generated interpolation data into the missing section.

本実施形態によれば、歩行の周期性に着目して生成された補間データを用いるため、センサデータの時系列データに含まれるデータの欠損区間の特徴を含めて、欠損区間のデータを補間できる。 According to this embodiment, interpolated data generated with a focus on the periodicity of walking is used, so that data for missing sections can be interpolated, including the characteristics of the missing sections of data contained in the time series data of the sensor data.

(ハードウェア)
ここで、本開示の各実施形態に係る制御や処理を実行するハードウェア構成について、図27の情報処理装置90を一例として挙げて説明する。なお、図27の情報処理装置90は、各実施形態の制御や処理を実行するための構成例であって、本開示の範囲を限定するものではない。
(Hardware)
Here, a hardware configuration for executing the control and processing according to each embodiment of the present disclosure will be described using an information processing device 90 in Fig. 27 as an example. Note that the information processing device 90 in Fig. 27 is an example configuration for executing the control and processing according to each embodiment, and does not limit the scope of the present disclosure.

図27のように、情報処理装置90は、プロセッサ91、主記憶装置92、補助記憶装置93、入出力インターフェース95、および通信インターフェース96を備える。図27においては、インターフェースをI/F(Interface)と略記する。プロセッサ91、主記憶装置92、補助記憶装置93、入出力インターフェース95、および通信インターフェース96は、バス98を介して、互いにデータ通信可能に接続される。また、プロセッサ91、主記憶装置92、補助記憶装置93、および入出力インターフェース95は、通信インターフェース96を介して、インターネットやイントラネットなどのネットワークに接続される。As shown in FIG. 27, the information processing device 90 includes a processor 91, a main memory device 92, an auxiliary memory device 93, an input/output interface 95, and a communication interface 96. In FIG. 27, the interface is abbreviated as I/F (Interface). The processor 91, the main memory device 92, the auxiliary memory device 93, the input/output interface 95, and the communication interface 96 are connected to each other via a bus 98 so as to be able to communicate data with each other. In addition, the processor 91, the main memory device 92, the auxiliary memory device 93, and the input/output interface 95 are connected to a network such as the Internet or an intranet via the communication interface 96.

プロセッサ91は、補助記憶装置93等に格納されたプログラムを、主記憶装置92に展開する。プロセッサ91は、主記憶装置92に展開されたプログラムを実行する。本実施形態においては、情報処理装置90にインストールされたソフトウェアプログラムを用いる構成とすればよい。プロセッサ91は、本実施形態に係る制御や処理を実行する。The processor 91 expands a program stored in the auxiliary storage device 93 or the like into the main storage device 92. The processor 91 executes the program expanded into the main storage device 92. In this embodiment, a configuration may be used in which a software program installed in the information processing device 90 is used. The processor 91 executes the control and processing related to this embodiment.

主記憶装置92は、プログラムが展開される領域を有する。主記憶装置92には、プロセッサ91によって、補助記憶装置93等に格納されたプログラムが展開される。主記憶装置92は、例えばDRAM(Dynamic Random Access Memory)などの揮発性メモリによって実現される。また、主記憶装置92として、MRAM(Magnetoresistive Random Access Memory)などの不揮発性メモリが構成/追加されてもよい。The main memory 92 has an area in which programs are expanded. Programs stored in the auxiliary memory 93 or the like are expanded in the main memory 92 by the processor 91. The main memory 92 is realized by a volatile memory such as a DRAM (Dynamic Random Access Memory). In addition, a non-volatile memory such as an MRAM (Magnetoresistive Random Access Memory) may be configured/added to the main memory 92.

補助記憶装置93は、プログラムなどの種々のデータを記憶する。補助記憶装置93は、ハードディスクやフラッシュメモリなどのローカルディスクによって実現される。なお、種々のデータを主記憶装置92に記憶させる構成とし、補助記憶装置93を省略することも可能である。The auxiliary storage device 93 stores various data such as programs. The auxiliary storage device 93 is realized by a local disk such as a hard disk or flash memory. Note that it is also possible to configure the main storage device 92 to store various data and omit the auxiliary storage device 93.

入出力インターフェース95は、規格や仕様に基づいて、情報処理装置90と周辺機器とを接続するためのインターフェースである。通信インターフェース96は、規格や仕様に基づいて、インターネットやイントラネットなどのネットワークを通じて、外部のシステムや装置に接続するためのインターフェースである。入出力インターフェース95および通信インターフェース96は、外部機器と接続するインターフェースとして共通化してもよい。The input/output interface 95 is an interface for connecting the information processing device 90 to peripheral devices based on standards and specifications. The communication interface 96 is an interface for connecting to external systems and devices through a network such as the Internet or an intranet based on standards and specifications. The input/output interface 95 and the communication interface 96 may be a common interface for connecting to external devices.

情報処理装置90には、必要に応じて、キーボードやマウス、タッチパネルなどの入力機器が接続されてもよい。それらの入力機器は、情報や設定の入力に使用される。なお、タッチパネルを入力機器として用いる場合は、表示機器の表示画面が入力機器のインターフェースを兼ねる構成としてもよい。プロセッサ91と入力機器との間のデータ通信は、入出力インターフェース95に仲介させればよい。 If necessary, input devices such as a keyboard, mouse, or touch panel may be connected to the information processing device 90. These input devices are used to input information and settings. When a touch panel is used as an input device, the display screen of the display device may also serve as an interface for the input device. Data communication between the processor 91 and the input devices may be mediated by the input/output interface 95.

また、情報処理装置90には、情報を表示するための表示機器を備え付けてもよい。表示機器を備え付ける場合、情報処理装置90には、表示機器の表示を制御するための表示制御装置(図示しない)が備えられていることが好ましい。表示機器は、入出力インターフェース95を介して情報処理装置90に接続すればよい。The information processing device 90 may also be equipped with a display device for displaying information. When a display device is equipped, it is preferable that the information processing device 90 is equipped with a display control device (not shown) for controlling the display of the display device. The display device may be connected to the information processing device 90 via the input/output interface 95.

また、情報処理装置90には、ドライブ装置が備え付けられてもよい。ドライブ装置は、プロセッサ91と記録媒体(プログラム記録媒体)との間で、記録媒体からのデータやプログラムの読み込み、情報処理装置90の処理結果の記録媒体への書き込みなどを仲介する。ドライブ装置は、入出力インターフェース95を介して情報処理装置90に接続すればよい。The information processing device 90 may also be equipped with a drive device. The drive device mediates between the processor 91 and a recording medium (program recording medium), reading data and programs from the recording medium, writing the processing results of the information processing device 90 to the recording medium, and the like. The drive device may be connected to the information processing device 90 via the input/output interface 95.

以上が、本発明の各実施形態に係る制御や処理を可能とするためのハードウェア構成の一例である。なお、図27のハードウェア構成は、各実施形態に係る制御や処理を実行するためのハードウェア構成の一例であって、本発明の範囲を限定するものではない。また、各実施形態に係る制御や処理をコンピュータに実行させるプログラムも本発明の範囲に含まれる。さらに、各実施形態に係るプログラムを記録したプログラム記録媒体も本発明の範囲に含まれる。記録媒体は、例えば、CD(Compact Disc)やDVD(Digital Versatile Disc)などの光学記録媒体で実現できる。記録媒体は、USB(Universal Serial Bus)メモリやSD(Secure Digital)カードなどの半導体記録媒体によって実現されてもよい。また、記録媒体は、フレキシブルディスクなどの磁気記録媒体、その他の記録媒体によって実現されてもよい。プロセッサが実行するプログラムが記録媒体に記録されている場合、その記録媒体はプログラム記録媒体に相当する。The above is an example of a hardware configuration for enabling the control and processing according to each embodiment of the present invention. The hardware configuration in FIG. 27 is an example of a hardware configuration for executing the control and processing according to each embodiment, and does not limit the scope of the present invention. In addition, a program that causes a computer to execute the control and processing according to each embodiment is also included in the scope of the present invention. Furthermore, a program recording medium on which a program according to each embodiment is recorded is also included in the scope of the present invention. The recording medium can be realized, for example, by an optical recording medium such as a CD (Compact Disc) or a DVD (Digital Versatile Disc). The recording medium may be realized by a semiconductor recording medium such as a USB (Universal Serial Bus) memory or an SD (Secure Digital) card. The recording medium may also be realized by a magnetic recording medium such as a flexible disk or other recording medium. When a program executed by a processor is recorded on a recording medium, the recording medium corresponds to a program recording medium.

各実施形態の構成要素は、任意に組み合わせてもよい。また、各実施形態の構成要素は、ソフトウェアによって実現されてもよいし、回路によって実現されてもよい。The components of each embodiment may be combined in any manner. The components of each embodiment may be realized by software or by circuitry.

以上、実施形態を参照して本発明を説明してきたが、本発明は上記実施形態に限定されるものではない。本発明の構成や詳細には、本発明のスコープ内で当業者が理解し得る様々な変更をすることができる。Although the present invention has been described above with reference to the embodiments, the present invention is not limited to the above-mentioned embodiments. Various modifications that can be understood by a person skilled in the art can be made to the configuration and details of the present invention within the scope of the present invention.

1、2 歩容計測システム
11、21 計測装置
12、22、32 補間装置
13、23 歩容計測装置
111 加速度センサ
112 角速度センサ
113 制御部
115 送信部
121、221 受信部
122、222、322 歩容情報処理部
123、223、323 欠損情報処理部
125、225 記憶部
126、226、326 補間部
128、228 送信部
131 取得部
132 検出部
133 歩容計測部
REFERENCE SIGNS LIST 1, 2 Gait measurement system 11, 21 Measurement device 12, 22, 32 Interpolation device 13, 23 Gait measurement device 111 Acceleration sensor 112 Angular velocity sensor 113 Control unit 115 Transmission unit 121, 221 Reception unit 122, 222, 322 Gait information processing unit 123, 223, 323 Missing information processing unit 125, 225 Storage unit 126, 226, 326 Interpolation unit 128, 228 Transmission unit 131 Acquisition unit 132 Detection unit 133 Gait measurement unit

Claims (9)

足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成するとともに、前記時系列データにおけるデータの欠損区間を特定する歩行情報処理手段と、
特定された前記欠損区間の歩行フェーズを計算する欠損情報処理手段と、
前記欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、前記欠損区間の前記歩行フェーズのデータを用いて、前記欠損区間を補間する補間データを生成し、生成された前記補間データを前記欠損区間に補間する補間手段と、を備え
前記補間手段は、
前記時系列データに含まれる前記欠損区間に関して、前記欠損区間に含まれる前記歩行フェーズのデータが欠損していない少なくとも一つの前記歩行波形から、前記欠損区間に含まれる前記歩行フェーズのデータを少なくとも一つ選択し、
選択された前記歩行フェーズのデータを用いて、前記補間データを生成する補間装置。
a gait information processing means for generating a gait waveform for each gait cycle using time series data of the sensor data relating to foot movement and for identifying a data missing section in the time series data;
A missing information processing means for calculating a walking phase of the identified missing section;
an interpolation means for generating interpolation data for interpolating the missing section using data of the walking phase of the missing section in a walking waveform having a walking cycle different from a walking cycle including the missing section, and for interpolating the generated interpolation data into the missing section ;
The interpolation means
With respect to the missing section included in the time series data, at least one piece of data of the walking phase included in the missing section is selected from at least one of the walking waveforms in which data of the walking phase included in the missing section is not missing;
An interpolation device that generates the interpolated data using data of the selected walking phase .
前記歩行情報処理手段は、
前記時系列データから一歩行周期分の波形を切り出し、
切り出された前記波形に含まれる前記歩行フェーズを正規化して前記歩行波形を生成する請求項1に記載の補間装置。
The walking information processing means includes:
A waveform for one walking cycle is extracted from the time series data,
The interpolation device according to claim 1 , wherein the walking waveform is generated by normalizing the walking phase included in the extracted waveform.
前記補間手段は、
複数の前記歩行周期の前記歩行波形から選択された前記歩行フェーズのデータの平均値を計算し、
算出された前記歩行フェーズのデータの平均値を用いて、前記補間データを生成する請求項1または2に記載の補間装置。
The interpolation means
Calculating an average value of data of the walking phase selected from the walking waveforms of the plurality of walking cycles;
3. The interpolation device according to claim 1, wherein the interpolated data is generated using an average value of the calculated data of the walking phase.
前記補間手段は、
遊脚相と比べてデータの変動が小さい立脚相に含まれる歩行ピリオドに関しては、解析的な手法で前記補間データを生成し、
前記立脚相と比べてデータの変動が大きい前記遊脚相に含まれる歩行ピリオドに関しては、前記欠損区間を含む前記歩行周期とは異なる前記歩行周期のデータを用いて前記補間データを生成する請求項1乃至のいずれか一項に記載の補間装置。
The interpolation means
For a walking period included in the stance phase in which data fluctuation is smaller than that of the swing phase , the interpolated data is generated using an analytical method;
4. The interpolation device according to claim 1, wherein, for a walking period included in the swing phase in which data fluctuation is larger than that of the stance phase, the interpolated data is generated using data of a walking cycle different from the walking cycle including the missing section.
前記補間手段は、
前記欠損区間の補間に用いられる複数の前記歩行周期の前記歩行波形のデータ群における、前記欠損区間に含まれる前記歩行フェーズごとのデータの代表値に基づいて、前記欠損区間の前記補間データを生成する請求項1乃至のいずれか一項に記載の補間装置。
The interpolation means
5. The interpolation device according to claim 1, wherein the interpolation data for the missing section is generated based on a representative value of data for each of the walking phases included in the missing section in a data group of the walking waveform for a plurality of the walking cycles used for interpolating the missing section.
前記補間手段は、
前記欠損区間に含まれる前記歩行フェーズごとのデータの偏差値を用いて、前記欠損区間の前記補間データを生成する請求項5に記載の補間装置。
The interpolation means
The interpolation device according to claim 5 , wherein the interpolation data for the missing section is generated using a deviation value of data for each of the walking phases included in the missing section.
請求項1乃至のいずれか一項に記載の補間装置と、
ユーザの履物に配置され、前記ユーザの歩行に応じて空間加速度および空間角速度を計測し、計測された前記空間加速度および前記空間角速度に基づくセンサデータを生成し、生成された前記センサデータを前記補間装置に出力する計測装置と、
前記センサデータの時系列データを用いて前記補間装置が生成した欠損区間を含まない歩行波形と、前記補間装置が欠損区間のデータを補間した歩行波形とを前記補間装置から取得し、取得した前記歩行波形から検出される歩行イベントに基づいて前記ユーザの歩容を計測し、計測された前記ユーザの歩容に関する情報を出力する歩容計測装置と、を備える歩容計測システム。
An interpolation device according to any one of claims 1 to 6 ,
a measurement device that is arranged in a user's footwear, measures spatial acceleration and spatial angular velocity in response to the user's walking, generates sensor data based on the measured spatial acceleration and spatial angular velocity, and outputs the generated sensor data to the interpolation device;
a gait measurement device that acquires, from the interpolation device, a gait waveform that does not include missing sections that is generated by the interpolation device using time-series data of the sensor data and a gait waveform in which the interpolation device interpolates data of the missing sections, measures a gait of the user based on walking events detected from the acquired gait waveform, and outputs information regarding the measured gait of the user.
コンピュータが、
足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成し、
前記時系列データにおけるデータの欠損区間を特定し、
特定された前記欠損区間の歩行フェーズを計算し、
前記欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、前記欠損区間の前記歩行フェーズのデータを用いて、前記欠損区間を補間する補間データを生成し、
生成された前記補間データを前記欠損区間に補間し、
前記補間において、
前記時系列データに含まれる前記欠損区間に関して、前記欠損区間に含まれる前記歩行フェーズのデータが欠損していない少なくとも一つの前記歩行波形から、前記欠損区間に含まれる前記歩行フェーズのデータを少なくとも一つ選択し、
選択された前記歩行フェーズのデータを用いて、前記補間データを生成する補間方法。
The computer
A walking waveform is generated for each walking cycle using time series data of the sensor data regarding the movement of the feet.
Identifying a data missing section in the time series data;
Calculate the walking phase of the identified missing section;
generating interpolation data for interpolating the missing section using data of the walking phase of the missing section in a walking waveform having a walking cycle different from the walking cycle including the missing section;
The generated interpolation data is interpolated into the missing section ;
In the interpolation,
With respect to the missing section included in the time series data, at least one piece of data of the walking phase included in the missing section is selected from at least one of the walking waveforms in which data of the walking phase included in the missing section is not missing;
An interpolation method for generating the interpolated data using data of the selected walking phase .
足の動きに関するセンサデータの時系列データを用いて歩行周期ごとの歩行波形を生成する処理と、
前記時系列データにおけるデータの欠損区間を特定する処理と、
特定された前記欠損区間の歩行フェーズを計算する処理と、
前記欠損区間を含む歩行周期とは異なる歩行周期の歩行波形における、前記欠損区間の前記歩行フェーズのデータを用いて、前記欠損区間を補間する補間データを生成する処理と、
生成された前記補間データを前記欠損区間に補間する処理と、
前記補間する処理において、
前記時系列データに含まれる前記欠損区間に関して、前記欠損区間に含まれる前記歩行フェーズのデータが欠損していない少なくとも一つの前記歩行波形から、前記欠損区間に含まれる前記歩行フェーズのデータを少なくとも一つ選択する処理と、
選択された前記歩行フェーズのデータを用いて、前記補間データを生成する処理と、をコンピュータに実行させるプログラム。
A process of generating a walking waveform for each walking cycle using time series data of the sensor data related to foot movement;
A process of identifying a data missing section in the time series data;
A process of calculating a walking phase of the identified missing section;
A process of generating interpolation data for interpolating the missing section using data of the walking phase of the missing section in a walking waveform having a walking cycle different from a walking cycle including the missing section;
A process of interpolating the generated interpolation data into the missing section;
In the interpolation process,
selecting at least one piece of data of the walking phase included in the missing section from at least one of the walking waveforms having no missing data of the walking phase included in the missing section, for the missing section included in the time series data;
and generating the interpolated data using data of the selected walking phase .
JP2023529225A 2021-06-21 2021-06-21 Interpolation device, gait measurement system, interpolation method, and program Active JP7626221B2 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/023444 WO2022269698A1 (en) 2021-06-21 2021-06-21 Interpolation device, gait measurement system, interpolation method, and recording medium

Publications (3)

Publication Number Publication Date
JPWO2022269698A1 JPWO2022269698A1 (en) 2022-12-29
JPWO2022269698A5 JPWO2022269698A5 (en) 2024-03-15
JP7626221B2 true JP7626221B2 (en) 2025-02-04

Family

ID=84545273

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2023529225A Active JP7626221B2 (en) 2021-06-21 2021-06-21 Interpolation device, gait measurement system, interpolation method, and program

Country Status (3)

Country Link
US (1) US20250134412A1 (en)
JP (1) JP7626221B2 (en)
WO (1) WO2022269698A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117021106A (en) * 2023-08-30 2023-11-10 深圳鹏行智能研究有限公司 Robot gait control method and robot

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016137228A (en) 2015-01-23 2016-08-04 村田機械株式会社 Walking measurement system
JP2019045967A (en) 2017-08-30 2019-03-22 富士通株式会社 Posture estimation device, method, and program
JP2020130335A (en) 2019-02-14 2020-08-31 日本電信電話株式会社 Time feature value calculation device, calculation method and its program

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016137228A (en) 2015-01-23 2016-08-04 村田機械株式会社 Walking measurement system
JP2019045967A (en) 2017-08-30 2019-03-22 富士通株式会社 Posture estimation device, method, and program
JP2020130335A (en) 2019-02-14 2020-08-31 日本電信電話株式会社 Time feature value calculation device, calculation method and its program

Also Published As

Publication number Publication date
JPWO2022269698A1 (en) 2022-12-29
US20250134412A1 (en) 2025-05-01
WO2022269698A1 (en) 2022-12-29

Similar Documents

Publication Publication Date Title
JP7327516B2 (en) Abnormality detection device, judgment system, abnormality detection method, and program
US20260088176A1 (en) Estimation device, estimation system, estimation method, and recording medium
JP7632676B2 (en) Gait measurement device, gait measurement system, gait measurement method, and program
JP7480851B2 (en) DETECTION APPARATUS, DETECTION SYSTEM, DETECTION METHOD, AND PROGRAM
JP7626221B2 (en) Interpolation device, gait measurement system, interpolation method, and program
US20240172966A1 (en) Harmonic index estimation device, estimation system, harmonic index estimation method, and recording medium
US20240277259A1 (en) Gait index calculation device, gait measurement system, gait index calculation method, and recording medium
JP7729406B2 (en) Dynamic balance estimation device, dynamic balance estimation system, dynamic balance estimation method, and program
US20250045354A1 (en) Feature quantity selection device, feature quantity selection method, and recording medium
JP7582318B2 (en) Estimation device, estimation system, estimation method, and program
JP7459965B2 (en) Discrimination device, discrimination system, discrimination method, and program
JP7704216B2 (en) Feature amount data generating device, gait measuring device, physical condition estimation system, feature amount data generating method, and program
JP7772091B2 (en) Muscle strength index estimation device, muscle strength index estimation system, muscle strength index estimation method, and program
JP7525052B2 (en) FEATURE GENERATION DEVICE, GATHERING MEASUREMENT SYSTEM, FEATURE GENERATION METHOD, AND PROGRAM
US20250114018A1 (en) Detection device, detection system, gait measurement system, detection method, and recording medium
JP7790451B2 (en) MOBILITY ABILITY ESTIMATION DEVICE, MOBILITY ABILITY ESTIMATION SYSTEM, MOBILITY ABILITY ESTIMATION METHOD, AND PROGRAM
JP7509229B2 (en) DETECTION APPARATUS, DETECTION SYSTEM, DETECTION METHOD, AND PROGRAM
US20240315600A1 (en) Gait information generation device, gait measurement system, gait information generation method, and recording medium
US20230397839A1 (en) Waist swinging estimation device, estimation system, waist swinging estimation method, and recording medium
JP7715212B2 (en) Static balance estimation device, static balance estimation system, static balance estimation method, and program
US20240138757A1 (en) Pelvic inclination estimation device, estimation system, pelvic inclination estimation method, and recording medium
US20240161921A1 (en) Biometric information processing device, information processing system, biometric information processing method, and recording medium
JP7494941B2 (en) Gait index calculation device, gait index calculation system, gait index calculation method, and program
JP2024140356A (en) KNEE PAIN RISK ESTIMATION DEVICE, PHYSICAL CONDITION ESTIMATION SYSTEM, KNEE PAIN RISK ESTIMATION METHOD, AND PROGRAM

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231213

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231213

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20241008

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20241206

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20241224

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20250106

R150 Certificate of patent or registration of utility model

Ref document number: 7626221

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150