[go: up one dir, main page]

WO2012036135A1 - Procédé de traitement de l'information, dispositif de traitement de l'information, dispositif de sortie, système de traitement de l'information, programme de traitement de l'information et support d'enregistrement lisible par ordinateur sur lequel le même programme est enregistré - Google Patents

Procédé de traitement de l'information, dispositif de traitement de l'information, dispositif de sortie, système de traitement de l'information, programme de traitement de l'information et support d'enregistrement lisible par ordinateur sur lequel le même programme est enregistré Download PDF

Info

Publication number
WO2012036135A1
WO2012036135A1 PCT/JP2011/070763 JP2011070763W WO2012036135A1 WO 2012036135 A1 WO2012036135 A1 WO 2012036135A1 JP 2011070763 W JP2011070763 W JP 2011070763W WO 2012036135 A1 WO2012036135 A1 WO 2012036135A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
phase
signal
evaluation coefficient
predetermined
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2011/070763
Other languages
English (en)
Japanese (ja)
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.)
Mitsubishi Chemical Corp
Original Assignee
Mitsubishi Chemical 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 Mitsubishi Chemical Corp filed Critical Mitsubishi Chemical Corp
Priority to JP2012534001A priority Critical patent/JPWO2012036135A1/ja
Publication of WO2012036135A1 publication Critical patent/WO2012036135A1/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4082—Diagnosing or monitoring movement diseases, e.g. Parkinson, Huntington or Tourette
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Measuring devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor or mobility of a limb
    • A61B5/112—Gait analysis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition

Definitions

  • the present invention relates to an information processing method, an information processing device, an output device, an information processing system, an information processing program, and a computer-readable recording medium on which the program is recorded.
  • the synchrony is an index indicating whether or not the ratio of the time from the right foot landing to the left foot landing and the time from the left foot landing to the right foot landing maintains a constant value (for example, 1) during walking. . That is, if the ratio of the time from the right foot landing to the left foot landing and the time from the left foot landing to the right foot landing is maintained at a constant value (1), it is judged that the synchrony is high, and the time from the right foot landing to the left foot landing is If the ratio from the left foot landing to the right foot landing does not maintain a constant value (1), it is determined that the synchrony is low.
  • a constant value for example, 1
  • the synchronism of walking is determined by determining the ratio of the time required for the rhythm representing walking to change by one cycle to the time required for the rhythm representing walking to change for the next one cycle. Conventionally, this ratio has been obtained from the distance between the feature points of the heel-off or toe-off.
  • the synchrony is evaluated using the interval between the feature points. For this reason, when the peak of a waveform is used as a feature point, a large error occurs if the above ratio is calculated from a waveform that does not have a clear peak, such as a split of the peak. That is, the evaluation of synchrony using the interval between feature points is not suitable for a waveform having no clear peak. Further, when the feature point interval is used, the ratio can be calculated only for each time corresponding to the peak position of the waveform. That is, the evaluation of synchrony using the interval between feature points is not suitable for continuous processing or real-time processing.
  • the present invention has been made in view of such problems, and an object thereof is to continuously evaluate the synchrony of body movement rhythm.
  • the present invention is not limited to the above-described object, and other effects of the present invention can be achieved by the functions and effects derived from the respective configurations shown in the embodiments for carrying out the invention which will be described later. Can be positioned as one of
  • the gist of the present invention is that a signal representing a repetitive voluntary movement of a living body, a time when the phase of the signal changes by a predetermined angle from a phase of the signal at a predetermined time (hereinafter referred to as a first time), and a predetermined A time calculating process for calculating a time (hereinafter referred to as a second time) when the phase of the signal has changed twice the predetermined angle from the phase of the signal at the time, and the first time calculated by the time calculating process And an evaluation coefficient calculation step of calculating an evaluation coefficient for evaluating the synchrony of the repetitive voluntary movement based on the second time.
  • the gist of the present invention is that a signal representing a repetitive voluntary movement of a living body, a time when the phase of the signal changes by a predetermined angle from a phase of the signal at a predetermined time (hereinafter referred to as a first time), and a predetermined A time calculating unit that calculates a time when the phase of the signal changes twice the predetermined angle from the phase of the signal at a time (hereinafter referred to as a second time), and the first time calculated by the time calculating unit. And an evaluation coefficient calculation unit that calculates an evaluation coefficient for evaluating the synchronization of the repetitive voluntary movement based on the second time.
  • the gist of the present invention is that the phase of the signal changes by a predetermined angle from the phase of the signal representing the repetitive voluntary movement of the living body at the predetermined time, and the phase of the signal from the phase of the signal at the predetermined time.
  • An output device including an output unit that outputs an evaluation coefficient for evaluating the synchronization of the repeated voluntary movement based on the time when the predetermined angle has been changed twice in a form in which the evaluation coefficient can be identified.
  • the gist of the present invention is that a signal representing a repetitive voluntary movement of a living body, a time when the phase of the signal changes by a predetermined angle from a phase of the signal at a predetermined time (hereinafter referred to as a first time), and a predetermined A time calculating unit that calculates a time when the phase of the signal changes twice the predetermined angle from the phase of the signal at a time (hereinafter referred to as a second time), and the first time calculated by the time calculating unit.
  • an evaluation coefficient calculating unit that calculates an evaluation coefficient for evaluating the synchrony of the repetitive voluntary movement based on the second time, and the evaluation coefficient calculated by the evaluation coefficient calculating unit is identified as the evaluation coefficient
  • the gist of the present invention is that a signal representing a repetitive voluntary movement of a living body, a time when the phase of the signal changes by a predetermined angle from a phase of the signal at a predetermined time (hereinafter referred to as a first time), and a predetermined A time calculating process for calculating a time (hereinafter referred to as a second time) when the phase of the signal has changed twice the predetermined angle from the phase of the signal at the time, and the first time calculated by the time calculating process And an evaluation coefficient calculation process for calculating an evaluation coefficient for evaluating the synchrony of the repetitive voluntary movement based on the second time.
  • the gist of the present invention is that a signal representing a repetitive voluntary movement of a living body, a time when the phase of the signal changes by a predetermined angle from a phase of the signal at a predetermined time (hereinafter referred to as a first time), and a predetermined A time calculating process for calculating a time (hereinafter referred to as a second time) when the phase of the signal has changed twice the predetermined angle from the phase of the signal at the time, and the first time calculated by the time calculating process
  • a computer-readable recording medium recording an information processing program for causing a computer to execute an evaluation coefficient calculation process for calculating an evaluation coefficient for evaluating the synchronization of the repetitive voluntary movement based on the second time Exist.
  • FIG. 1 is a diagram illustrating a configuration of a system as an example of an embodiment.
  • FIG. 2 is a diagram illustrating an exercise rhythm as an example of the embodiment.
  • FIG. 3 is a diagram for explaining the operation of the determination unit as an example of the embodiment.
  • FIG. 4 is a flowchart for explaining the operation of the system as an example of the embodiment.
  • 5 (A) and 5 (B) are diagrams showing body motion signals and exercise rhythms as an example of the embodiment, and FIG. 5 (A) shows a body motion signal detection device at the center of the abdomen of a healthy subject.
  • FIG. 5B shows an example of an exercise rhythm extracted from the body motion signal.
  • FIG. 5B shows a part of the 15-minute walking acceleration signal measured by the body motion signal detection device at 100 Hz sampling.
  • FIG. 6 is a flowchart for explaining the operation of the time constant determination unit as an example of the embodiment.
  • FIG. 7 is a diagram illustrating a change in dispersion with respect to a change in time constant as an example of an embodiment.
  • FIG. 8 is a flowchart for explaining the operation of the rhythm extraction unit as an example of the embodiment.
  • FIG. 9 is a flowchart for explaining the operation of the evaluation coefficient determination unit as an example of the embodiment.
  • FIGS. 10A to 10C are diagrams for explaining a pattern matching method as an example of the embodiment.
  • FIG. 10A shows a part of the waveform from the data of FIG. FIG.
  • FIG. 10 (B) shows the result of calculating the autocorrelation coefficient by selecting a reference wave with a width of 0.4 seconds around the time indicated by * in FIG. 10 (A).
  • FIG. 10C shows a result of calculating an autocorrelation coefficient by selecting a reference wave having a width of 0.4 seconds around the time indicated by a circle in FIG. 10A.
  • FIG. 11 is a diagram illustrating evaluation coefficients as an example of the embodiment.
  • FIG. 12 is a flowchart for explaining the operation of the determination unit as an example of the embodiment.
  • FIG. 13 is a flowchart for explaining the operation of the determination unit as an example of the embodiment.
  • FIG. 14 is a diagram illustrating evaluation coefficients as an example of the embodiment.
  • FIG. 15 is a diagram illustrating an exercise rhythm as an example of the embodiment.
  • FIGS. 16A to 16C are diagrams for explaining the evaluation coefficient obtained by the system as an example of the embodiment.
  • FIG. 16A is detected by the body motion signal detection device 10.
  • FIG. 16B is a diagram showing an acceleration signal after second order integration, that is, a motion trajectory
  • FIG. 16C is an evaluation obtained from the acceleration signal after second order integration. It is a figure which shows a coefficient.
  • FIG. 17A to FIG. 17C are diagrams for explaining evaluation coefficients obtained by the system as an example of the embodiment.
  • FIG. 17A is detected by the body motion signal detection device 10.
  • FIG. 17B is a diagram showing an acceleration signal after second-order integration, that is, a motion trajectory, and FIG.
  • FIG. 17C is an evaluation obtained from the acceleration signal after second-order integration. It is a figure which shows a coefficient.
  • 18A to 18C are diagrams for explaining the evaluation coefficient obtained by the system as an example of the embodiment.
  • FIG. 18A is detected by the body motion signal detection device 10.
  • 18 (B) is a diagram showing the acceleration signal and angular velocity signal after the second integration
  • FIG. 18 (C) is a graph showing the acceleration signal and angular velocity after the second integration. It is a figure which shows the evaluation coefficient each calculated
  • FIG. 19A is a diagram showing the time change of the walking index after the smoothing process is performed
  • FIG. 19B is a diagram showing the time change of the body motion index after the smoothing process is performed. It is.
  • FIG. 1 is a diagram illustrating a configuration of a system according to an example of an embodiment.
  • the system 1 includes a body motion signal detection device 10 and an information processing device 20.
  • the body motion signal detection device 10 is communicably connected to the information processing device 20 via, for example, wired such as the Internet or wireless such as a wireless LAN (Local Area Network) or Bluetooth (registered trademark).
  • wired such as the Internet
  • wireless such as a wireless LAN (Local Area Network) or Bluetooth (registered trademark).
  • the body motion signal detection device 10 detects (measures) a repetitive rhythmic motion accompanying a repetitive voluntary movement of a subject (living body) (hereinafter sometimes simply referred to as a voluntary movement) noninvasively and continuously, for example.
  • a repetitive rhythmic motion accompanying a repetitive voluntary movement of a subject (hereinafter sometimes simply referred to as a voluntary movement) noninvasively and continuously, for example.
  • the repeated rhythmic movement accompanying the voluntary movement may be simply referred to as rhythmic movement.
  • the voluntary exercise is, for example, walking, jogging, running, cycling, swimming, gymnastics, weight training, physical strength measurement (stepping up / down, repetitive side jumping), juggling (beading a beanbag, lifting a soccer ball) and the like.
  • the repetitive rhythmic movement accompanying the voluntary movement includes, for example, the rhythmic movement of the walking itself when the voluntary movement is walking.
  • the voluntary movement is a movement that is repeated regularly, the rhythm of the voluntary movement itself can be extracted more accurately.
  • movement repeated regularly includes not only that the completely same exercise
  • Non-invasive means, for example, that the body of the subject is not damaged or that the subject is not burdened.
  • the body motion signal detection device 10 is configured to be portable, for example. In addition, if the attachment position to the test subject of the body motion signal detection apparatus 10 is a site
  • the body motion signal detection device 10 includes, for example, a body motion signal detection unit 11, a storage unit 12, and an interface unit 13. The body motion signal detection unit 11, the storage unit 12, and the interface unit 13 are connected to be communicable with each other.
  • the body motion signal detection unit 11 detects (measures), for example, a repetitive rhythm motion associated with a voluntary motion as a body motion signal (a signal based on a repetitive rhythm motion associated with a voluntary motion of a living body). That is, the body motion signal detection unit 11 detects a signal based on a repetitive rhythm movement accompanying a voluntary movement of a living body as a body motion signal. From a different point of view, one body motion signal detection unit 11 detects, for example, a repetitive rhythm motion accompanying a voluntary motion of a living body as a body motion signal.
  • the subject's rhythmic movements for example, force changes, spatial body position changes, sounds emitted from the body, waves such as electromagnetic waves or fine energy changes, or field changes around the body, etc. Detect (measure) as a body motion signal. That is, the body motion signal detection unit 11 measures a signal based on repetitive rhythm movement accompanying voluntary movement.
  • the body motion signal detection unit 11 is realized by an inertial sensor such as an acceleration sensor, a speed sensor, or a gyro sensor, for example.
  • an inertial sensor such as an acceleration sensor, a speed sensor, or a gyro sensor, for example.
  • the inertial sensor which detects the said body motion signal it selects suitably according to the kind of signal to detect, for example.
  • an acceleration sensor that measures acceleration of body movement is preferably used, but is not limited to an acceleration sensor.
  • the acceleration sensor a single-axis to three-axis sensor can be arbitrarily used, but the three axes for detecting acceleration acting in three directions of the vertical direction, the horizontal front-rear direction, and the horizontal left-right direction during walking.
  • an acceleration sensor is preferably used, it is not limited to a triaxial acceleration sensor.
  • the body motion signal detection unit 11 measures the body motion signal at a predetermined sampling frequency (for example, 100 Hz), for example.
  • the storage unit 12 is a storage device capable of storing various information such as RAM (Random Access Memory), HDD (Hard Disk Drive), SSD (Solid State Drive), flash memory, and the like. Specifically, the storage unit 12 stores the body motion signal obtained by the body motion signal detection unit 11, for example. Moreover, the memory
  • the interface unit 13 is an interface connected to the information processing apparatus 20 so as to be communicable, for example.
  • the interface unit 13 includes an antenna.
  • the interface unit 13 is a connection terminal that can be connected to a wire.
  • the storage unit 12 is detachably provided, for example, in a slot provided in the body motion signal detection device 10, and the storage unit 12 is detached from the body motion signal detection device 10 and connected to the information processing device 20. In this case, the interface unit 13 may not be provided in the body motion signal detection device 10.
  • the information processing apparatus 20 is, for example, a PC (Personal Computer), and calculates an evaluation coefficient to be described later from the body motion signal obtained by the body motion signal detection apparatus 10.
  • the information processing apparatus 20 includes, for example, a central processing unit 21, a storage unit 22, an output unit 23, and an interface unit 24.
  • the central processing unit 21, the storage unit 22, the output unit 23, and the interface unit 24 are connected to be communicable with each other.
  • the storage unit 22 is a storage device that can store application programs and data, such as RAM, HDD, and SSD, and stores various types of information.
  • the central processing unit 21 is a processing device that performs various calculations or controls by executing various application programs stored in the storage unit 22, thereby realizing various functions.
  • the central processing unit 21 executes an information processing program stored in the storage unit 22 to thereby execute a time constant determination unit 211, a rhythm extraction unit 212, an evaluation coefficient determination unit 213, a determination unit 214, and an output control unit 215.
  • the information processing program is a program that causes the central processing unit 21 to function as the time constant determination unit 211, the rhythm extraction unit 212, the evaluation coefficient determination unit 213, the determination unit 214, and the output control unit 215.
  • the time constant determination unit 211 acquires, for example, the body motion signal obtained by the body motion signal detection unit 11 from the storage unit 12, and based on the acquired body motion signal, a signal (hereinafter, referred to as a rhythm of voluntary exercise itself).
  • the time constant used to extract the motion rhythm is sometimes calculated.
  • the movement rhythm corresponds to a signal representing repeated voluntary movement of the living body.
  • the body motion signal including the motion rhythm is a signal detected by one body motion signal detection unit 11, a signal representing repeated voluntary movement of the living body is generated by one body motion signal detection unit 11 (detection unit). It is a detected signal.
  • the time constant determination unit 211 determines the time constant by performing the following processing.
  • the time constant determining unit 211 applies the body motion signal X to a high-pass filter or a band-pass filter characterized by a certain time constant A.
  • a high-pass filter when the process of smoothing the body motion signal X with a zero phase moving average filter having a time width (time constant) A is described as F (X, A), the output waveform (signal)
  • F (X, A) the process of smoothing the body motion signal X with a zero phase moving average filter having a time width (time constant) A
  • F (X, A) the output waveform
  • a / 2.5 is performed.
  • the zero phase moving average filter refers to a moving average filter whose phase shift is zero.
  • the zero phase moving average filter can be realized by using various known methods, and detailed description thereof is omitted.
  • An example of the time constant of the bandpass filter is A / 2.5, but the present invention is not limited to this.
  • the time constant determination unit 211 quantifies the regularity of the output waveform Y.
  • the time constant determination unit 211 obtains regularity of the output waveform Y when the time constant A is changed.
  • the time constant determination unit 211 graphs the regularity of the output waveform Y when the time constant T is changed, for example.
  • the time constant determination unit 211 calculates a change in the regularity of the output waveform Y with respect to the time constant A by, for example, calculating the CV at each time constant A when the time constant A is changed.
  • the example using the absolute value CV at the peak of the output waveform Y has been described as the quantification of the regularity of the output waveform Y in the time constant determination unit 211, it is not limited to this example.
  • the time constant determination unit 211 obtains a minimum point from the change in regularity (change in CV) of the output waveform Y with respect to the time constant A obtained in the process of (3) above. Then, the time constant determination unit 211 separates the motion rhythm from the body motion signal, for example, the time constant A corresponding to the minimum point with the smaller time constant A among the obtained minimum points (for example, two minimum points). Thus, it is determined as a time constant for extraction (hereinafter sometimes referred to as a first time constant). That is, the time constant determination unit 211 selects a time constant that enables separation and extraction of a waveform having regularity as much as possible.
  • the time constant determination unit 211 selects a plurality of minimum points (for example, two points) and determines the time constant corresponding to the minimum point with the smaller time constant A as the first time constant. That is, the time constant determination unit 211 determines the first time constant based on the output when the band pass filter is applied to the body motion signal.
  • N 1, 2, 3,
  • the order of the integration operation and the filtering process may be any order as long as the following conditions are satisfied, for example. ⁇ Do not perform the filter process more than once in succession. ⁇ At least once in the filter process, the integration process is completed. -When the filter process is performed twice or more, the filter is a high-pass filter except for the last one.
  • the rhythm extraction unit 212 acquires the body motion signal obtained by the body motion signal detection unit 11 from the storage unit 12, and performs filtering using the first time constant (A1) determined by the time constant determination unit 211. By processing, the movement rhythm is extracted from the body movement signal. Specifically, the rhythm extraction unit 212 extracts an exercise rhythm from the body motion signal by performing the following process (filtering process), for example. In addition, the following process is an illustration and is not limited to this example.
  • the rhythm extraction unit 212 performs integration on the signal after the high-pass filtering.
  • the rhythm extraction unit 212 performs second-order integration on the signal after the high-pass filtering.
  • the rhythm extraction unit 212 performs the same processing as (10) on the signal X2 after integration in the processing (11). That is, the rhythm extraction unit 212 performs processing represented by F (X2, A1) on the integrated body motion signal. In other words, the rhythm extraction unit 212 performs high-pass filtering on the integrated signal that is the signal after integration, using the first time constant. (8) Furthermore, the rhythm extraction unit 212 creates an envelope connecting the maximum values of the signal obtained by the processing of (12) and an envelope connecting the minimum values, and the two envelopes A process of subtracting the average signal from the signal obtained in the process (12) is performed. Since the amplitude of the signal obtained by this processing changes with the amplitude value 0 interposed therebetween, the phase can be accurately obtained when the Hilbert transform method described later is used. This process may be omitted.
  • the rhythm extraction unit 212 extracts an exercise rhythm from the body motion signal. That is, the rhythm extraction unit 212 functions as a rhythm extraction unit that extracts a rhythm (voluntary rhythm) of voluntary exercise from the body motion signal by performing filtering processing on the body motion signal. Therefore, for example, when evaluating the synchrony of the walking rhythm from the body motion signal, the walking area may be extracted in advance, or such preprocessing may not be performed.
  • the body motion signal is second-order integrated.
  • the process (filter process) of (5) [or (7)] N times or more it is preferable to perform the process (filter process) of (5) [or (7)] N times or more.
  • the order of the integration operation and the filtering process may be any order as long as the following conditions are satisfied, for example. ⁇ Do not perform the filter process more than once in succession. ⁇ At least once in the filter process, the integration process is completed. -When the filter process is performed twice or more, the filter is a high-pass filter except for the last one.
  • the evaluation coefficient determination unit 213 determines an evaluation coefficient for evaluating the synchronization of the motor rhythm based on, for example, the phase of the motor rhythm.
  • the phase is an index representing how many rotation angles each point on the rhythm waveform corresponds to, assuming that a rhythm that repeats in time is a movement rotating on the circumference. For example, there is a difference of 360 degrees between the phases of two adjacent peak points in a sine wave.
  • FIG. 2 is a diagram illustrating an example of an exercise rhythm, and the peak in the figure corresponds to, for example, the landing of the subject's right foot or the left foot.
  • an example of the evaluation coefficient is, for example, the peak corresponding to one foot and the peak corresponding to the other foot between the peaks corresponding to one foot (for example, the right foot) in the example of the exercise rhythm shown in FIG.
  • the evaluation coefficient for evaluating the synchrony of the exercise rhythm is not limited to this evaluation coefficient.
  • the value of the evaluation coefficient when the first time T1 and the second time T2 are equal is 0, but is not limited to 0.
  • an evaluation coefficient calculation formula may be used so that the value of the evaluation coefficient when the first time T1 and the second time T2 are equal is a value other than 0 such as 1 or the like. That is, the calculation formula for the evaluation coefficient is not limited to the following formula (1).
  • the following formula (1) indicates that the difference between the time from the landing of one foot to the landing of the other foot and the time from the landing of the other foot to the landing of one foot, and the landing of one foot Represents the ratio of the time from the landing of the other foot to the sum of the time from the landing of the other foot to the landing of the one foot.
  • the time intervals between the peak corresponding to one foot and the peak corresponding to the other foot between the peaks corresponding to one foot are respectively the first time T1 and the second time.
  • T2 it is not limited to the time interval between peaks.
  • the first time T1 may be a time until the phase changes by 360 degrees from a predetermined position that does not correspond to the peak of the movement rhythm.
  • the second time T2 may be a time from a position where the phase has changed 360 degrees from a predetermined position not corresponding to the peak of the movement rhythm to a further change in phase by 360 degrees.
  • the first time T1 is the time until the phase changes 360 degrees from the peak of the movement rhythm or a predetermined position not corresponding to the peak, and the second time T2 does not correspond to the peak or peak of the movement rhythm. Although it is the time from the position where the phase has changed 360 degrees from the predetermined position until the phase has changed 360 degrees, the present invention is not limited to this.
  • the first time T1 is a time until the phase changes by -360 degrees from a predetermined position that does not correspond to the peak of the movement rhythm or a peak (a time that goes back 360 degrees from the predetermined position)
  • the second time T2 Is the time from when the phase changes by -360 degrees from a predetermined position that does not correspond to the peak of the motor rhythm or when the phase changes by -360 degrees (time that goes back from the predetermined position by 720 degrees). May be.
  • the first time T1 is a time until the phase changes by 360 degrees from a predetermined position not corresponding to the peak of the movement rhythm or the peak
  • the second time T2 is a peak or peak of the movement rhythm. It may be a time until the phase changes by -360 degrees from a predetermined position that does not correspond (a time that is 360 degrees backward from the predetermined position).
  • the first time T1 when the first time T1 is set as a time that goes back 360 degrees from the predetermined position, or when the second time T2 is set as a time that goes back 360 degrees or 720 degrees from the predetermined position, the first time T1 and The value of the second time T2 may be negative. In such a case, it is necessary to pay attention to the first time T1 and the second time T2 that are substituted into Equation 1.
  • the evaluation coefficient determination unit 213 includes a time calculation unit 223 and an evaluation coefficient calculation unit 233.
  • the time calculation unit 223 calculates the movement rhythm phase, thereby changing the phase from the movement rhythm phase at a predetermined time (arbitrary time) by a predetermined angle (hereinafter sometimes referred to as a first angle). (First time) and a time (second time) at which the phase changes twice the first angle from the phase of the motion rhythm at a predetermined time.
  • the time calculation unit 223 can calculate, for example, a time from a predetermined time to the first time (first time: T1). In addition, the time calculation unit 223 can calculate, for example, a time from the first time to the second time (second time: T2).
  • the first angle is, for example, a positive value of 180 degrees or 360 degrees.
  • 360 degrees includes a case of strictly 360 degrees and a case of approximately 360 degrees.
  • 180 degrees includes a case of strictly 180 degrees and a case of approximately 180 degrees. Including.
  • the time calculation unit 223 calculates, for example, the phase of the movement rhythm, for example, the time that the phase goes back the first angle from the phase of the movement rhythm at a predetermined time (arbitrary time) And the time when the phase has changed twice the first angle from the phase of the movement rhythm at the predetermined time (the time that is double the first angle is also included in the second time) May be calculated. That is, the first angle may be a negative value of ⁇ 180 degrees or ⁇ 360 degrees, for example.
  • the time calculation unit 223 can calculate the time from the first time to the predetermined time (first time: T1).
  • the time calculation unit 223 can calculate the time from the second time to a predetermined time (second time: T2).
  • the time calculation unit 223 has a phase (first time) when the phase changes from the phase of the motion rhythm at a predetermined time (arbitrary time) and a phase from the phase of the motion rhythm at the predetermined time. You may make it calculate the time which changed -1 time of the angle. In this case, the time calculation unit 223 can calculate the time from the predetermined time to the first time (first time: T1). In addition, the time calculation unit 223 can calculate the time from the second time to a predetermined time (second time: T2).
  • the evaluation coefficient calculation unit 233 calculates the first time (the first time is larger than the predetermined time) and the second time (the second time is larger than the first time) will be described. .
  • the calculation method is the same in the case of calculating the (large) and second time (the second time is smaller than the predetermined time).
  • the time calculation unit 223 calculates the phase of the motion rhythm using, for example, the Hilbert transform method or the pattern matching method.
  • the method for calculating the phase is not limited to the Hilbert transform method or the pattern matching method.
  • the Hilbert transform method is a method for specifically calculating the phase at a predetermined position of the waveform
  • the pattern matching method is not a method for directly specifying the phase, but is a method for finding points where the phases are mutually shifted by 360 degrees. is there.
  • the pattern matching method is a method of finding a point whose phase is shifted by an integral multiple of 360 degrees from a predetermined position of the waveform.
  • the pattern matching method calculates a relative phase (for example, 360 degrees) with respect to a predetermined position of the waveform.
  • the evaluation coefficient calculation unit 233 calculates the evaluation coefficient based on the first time and the second time calculated by the time calculation unit 223, for example. For example, the evaluation coefficient calculation unit 233 calculates the first time and the second time based on the first time and the second time calculated by the time calculation unit 223, and sets the first time and the second time to the above (1). An evaluation coefficient is calculated by substituting it into the equation. That is, the evaluation coefficient calculation unit 233 calculates an evaluation coefficient for evaluating the synchrony of repeated voluntary movement based on the first time and the second time.
  • the Hilbert transform is a mathematical method for deriving a corresponding imaginary part Y (t) from an arbitrary real time series signal X (t).
  • the time calculation unit 223 can directly obtain the phase ⁇ (t) of X (t) from the following equation (2). That is, the phase calculation unit 224 can obtain the phase of the motion rhythm at a predetermined time (predetermined position) from the following equation (2).
  • evaluation coefficient calculation part 233 calculates an evaluation coefficient from the following (3) formula, for example using the time t1 and the time t2 which were calculated by the time calculation part 223, for example.
  • the pattern matching method is a method for quantifying the similarity between two signals. There are many methods for defining or calculating the degree of similarity. Specifically, for example, a method described in “Image processing engineering (Ryoichi Suematsu, Hirohisa Yamada, Corona)” or the like is used. The most representative is the autocorrelation coefficient. For example, the following calculation is performed for a three-dimensional body motion signal.
  • the autocorrelation coefficient is calculated by the following equation (5).
  • the same calculation is performed for a one-dimensional signal.
  • the pattern matching method it is possible to easily obtain the time t1 when the phase of the rhythm waveform at an arbitrary time t is shifted by 360 degrees, the time t2 when the phase is shifted by 720 degrees, or the like.
  • the time calculation unit 223 can obtain a time whose phase is shifted by an integer multiple of 360 degrees from the phase of the motion rhythm at the time t.
  • the evaluation coefficient calculation unit 233 calculates the evaluation coefficient from the above equation (3), for example, using the time t1 and the time t2 calculated by the time calculation unit 223, for example. Note that the evaluation coefficient calculation unit 233 may calculate the evaluation coefficient from the equation (1) using the first time T1 and the second time T2 calculated based on the time t1 and the time t2.
  • Determining unit 214 evaluates (determines) exercise rhythm synchrony based on the evaluation coefficient obtained by evaluation coefficient determining unit 213, for example. That is, the determination unit 214 evaluates the synchronization or balance of the exercise rhythm based on the evaluation coefficient. Specifically, for example, when the voluntary movement is walking or the like, the determination unit 214 evaluates the synchronization of left and right footsteps, that is, the synchronization or balance of left and right footsteps. Here, whether the left and right footing is synchronized or not is determined based on whether the evaluation coefficient maintains a constant value regardless of time or whether the evaluation coefficient changes periodically.
  • the determination unit 214 determines that the left and right footsteps are synchronized when the evaluation coefficient maintains a constant value regardless of time or when the evaluation coefficient changes periodically. Further, whether the right / left footing balance is good or bad is determined based on whether or not the evaluation coefficient is close to 0 when the evaluation coefficient is determined using the equation (1). For example, the determination unit 214 determines that the left and right footing balance is better as the evaluation coefficient is closer to zero.
  • the determination unit 214 calculates, for example, the following three indices R, S, and C from the change in the evaluation coefficient over a predetermined time.
  • the predetermined time is a time required for the subject to walk about 10 steps (for example, 5 seconds), but is not limited to this time.
  • the determination unit 214 evaluates synchrony based on these three indices R, S, and C.
  • FIG. 3 is a diagram for explaining an example of the operation of the determination unit 214.
  • the determination unit 214 has a very good balance between left and right footing (walking rhythm) and right and left footing (in FIG. 3). Judgment). That is, the determination unit 214 functions as a determination unit that determines whether the evaluation coefficient is within a predetermined range (for example, 0 ⁇ 0.02) over a predetermined time (for example, 5 seconds). Further, for example, when R is not less than 0.02 but S is not more than 0.01, the determination unit 214 has a good balance between right and left footsteps and a right and left footstep balance (FIG. 3). Judgment)
  • the determination unit 214 synchronizes left and right footsteps and right and left footsteps based on C. Assess balance. For example, when C is 0.5 or more, the determination unit 214 determines that the left and right footsteps are synchronized, but the left and right footsteps are not well balanced (see x in FIG. 3), and C is 0. If it is smaller than .5 and greater than 0.2, it is determined that the left and right footsteps are normally synchronized (see ⁇ in FIG. 3), but the left and right footsteps are not well balanced (see ⁇ in FIG. 3). Furthermore, for example, when C is less than 0.2, the determination unit 214 determines that the left and right footsteps are not synchronized (see “X” in FIG. 3), and the right and left footstep balance cannot be determined. .
  • the numerical values (0.02, 0.01, 0.5, and 0.2) in FIG. 3 are examples, and are not limited to these numerical values. Further, in FIG. 3, five patterns are determined using three indexes R, S, and C, but the determination method is not limited to this, and finer determination may be performed. A rough evaluation may be performed. For example, two patterns may be determined based only on the value of R.
  • the determination unit 214 can directly index the synchronization of left and right footsteps and the balance of left and right footsteps from the change in the evaluation coefficient over a predetermined time, for example, as follows. For example, when the voluntary movement is walking, the predetermined time is a time required for the subject to walk about 10 steps (for example, 5 seconds), but is not limited to this time.
  • H (t) the evaluation coefficient calculated from the phase of the movement rhythm is denoted as H (t) as a function of time t.
  • t2 be the time when the phase of the motion rhythm has changed by 720 degrees starting from an arbitrary time t1.
  • S1 The standard deviation of H (t) between time t1 ⁇ t ⁇ time t2 is obtained.
  • the time change of the standard deviation is calculated by sequentially changing the time t1 over a predetermined time. With respect to the time series data of the standard deviation thus obtained, an average value within a predetermined time is obtained and used as an index S1 for the balance of left and right foot travel.
  • the phase obtained by the time calculation unit 223 or the index obtained by the determination unit 214 can be used as the quantification of the regularity of the output waveform Y in the time constant determination unit 211 described above.
  • the CV of the phase period or the synchronization index S2 of the left and right footing is considered to reflect the regularity of the output waveform Y.
  • the time constant A of the filter is changed (for example, between 0 and 1 second in the case of walking rhythm)
  • the CV or index S2 of the phase period is obtained and the value takes the minimum value.
  • the constant A is determined as the optimal time constant.
  • the output control unit 215 controls the output unit 23.
  • the output control unit 215 controls the display state of the output unit 23 to display various information on the output unit 23.
  • the output control unit 215 causes the output unit 23 to display the evaluation coefficient obtained by the exercise rhythm extracted by the rhythm extraction unit 212 or the evaluation coefficient determination unit 213.
  • the evaluation coefficient may be displayed as a graph on the output unit 23 or may be displayed as a message such as “Evaluation coefficient is **”.
  • the output control unit 215 may cause the output unit 23 to display the determination result by the determination unit 214 in addition to the display of the evaluation coefficient or in place of the display of the evaluation coefficient.
  • the output control unit 23 warns, for example, by an alarm or vibration, the output control unit when the determination unit 214 determines that the evaluation coefficient is not a value within a predetermined range, for example. 215 controls the output unit 23 to cause the output unit 23 to issue a warning.
  • the output control unit 215 is, for example, an evaluation coefficient obtained by the exercise rhythm or evaluation coefficient determination unit 213 extracted by the rhythm extraction unit 212.
  • the display state of the output unit 23 is controlled by transmitting various information such as the above to the output unit 23 via wireless or wired communication.
  • the output unit 23 outputs various types of information under the control of the output control unit 215.
  • the output unit 23 is a display, and when the evaluation coefficient obtained by the evaluation coefficient determination unit 213 or the determination unit 214 determines that the evaluation coefficient is not a value within a predetermined range, for example. , Display a warning.
  • the output unit 23 may be a unit that notifies about a change in state or a sudden abnormality by, for example, an alarm or vibration, and the evaluation unit 214 does not have a value within a predetermined range, for example. If it is determined, a warning is given by an alarm or vibration.
  • the output unit 23 may not be provided in the information processing apparatus 20 but may be provided outside the information processing apparatus 20.
  • the case where the output unit 23 is provided outside the information processing device 20 is, for example, the case where the output unit 23 is included in the body motion signal detection device 10, or the body motion signal detection device 10 and the information processing. This is a case where the device 20 is not included in either.
  • the information processing apparatus 20 is connected to the output unit 23 via wireless or wired connection. That is, the output unit 23 functions as an output unit that outputs a result (for example, an evaluation coefficient) obtained by the information processing apparatus.
  • the output unit 23 functions as an output device including an output unit that outputs a result (for example, an evaluation coefficient) obtained by the information processing device.
  • the interface unit 24 is, for example, an interface that is communicably connected to the body motion signal detection device 10.
  • the interface unit 24 includes an antenna.
  • the interface unit 24 is a wired connection terminal.
  • the storage unit 12 is provided detachably with respect to the body motion signal detection device 10, and when the storage unit 12 is detached from the body motion signal detection device 10 and connected to the information processing device 20, The interface unit 24 also functions as a slot to which the storage unit 12 is connected, for example.
  • steps A1 to A6 The operation of the system 1 as an example of the embodiment configured as described above will be described with reference to the flowchart (steps A1 to A6) shown in FIG.
  • steps A1 to A6 an example of a specific method for evaluating the synchronization of movement rhythm based on the integral signal will be described in detail for the case of using an acceleration signal during walking. The same processing can be applied to the above signal.
  • the body motion signal detection unit 11 detects a repetitive rhythm motion accompanying a voluntary motion as a body motion signal (step A1).
  • FIG. 5A shows that the body motion signal detection device 10 (body motion signal detection unit 11) performs 100 Hz sampling with the body motion signal detection device 10 (triaxial acceleration sensor) attached to the center of the abdomen of a healthy subject. It is a part of the acceleration signal during walking for 15 minutes measured by. In FIG. 5A, only the acceleration change in the vertical direction among the three axes is shown.
  • the time constant determination unit 211 determines a time constant for extracting an exercise rhythm (step A2).
  • the rhythm extraction unit 212 extracts an exercise rhythm from the body motion signal (step A3).
  • FIG. 5B is a diagram illustrating an example of an exercise rhythm extracted from a body motion signal.
  • the evaluation coefficient determination unit 213 determines an evaluation coefficient from the extracted exercise rhythm (step A4). Based on this evaluation coefficient, the determination unit 214 determines, for example, the synchronization of the exercise rhythm (step A5: determination process). Next, based on the determination result of the determination unit 214, the output control unit 215 causes the output unit 23 to output the determination result (step A6). Next, details of the time constant determination unit 211, that is, detailed operations of step A2 in FIG. 2 will be described with reference to the flowchart (steps A21 to A24) shown in FIG.
  • the time constant determination unit 211 performs filtering on the body motion signal detected by the body motion signal detection unit 11 using a predetermined time constant A (step A21).
  • the time constant determination unit 211 obtains the absolute value CV of the values at the maximum and minimum peaks of the output waveform Y, for example (step A22). Then, the time constant determining unit 211 obtains a change in CV with respect to a change in the time constant A by changing the time constant A and performing the processes of steps A21 and A22, for example (step A23).
  • FIG. 7 is a diagram illustrating an example of a change in CV (dispersion) with respect to a change in time constant (filter time) A obtained in step A23.
  • FIG. 7 is a diagram illustrating an example of a change in CV (dispersion) with respect to a change in time constant (filter time) A when the body motion signal includes a respiratory rhythm.
  • the time constant determination unit 211 obtains a minimum point from the change in CV with respect to the change in the time constant T obtained in step A23. For example, as shown in FIG. 7, the time constant determination unit 211 obtains two local minimum points (marks ⁇ and ⁇ in FIG. 7) having a small CV value. Then, the time constant determination unit 211 calculates, for example, a time constant value (for example, 0.4) in the vicinity of the time constant A corresponding to the minimum point (marked with ⁇ in FIG. 7) having the smaller time constant A. It is determined as a time constant (step A24). As described above, the time constant determination unit 211 may determine the time constant in the vicinity of the minimum point as the first time constant instead of setting the time constant corresponding to the minimum point itself as the first time constant. Even in this case, there is no great difference in the results obtained.
  • a time constant value for example, 0.4
  • step A37 the rhythm extraction unit 212 performs high-pass filtering on the body motion signal detected by the body motion signal detection unit 11 using the first time constant (step A37).
  • the rhythm extraction unit 212 performs N (for example, 2) order integration (step A38).
  • step A39 high-pass filtering is again performed on the signal after integration using the first time constant (step A39).
  • FIG. 5B is a diagram illustrating an example of the signal obtained in step A39. This waveform corresponds to the voluntary movement of the subject who is walking, for example. More specifically, the waveform obtained in step 39 corresponds to the relative motion trajectory of each step by walking. That is, the waveform shown in FIG. 5B corresponds to a waveform indicating a movement rhythm.
  • the time calculation unit 223 selects an arbitrary point at an arbitrary time t of the exercise rhythm (step A41).
  • the time calculation unit 223 calculates the phase using, for example, Hilbert transform or pattern matching, so that the time t1 at which the phase is shifted 360 degrees from the arbitrary point and the phase is shifted 720 degrees from the arbitrary point
  • the point time t2 is calculated (step A42: time calculation process).
  • the evaluation coefficient calculation unit 233 calculates the evaluation coefficient from, for example, the above equation (1) using the times t1 and t2 (step A43: evaluation coefficient calculation process). That is, in step A42, from the signal representing the repetitive voluntary movement of the living body, a time (first time) when the phase of the signal changes by a predetermined angle (first angle) from the phase of the signal at a predetermined time, It is an example of the time calculation process which calculates the time (2nd time) when the phase of the signal changed twice the 1st angle from the phase of the signal in time.
  • Step A43 is an example of a second calculation process for calculating an evaluation coefficient for evaluating the synchrony of repeated voluntary movement based on the first time and the second time.
  • FIGS. 10A to 10C are diagrams for explaining an example of a signal obtained when the pattern matching method is used.
  • An example of detailed operation of the evaluation coefficient determination unit 213 when pattern matching is used will be described with reference to FIGS. 10 (A) to 10 (C).
  • FIG. 10A is a partial waveform extracted from the data of FIG.
  • FIG. 10B shows the result of the time calculation unit 223 selecting the reference wave having a width of 0.4 seconds around the time indicated by * and calculating the autocorrelation coefficient.
  • FIG. 10C shows an autocorrelation coefficient obtained from a reference wave having a width of 0.4 seconds around the circle mark in FIG.
  • the positions of three points whose phases are shifted by 360 degrees are indicated by dotted lines.
  • the evaluation coefficient calculation unit 233 calculates the first time T1 and the second time T2 from the waveforms shown in FIG. 10B or FIG. 10C, and sets the times T1 and T2 to the above formula (1).
  • the evaluation coefficient is calculated by substituting for.
  • the result of an example in which the evaluation coefficient is obtained every time the phase of the motion rhythm increases by 3.6 degrees using the Hilbert transform method is shown by the solid line in FIG. This is equivalent to obtaining an evaluation coefficient every 0.05 seconds on average when converted to a time interval. From FIG. 11, it can be seen that the evaluation coefficient fluctuates periodically around 0.
  • the evaluation coefficient determination unit 213 calculates an evaluation coefficient over a predetermined time (step A51), that is, the above steps A41 to A44 are repeated over a predetermined time.
  • the determination unit 214 calculates three indexes R, S, and C using the evaluation coefficient calculated over a predetermined time (step A52). Then, the determination unit 214 determines the synchronization of the exercise rhythm based on these three indicators (step A53).
  • step A531 determines whether or not R is 0.02 or less (step A531).
  • step A532 determines that the synchronization and balance of the exercise rhythm are very good, assuming that the evaluation coefficient is substantially zero over a predetermined time (step S531).
  • step A531 is an example of a determination process for determining whether the evaluation coefficient is within a predetermined range (for example, 0 ⁇ 0.02) over a predetermined time.
  • the determination unit 214 determines whether S is 0.01 or less (Step A533). When S is 0.01 or less (see the Yes route in step A533), the determination unit 214 determines that the synchronization and balance of the exercise rhythm are good, assuming that the evaluation coefficient is not substantially 0 but remains constant. (Step A534). On the other hand, when S is larger than 0.01 (see No route in step A533), the determination unit 214 determines whether C is 0.5 or more (step A535).
  • step A535 When C is equal to or greater than 0.5 (see the Yes route in step A535), it is determined that the evaluation coefficient fluctuates periodically, and the determination unit 214 determines that the exercise rhythm is synchronized but the balance is poor (step A536). ). On the other hand, when C is less than 0.5 (see No route in step A535), the determination unit 214 determines whether C is 0.2 or more (step A537). When C is 0.2 or more (see the Yes route in step A537), the determination unit 214 determines that the synchronization of the exercise rhythm is normal and the balance is bad (step A538). On the other hand, when C is less than 0.2 (see No route of step A537), the determination unit 214 determines that the synchronization of the exercise rhythm is poor and the balance cannot be determined (step A539).
  • the problem is the effect of data boundaries. That is, when signal processing such as spectrum analysis or Hilbert transform is performed, the error increases near the end points (measurement start point and end point) of the data. That is, the evaluation coefficient at the current time may include an error based on signal processing.
  • the result of calculating the evaluation coefficient from the phase change by the Hilbert transform method using only the data up to this time as the current time for the acceleration data of FIG. This is indicated by a broken line.
  • solid line in FIG. 14 Compared with the result of the accurate calculation method using the entire 15-minute data (solid line in FIG. 14), a value with almost no error is obtained in a time that goes back about 1.3 seconds from the current time. This is sufficient performance for practical real-time processing.
  • the evaluation coefficient for evaluating the synchrony of the movement rhythm is calculated based on the change of the phase by paying attention to the phase, a clear peak An evaluation coefficient can be calculated from data that does not have rhythm, and the synchronization of the motor rhythm can be evaluated. That is, according to the system 1 in the example of the present embodiment, since attention is paid to the phase, the evaluation coefficient can be obtained and the synchronization of the exercise rhythm can be evaluated without paying attention to the peak of the exercise rhythm.
  • the evaluation coefficient can be obtained almost continuously at an arbitrary time, and further, the evaluation coefficient can be obtained in real time. Can be requested. That is, according to the system 1 in the example of the present embodiment, the synchronization of the movement rhythm can be continuously evaluated, and further, the synchronization of the movement rhythm can be evaluated in real time.
  • the waveform after autocorrelation has less noise than the original body motion signal and thus has a clear peak. Since the peak position can be accurately specified, the evaluation coefficient can be obtained with high accuracy.
  • the signal is not limited to the signal from only the acceleration change in the vertical direction.
  • an evaluation coefficient may be calculated by extracting a motion rhythm from an acceleration signal in the front-rear direction or the left-right direction.
  • a point where the phase is shifted by 360 degrees from the predetermined position in the range where the phase changes from the predetermined position by 720 degrees is searched, and the evaluation coefficient is calculated.
  • it is not limited to 360 degrees.
  • the adjacent acceleration signal in the left and right direction corresponds to the landing of the same foot as shown in FIG.
  • the point where the phase is shifted from the predetermined position by 180 degrees is searched, and the evaluation coefficient is calculated with the time interval between them as the first time T1 and the second time T2, respectively.
  • the same processing as when calculating the evaluation coefficient from the lateral acceleration signal is performed.
  • the body motion signal detection device 10 does not include the central processing unit 21, but is not limited to this configuration, and the body motion signal detection device 10 includes the central processing unit 21. May be.
  • the central processing unit 21 executes an information processing program stored in a storage device (for example, the storage unit 12) in the body motion signal detection device 10 or a storage unit (not shown) outside the body motion signal detection device 10.
  • a storage device for example, the storage unit 12
  • a storage unit not shown
  • the body motion signal detection device 10 includes the body motion signal detection unit 11 and the storage unit 12, but is not limited to this configuration.
  • the body motion signal detection device 10 may include the body motion signal detection unit 11 but may not include the storage unit 12.
  • the body motion signal detection device 10 and the storage unit 12 are connected by wire or wirelessly, and the body motion signal detected by the body motion signal detection unit 11 is transferred to the storage unit 12 via wire or wirelessly.
  • the storage unit 12 and the information processing device 20 are connected by wire or wirelessly, and the information processing device 20 acquires a body motion signal from the storage unit 12.
  • the body motion signal detection device 10 detects the body detected by the body motion signal detection unit 11. It has a function of transmitting a moving signal to the information processing apparatus 20 (or the storage unit 12) via the interface unit 13 that is an antenna. This transmission function is realized by executing a program stored in a storage unit (not shown) by a processing unit (not shown) provided by the body motion signal detection device 10 such as a central processing unit.
  • the information processing apparatus 20 determines the first time constant from the body motion signal and extracts the exercise rhythm based on the first time constant, but is limited to this extraction method. It is not a thing.
  • the frequency characteristic of the body motion signal is obtained by spectrum analysis such as FFT (Fast Fourier Transform) or wavelet analysis, and the frequency range of the motion rhythm is specified from the result of the spectrum analysis. Then, for example, a movement rhythm may be extracted by applying a bandpass filter corresponding to the specified frequency range to the body motion signal.
  • the body motion signal is decomposed into each mode waveform (Intrinsic Mode Function) by EMD (Empirical Mode Decomposition) or EEMD (Ensemble Empirical Mode Decomposition), and from this result, the mode waveform corresponding to the exercise rhythm (for example, strong intensity) (Waveform) may be selected.
  • the rhythm extraction unit 212 extracts the body motion rhythm from the body motion signal by the filtering process, but this is not limitative. It is not a thing. For example, in the process (3), all output waveforms Y when the time constant T is changed are stored in the storage unit 22. The rhythm extraction unit 212 may select a waveform corresponding to the first time constant determined in the process (4) from all the output waveforms Y.
  • An example of this embodiment is applicable not only to people but also to animals such as pets, livestock, and horses.
  • an example of the present embodiment mainly describes the case where the voluntary movement is walking, the present invention is not limited to walking, and the example of the present embodiment can be applied to other voluntary movements. For example, if the voluntary exercise is juggling, it is only necessary to focus on the exercise rhythm of the hand, not the foot.
  • Various application programs for realizing each function of the central processing unit 21 provided in the information processing apparatus 20 are, for example, a flexible disk, a CD (CD-ROM, CD-R, CD-RW, etc.), DVD (DVD -Recorded in a computer-readable recording medium such as a ROM, DVD-RAM, DVD-R, DVD + R, DVD-RW, DVD + RW, HD DVD), Blu-ray disc, magnetic disc, optical disc, magneto-optical disc, etc.
  • the computer reads the program from the recording medium, transfers it to the internal storage device or the external storage device, and uses it.
  • the program may be recorded in a storage device (recording medium) such as a magnetic disk, an optical disk, or a magneto-optical disk, and provided from the storage device to the computer via a communication path.
  • Various application programs for realizing each function of a central processing unit (not shown) provided in the body motion signal detection apparatus 10 are, for example, a flexible disk, a CD (CD-ROM, CD-R, CD-RW, etc.), Recorded on computer-readable recording media such as DVD (DVD-ROM, DVD-RAM, DVD-R, DVD + R, DVD-RW, DVD + RW, HD DVD, etc.), Blu-ray disc, magnetic disc, optical disc, magneto-optical disc, etc. Provided in different forms. Then, the computer reads the program from the recording medium, transfers it to the internal storage device or the external storage device, and uses it. The program may be recorded in a storage device (recording medium) such as a magnetic disk, an optical disk, or a magneto-optical disk, and provided from the storage device to the computer via a communication path.
  • a storage device recording medium
  • FIG. 16A to FIG. 16C are diagrams showing walking results when on.
  • FIG. 16A is a diagram illustrating acceleration detected by the body motion signal detection device 10.
  • FIG. 16B is a diagram showing an acceleration signal after second-order integration, that is, a motion trajectory.
  • 16C is a diagram showing the evaluation coefficient obtained from the acceleration signal after second-order integration.
  • the evaluation coefficient since the evaluation coefficient has a regular fluctuation consistently, it can be said that there is a left-right difference in the walking step, but the left-right footing synchronization is high.
  • FIGS. 17 (A) to 17 (C) are diagrams showing walking results when the walking is relatively stable during the time of 3 to 17 seconds, but sudden progress is recognized before and after the walking. It is.
  • FIG. 17A is a diagram illustrating acceleration detected by the body motion signal detection device 10.
  • FIG. 17B is a diagram illustrating an acceleration signal after the second-order integration, that is, a motion trajectory.
  • FIG. 17C is a diagram showing the evaluation coefficient obtained from the acceleration signal after second-order integration.
  • the solid line indicates the evaluation coefficient calculated using the Hilbert transform method
  • the broken line indicates the evaluation coefficient calculated using the pattern matching method. Comparing FIG. 16 (A) to FIG. 16 (C) and FIG. 17 (A) to FIG. 17 (C), in FIG.
  • a small wireless hybrid sensor WAA-006 manufactured by Wireless Technology as a body motion signal detection device 10 in the center of the abdomen of a healthy subject, sampling a 3-axis acceleration signal and 3-axis angular velocity signal when walking in the city with a heavy bag on the shoulder Simultaneous measurement at a frequency of 200 Hz.
  • the acceleration signal in the vertical direction and the angular velocity signal around the vertical axis are: • High-pass filtering using time constant 1 • Second-order integration • Hilbert transform method after extracting motion rhythm by applying high-pass filtering using time constant 1 To obtain the phase. Evaluation coefficients were calculated from two points where the phase difference was 360 degrees for the acceleration signal and two points where the phase difference was 180 degrees for the angular velocity signal.
  • FIG. 18A is a diagram showing acceleration and angular velocity detected by the body motion signal detection device 10.
  • FIG. 18B is a diagram showing an acceleration signal and an angular velocity signal after second-order integration.
  • FIG. 18C is a diagram showing evaluation coefficients obtained from the acceleration signal and the angular velocity after the second-order integration.
  • the solid line indicates the acceleration signal
  • the broken line indicates the angular velocity signal.
  • the evaluation coefficient shows almost the same fluctuation in strength. It can also be seen that the balance between the left and right is getting worse, reflecting the heavy bag on his shoulder.
  • the synchrony of body movement rhythm can be evaluated regardless of the type of inertial sensor.
  • an acceleration recorder “Watching Gate” manufactured by Mitsubishi Chemical Corporation was put on a dedicated belt and wound around the abdomen of a Parkinson's disease patient, and body motion signals were continuously sampled at a sampling frequency of 100 Hz for 38 hours. At the same time, the subjects were asked to fill in a diary of the degree of ease of movement and the time when they fell.
  • FIG. 19A shows the time change of the walking index after the smoothing process is performed.
  • the subject's self-report is an intermediate evaluation of “not easy to move” or “easy to move” — “not easy to move” at times ta and tb. A fall is happening.
  • the walking index is also low. From 17:00 to 18:00 on the first day, from 10:00 to 11:00 on the second day, the subject's self-report is “easy to move”, but the gait index is lower than the time zone before and after, A fall occurred at times tc and td. Moreover, it can be seen that the walking index tends to decrease considerably before the fall time.
  • the synchrony of walking rhythm is an index for predicting the risk of falls.
  • FIG. 19B shows a time change of the body motion index after the smoothing process is performed.
  • the subject's self-reported five-step evaluation of ease of movement (movable, difficult to move, unable to move in three steps, easy to move-intermediate evaluation of difficult to move, and difficult to move-unmovable 5 grades that include 2 grades of the intermediate assessment)), and the actual fall times are indicated as ta, tb, tc, and td.
  • the transition of the body motion index is indicated by a dotted line
  • the transition of the subject's self-report is indicated by a solid line.
  • FIG. 19 (B) in the vicinity of the fallen area, the same behavior as in FIG. 19 (A) is shown. That is, between 9 am and 10 am on the first day, the subject's self-report is an intermediate evaluation of “not easy to move” or “easy to move” — “not easy to move”, and a fall occurred at times ta and tb. In accordance with this, the body motion index is also low. From 17:00 to 18:00 on the first day, from 10:00 to 11:00 on the second day, the subject's self-report is “easy to move”, but the body motion index is lower than before and after, A fall occurred at times tc and td. Moreover, it can be seen that the body motion index tends to decrease considerably before the fall time.
  • the synchrony of body movement rhythm is an index for predicting the risk of falls.
  • the synchronization of the movement rhythm is continuously evaluated in real time. There is an effect that can be done.

Landscapes

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

Abstract

Le synchronisme des rythmes des mouvements du corps est continuellement évalué. Le procédé comprend : un premier procédé de calcul pour calculer le temps auquel une phase a changé d'un angle prédéterminé par rapport à la phase d'un signal à un temps prédéterminé, et calculer le temps auquel la phase a changé de deux fois l'angle prédéterminé par rapport à la phase du signal au temps prédéterminé ; et un second procédé de calcul pour calculer un coefficient d'évaluation sur la base de ces temps.
PCT/JP2011/070763 2010-09-17 2011-09-12 Procédé de traitement de l'information, dispositif de traitement de l'information, dispositif de sortie, système de traitement de l'information, programme de traitement de l'information et support d'enregistrement lisible par ordinateur sur lequel le même programme est enregistré Ceased WO2012036135A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012534001A JPWO2012036135A1 (ja) 2010-09-17 2011-09-12 情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010209112 2010-09-17
JP2010-209112 2010-09-17

Publications (1)

Publication Number Publication Date
WO2012036135A1 true WO2012036135A1 (fr) 2012-03-22

Family

ID=45831594

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2011/070763 Ceased WO2012036135A1 (fr) 2010-09-17 2011-09-12 Procédé de traitement de l'information, dispositif de traitement de l'information, dispositif de sortie, système de traitement de l'information, programme de traitement de l'information et support d'enregistrement lisible par ordinateur sur lequel le même programme est enregistré

Country Status (2)

Country Link
JP (1) JPWO2012036135A1 (fr)
WO (1) WO2012036135A1 (fr)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014193349A (ja) * 2013-03-19 2014-10-09 Avita Corp 生理状態を監視するための装置および方法
WO2014181604A1 (fr) * 2013-05-10 2014-11-13 オムロンヘルスケア株式会社 Dispositif de mesure de posture de marche et programme associé
JP2016059729A (ja) * 2014-09-22 2016-04-25 カシオ計算機株式会社 測定装置、測定方法及び測定プログラム
JP2016112053A (ja) * 2014-12-11 2016-06-23 国立研究開発法人産業技術総合研究所 歩行状態判定方法、プログラム、及び装置
US11607138B2 (en) * 2012-06-01 2023-03-21 Vital Connect, Inc. Respiratory rate detection using decomposition of ECG
JP2023549625A (ja) * 2020-08-28 2023-11-29 ローベルト ボツシユ ゲゼルシヤフト ミツト ベシユレンクテル ハフツング 水泳ストロークを決定するためのコントローラ及び方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007244495A (ja) * 2006-03-14 2007-09-27 Sony Corp 体動検出装置、体動検出方法および体動検出プログラム
WO2008035611A1 (fr) * 2006-09-19 2008-03-27 Mitsubishi Chemical Corporation Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données
WO2010058535A1 (fr) * 2008-11-18 2010-05-27 オムロンヘルスケア株式会社 Dispositif, programme et procédé de détection d’équilibre lors de mouvements corporels, et procédé pour diagnostiquer l’équilibre lors de mouvements corporels

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007244495A (ja) * 2006-03-14 2007-09-27 Sony Corp 体動検出装置、体動検出方法および体動検出プログラム
WO2008035611A1 (fr) * 2006-09-19 2008-03-27 Mitsubishi Chemical Corporation Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données
WO2010058535A1 (fr) * 2008-11-18 2010-05-27 オムロンヘルスケア株式会社 Dispositif, programme et procédé de détection d’équilibre lors de mouvements corporels, et procédé pour diagnostiquer l’équilibre lors de mouvements corporels

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
HIDEYA TAKANASHI: "Co-emergence Robot Walk- Mate and Its Support for Elderly People", TRANSACTIONS OF THE SOCIETY OF INSTRUMENT AND CONTROL ENGINEERS, vol. 39, no. 1, January 2003 (2003-01-01), pages 74 - 81 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11607138B2 (en) * 2012-06-01 2023-03-21 Vital Connect, Inc. Respiratory rate detection using decomposition of ECG
JP2014193349A (ja) * 2013-03-19 2014-10-09 Avita Corp 生理状態を監視するための装置および方法
WO2014181604A1 (fr) * 2013-05-10 2014-11-13 オムロンヘルスケア株式会社 Dispositif de mesure de posture de marche et programme associé
JP2014217694A (ja) * 2013-05-10 2014-11-20 オムロンヘルスケア株式会社 歩行姿勢計およびプログラム
CN105208932A (zh) * 2013-05-10 2015-12-30 欧姆龙健康医疗事业株式会社 步态仪和程序
US10898112B2 (en) 2013-05-10 2021-01-26 Omron Healthcare Co., Ltd. Gait posture meter and program
JP2016059729A (ja) * 2014-09-22 2016-04-25 カシオ計算機株式会社 測定装置、測定方法及び測定プログラム
JP2016112053A (ja) * 2014-12-11 2016-06-23 国立研究開発法人産業技術総合研究所 歩行状態判定方法、プログラム、及び装置
JP2023549625A (ja) * 2020-08-28 2023-11-29 ローベルト ボツシユ ゲゼルシヤフト ミツト ベシユレンクテル ハフツング 水泳ストロークを決定するためのコントローラ及び方法
JP7521115B2 (ja) 2020-08-28 2024-07-23 ローベルト ボツシユ ゲゼルシヤフト ミツト ベシユレンクテル ハフツング 水泳ストロークを決定するためのコントローラ及び方法

Also Published As

Publication number Publication date
JPWO2012036135A1 (ja) 2014-02-03

Similar Documents

Publication Publication Date Title
US11298036B2 (en) Wearable device including PPG and inertial sensors for assessing physical activity and biometric parameters
CN102548474B (zh) 针对体动信号的信息处理方法、信息处理系统及信息处理装置
JP6403696B2 (ja) 身体的活動のモニタリングデバイス及びその方法
EP2255209B1 (fr) Procédé et appareil pour déterminer la position de fixation d'un appareil de détection de mouvement
WO2012036135A1 (fr) Procédé de traitement de l'information, dispositif de traitement de l'information, dispositif de sortie, système de traitement de l'information, programme de traitement de l'information et support d'enregistrement lisible par ordinateur sur lequel le même programme est enregistré
US11504030B2 (en) Method and system for heterogeneous event detection
JP6771162B2 (ja) 認知機能評価装置、認知機能評価方法及びプログラム
US10264997B1 (en) Systems and methods for selecting accelerometer data to store on computer-readable media
JP5659644B2 (ja) 情報処理方法，情報処理システム，情報処理装置，情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体
CN114341947B (zh) 用于使用可穿戴设备的锻炼类型辨识的系统和方法
JP4504071B2 (ja) 身体運動解析装置
Stamatakis et al. Gait feature extraction in Parkinson's disease using low-cost accelerometers
KR20110127856A (ko) 운동패턴 분석방법 및 이를 이용한 운동량 산출장치
JP5742146B2 (ja) 情報処理方法，情報処理システム，情報処理装置，表示装置，情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2020120807A (ja) 転倒リスク評価装置、転倒リスク評価方法及び転倒リスク評価プログラム
CN103759738A (zh) 一种计步器
JP6638860B2 (ja) 情報処理システム、情報処理装置、および情報処理方法
JP6629485B2 (ja) 被検者の歩行運動を決定する処理装置及び方法
WO2020049621A1 (fr) Programme de détermination de statut de déambulation, procédé de détermination statut de déambulation, et dispositif de traitement d'informations
Kuruvithadam An Intelligent In-Shoe System for Real-Time Gait Monitoring and Analysis
Christ et al. An approach for determining linear velocities of athletes from acceleration measurements using a neural network
Sher Automating gait analysis using a smartphone
Ye et al. Inertial sensor based post fall analysis for false alarming reduction
HK1248507B (zh) 确定人类运动活动类型的方法及其实施装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11825137

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2012534001

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 11825137

Country of ref document: EP

Kind code of ref document: A1