WO2013171799A1 - 生体リズム推定装置 - Google Patents

生体リズム推定装置 Download PDF

Info

Publication number
WO2013171799A1
WO2013171799A1 PCT/JP2012/003258 JP2012003258W WO2013171799A1 WO 2013171799 A1 WO2013171799 A1 WO 2013171799A1 JP 2012003258 W JP2012003258 W JP 2012003258W WO 2013171799 A1 WO2013171799 A1 WO 2013171799A1
Authority
WO
WIPO (PCT)
Prior art keywords
period
biological rhythm
data
detection system
sleep
Prior art date
Application number
PCT/JP2012/003258
Other languages
English (en)
French (fr)
Inventor
健 河本
栗山　裕之
美蘭 小笠原
奈津美 佐藤
Original Assignee
株式会社日立製作所
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 株式会社日立製作所 filed Critical 株式会社日立製作所
Priority to PCT/JP2012/003258 priority Critical patent/WO2013171799A1/ja
Priority to JP2014515355A priority patent/JPWO2013171799A1/ja
Publication of WO2013171799A1 publication Critical patent/WO2013171799A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B10/00—Other methods or instruments for diagnosis, e.g. instruments for taking a cell sample, for biopsy, for vaccination diagnosis; Sex determination; Ovulation-period determination; Throat striking implements
    • A61B10/0012—Ovulation-period determination
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording 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, mobility of a limb
    • A61B5/1118—Determining activity level
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/48—Other medical applications
    • A61B5/4857—Indicating the phase of biorhythm
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A61B5/681—Wristwatch-type devices
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches

Definitions

  • the present invention relates to a measurement method for measuring a subject's biological rhythm for a long period of time, an output device for outputting the subject's biological rhythm, and a control method for the output device.
  • Detecting the menstrual cycle of women has high needs for health management and QOL improvement of a wide range of women. For example, predicting ovulation dates is useful information for those who want children and for contraception, especially for those who are prone to disruption of the menstrual cycle. Also become. Furthermore, since appropriate makeup and skin care methods differ depending on the stage of the menstrual cycle, if the menstrual cycle can be known accurately, the cosmetic effect of the makeup can be optimized by itself. Menopause, which plagues many women, is known to have irregular physiology as a sign, and it may be possible to reduce symptoms by early intervention such as hormonal treatment by detecting irregular physiology.
  • Patent Document 1 Based on the above background, for example, a technique based on pulse measurement during sleep as in Patent Document 1 has been proposed as means for easily detecting the menstrual cycle. This is because the average pulse rate during sleep increases with the onset of menstruation, and the pulse during sleep is measured by a mattress-type pulse wave sensor. This is a method of presenting a day with a high average pulse rate as a menstruation start date candidate.
  • Patent Document 1 The technique based on pulse measurement represented by Patent Document 1 cannot be used by everyone because it lacks low cost, reliability, and minimal invasiveness. This is because each of the following current pulse measurement methods has its own simplicity and invasiveness:
  • Mattress type sensor A mattress where a person sleeps is filled with liquid, and a pressure sensor is used to extract minute vibrations such as human body movements and pulse waves. It is expensive if it is a large-scale device that is expensive and does not sleep on the bed it is installed on (such as someone else's house or sleeping on the sofa watching TV), or if there are multiple people on the mattress. There was a problem that it could not be done. For accurate detection of the menstrual cycle, it is important to be able to measure without missing daily conditions, and therefore, the number of subjects who can effectively use the mattress sensor for menstruation is limited.
  • Chest electrocardiograph An electrocardiograph that is wound around the chest mainly for running, etc. is on the market, but because it is used during exercise, it has a strong feeling of pressure and is uncomfortable due to stuffiness. It is difficult to put on.
  • Wristband type heart rate monitor Although the mainstream method is to measure the amount of infrared reflection from the vein, it requires a lot of time to charge because it consumes more electricity than an electrocardiograph, and it cannot be measured unless it is in strong contact with the arm. High degree of invasiveness.
  • Patent Document 2 No simple means for measuring respiratory rate is known for Patent Document 2, and how the respiratory rate during sleep specifically does not have been known in the past. It has never been possible to detect the period.
  • pulse measurement does not satisfy the need for anyone to detect menstruation easily, inexpensively, with high reliability, and with minimal invasiveness.
  • a bracelet-type acceleration sensor that is worn on the wrist is known as a means for simply measuring the sleep state in daily life.
  • the bracelet-type acceleration sensor is very inexpensive, has a small burden on the body, and in recent years has become widespread as a sleep meter that can measure an activity meter and sleep / wake rhythm.
  • accelerometers have traditionally measured body movement, that is, walking and turning over, it has not been possible to detect breathing or menstrual cycles using bracelet-type acceleration sensors.
  • a body motion sensor that acquires sensing data indicating body movement, a target period calculation unit that calculates an analysis target period from the sensing data, and sensing data for each predetermined period within the analysis target period
  • a main period component detecting unit for detecting main period components
  • a representative value calculating unit for calculating a main period component representative value from a group of a plurality of main period components calculated within the target period
  • a calculated main value There is a biological rhythm detection system including a biological rhythm detection unit that detects a biological rhythm period from a periodic component representative value.
  • the menstrual cycle can be accurately detected from data obtained from an acceleration sensor worn on the wrist.
  • the example of the whole figure which shows an example of a structure of the menstrual cycle detection system to which this invention is applied.
  • mold sensor node 1 which comprises the life visualization system of a present Example.
  • the example of the flowchart figure which shows the whole flow of a data process.
  • Explanatory drawing which shows the example of a format of a respiration estimation data table.
  • the present inventors have been researching the bracelet-type acceleration sensor disclosed in Non-Patent Document 1, and the low-frequency component of the acceleration data measured by the bracelet-type acceleration sensor includes the trunk of the trunk derived from respiratory motion. It was discovered that periodic movements may be reflected and the respiratory rate during sleep can be estimated from a bracelet-type acceleration sensor. Furthermore, since we discovered a specific tendency that the average respiratory rate during sleep rises immediately after ovulation and falls shortly before the start of menstruation, we can detect the menstrual cycle by measuring long-term fluctuations in respiratory rate. showed that. With these two findings, we have realized a technology that detects the menstrual cycle simply, inexpensively, with high reliability, and with minimal invasiveness using a bracelet-type acceleration sensor. This is shown below.
  • FIG. 17A shows a frequency spectrum of a scalar value of an acceleration waveform in actual sleep data of less than 7 minutes. In this example, 0.22 Hz (one The case where there is a strong peak in the respiration frequency of 13.2 times per minute) is shown. As shown in FIG. 18A, this phenomenon occurs when the wrist on which the bracelet type acceleration sensor 1 is mounted is in a position where it is affected by the respiratory motion of the trunk. Therefore, for example, when the wrist is separated from the trunk as shown in FIG. 18B, the respiration cycle component does not appear in the acceleration data as shown in FIG. 17B.
  • Patent Document 1 discloses means for estimating a pulse and a respiratory rate from a mattress type sensor.
  • breathing motion is reflected in the sensor data regardless of the posture of the body except during the body motion period. It is possible to accurately estimate the respiratory frequency simply by finding the frequency component.
  • the influence of noise (temperature change, sensor measurement noise, etc.) on the sensor value is dominant, and simply finding the main frequency component has no meaning not related to breathing. The value will be calculated.
  • FIG. 19 (A) shows the actual respiratory frequency data for one minute during the night measured by the nostril breath sensor
  • FIG. 19 (B) shows the bracelet type acceleration sensor worn on the wrist during the same period as the breath sensor.
  • FIG. 20 (A) shows a scatter diagram in which the breath sensor and the wristband type acceleration sensor are measured simultaneously for 21 people each night, and the actually measured respiration frequency and the estimated respiration frequency per minute are calculated. In this scatter diagram, it can be seen that the estimated respiratory frequency and the measured respiratory frequency vary greatly from each other.
  • an index is defined that evaluates that the power value of the main frequency component is sufficiently prominent compared with the power value of the surrounding frequency region, and this exceeds the threshold value (ie, sufficiently).
  • the detected peak is used only when the peak is a strong peak, and the respiratory frequency is accurately detected from the acceleration data of the wrist where the respiratory motion may not be reflected by adopting the non-respiratory detection otherwise.
  • FIG. 19C although there are many missing values, it is possible to estimate the respiration frequency following the measured values with high accuracy.
  • FIG. 20B shows that the respiration frequency estimated by the proposed method can closely approximate the actually measured respiration frequency even in the data of 21 persons.
  • Fig. 21 (A) shows the average daily respiratory rate data obtained as a result of a woman wearing a bracelet type accelerometer during sleep for a long period of time, and the recorded start date of menstruation is indicated by a vertical line. ing.
  • the respiratory rate for each day varies widely, but by passing through a band-pass filter with a period of 20 days or more and 40 days or less as a band pass, data is extracted along the period starting date as shown in FIG. can do. Since the normal menstrual cycle is in the range of 25 to 38 days, this phenomenon can be reproduced by using a filter including this range.
  • the menstrual cycle detection method based on pulse data disclosed in Patent Document 1 or the like cannot be estimated with high accuracy.
  • the respiratory rate during sleep strongly varies depending on the room temperature in the bed environment and the mental and physical condition of the person, so it varies more than the pulse rate. Therefore, in a method in which a day that exceeds two days or more from the average is set as the start date of menstruation, false detection increases in respiration.
  • a day that exceeds two days or more from the average is set as the start date of menstruation, false detection increases in respiration.
  • the respiratory rate tends to decrease during cold periods and increase during summer. Therefore, in comparison with the average value of last month as in the above-mentioned document, since the baseline of respiratory rate itself has changed, accurate detection cannot be performed.
  • FIG. 22 shows the menstrual cycle estimated using the above method and the actual menstrual cycle.
  • the vertical axis shows the length of the menstrual cycle and represents how many days have passed since the previous menstruation start date and the next menstruation began.
  • the data indicated by the asterisk is the menstrual cycle estimated by the above method, and the data indicated by the circle indicates the menstrual cycle based on the person's record.
  • the average absolute error is within 1.57 days, although there is a gap depending on the month.
  • the acceleration data is measured by a sensor attached to a person (user), the respiratory frequency during sleep is calculated, the start date of menstruation is detected, and the menstruation presented to the user A period detection system is shown.
  • FIG. 1 is an overall view showing an example of the configuration of a menstrual cycle detection system to which the present invention is applied.
  • the respiration frequency estimation system of the present embodiment uses the bracelet type sensor node 1 having an acceleration sensor as a sensor for detecting the operation (or state) of the user of the system and detects the acceleration of the arm as biological information. Indicates.
  • the bracelet type sensor node 1 is mounted on the user's arm and detects acceleration, and the acceleration data detected at a predetermined cycle is transmitted as sensing data to the base station 102 via the antenna 101 and is transmitted to the client computer (PC). ) To 103. If wired communication is possible, the data is transmitted directly to the PC 103 via a USB connection or the like. Regardless of the method of sending the sensed data sequentially regardless of wireless or wired, a method of sending the data accumulated in the sensor node at regular intervals may be used.
  • the PC 103 communicates with one or a plurality of bracelet type sensor nodes 1, receives sensing data corresponding to the movement of the user from each bracelet type sensor node 1, analyzes the received sensing data, and displays it. Output data.
  • the output display data can be browsed by the client computer (PC) 103 operated by the user.
  • FIG. 2 is a diagram showing an example of the bracelet type sensor node 1 constituting the life visualization system of the present embodiment
  • FIG. 2 (A) is a schematic view seen from the front of the bracelet type sensor node 1
  • FIG. ) Is a cross-sectional view of the bracelet type sensor node 1 viewed from the side.
  • This bracelet type sensor node 1 mainly measures the movement of the user.
  • the bracelet type sensor node 1 includes a case 11 for storing a sensor and a control device, and a band 12 for mounting the case 11 on a human arm.
  • a substrate 10 including a microcomputer (not shown), a sensor 6, etc. is stored as shown in FIG.
  • a sensor 6 that measures the movement of the human body (living body)
  • an acceleration sensor that measures the three-axis accelerations XYZ in the drawing is shown.
  • the bracelet type sensor node 1 may be provided with a pressure sensor or a capacitance sensor (not shown), and whether or not the user wears the wristwatch type sensor node 1 may be output as a wearing state.
  • FIG. 3 shows a block diagram of an electronic circuit attached to the substrate 10 of the bracelet type sensor node 1.
  • a substrate 10 includes a wireless communication unit (RF) 2 including an antenna 5 that communicates with a base station 102, a USB communication unit 39 that is wired to a PC 103, a sensor 6 that is an acceleration sensor, and wireless communication.
  • a microcomputer 3 that controls the unit 2; a real-time clock (RTC) 4 that functions as a timer for intermittently starting the microcomputer 3; a battery 7 that supplies power to each unit; and a power supply to the sensor 6
  • a switch 8 is disposed to control the above. Further, a bypass capacitor C1 is connected between the switch 8 and the sensor 6 to eliminate noise and reduce charge / discharge speed to prevent wasteful power consumption.
  • the microcomputer 3 includes a CPU 34 that performs arithmetic processing, a ROM 33 that stores programs to be executed by the CPU 34, a RAM 32 that stores data, and an interrupt that interrupts the CPU 34 based on a signal (timer interrupt) from the RTC 4.
  • a serial communication interface (SCI) 36 that transmits and receives signals as serial signals between the control unit 35, the A / D converter 31 that converts the analog signal output from the sensor 6 into a digital signal, and the wireless communication unit 2.
  • a parallel interface (PIO) 37 that controls the wireless communication unit 2 and the USB communication unit 39 and the switch 8, and an oscillation unit (OSC) 30 that supplies a clock to each of the units in the microcomputer 3. .
  • PIO parallel interface
  • OSC oscillation unit
  • the above-described units in the microcomputer 3 are connected via a system bus 38.
  • the RTC 4 outputs an interrupt signal (timer interrupt) at a predetermined period set in advance in the interrupt control unit 35 of the microcomputer 3 and outputs a reference clock to the SCI 36.
  • the PIO 37 controls ON / OFF of the switch 8 in accordance with a command from the CPU 34 and controls power supply to the sensor 6.
  • the bracelet type sensor node 1 starts the microcomputer 3 at a predetermined period (for example, 1 second), acquires sensing data from the sensor 6, and an identifier for identifying the bracelet type sensor node 1 in the acquired sensing data. A time stamp is given and transmitted to the base station 102.
  • FIG. 4 is a block diagram showing each component of the menstrual cycle detection system to which the present invention shown in FIG. 1 is applied.
  • the sensing data transmitted by the bracelet type sensor node 1 is accumulated in the sensing data table 1150 of the recording device 1100 of the client (PC) 103 via the base station 102. Or you may communicate with PC103 directly via wired communication which is not illustrated.
  • the PC 103 includes a display device 1031 that displays various types of information and an input device 1032 that allows various information to be input by user operations.
  • the display device 1031 may be a display terminal such as a liquid crystal display or a CRT display, or a printer or an image file output.
  • the input device 1032 is an input device such as a keyboard and a mouse.
  • the display device 1031 and the input device 1032 may be a single device having both functions, such as a touch panel display.
  • the PC 103 includes a processor 107, a memory 108, and a recording device 1100.
  • the recording device 1100 records various programs and various data tables, which will be described later, such as a hard disk drive, a CD-ROM drive, and a flash memory. Various programs and various data tables may be divided and recorded in a plurality of recording devices.
  • the processor 107 implements various functions by reading various programs recorded in the recording device 1100 into the memory 108 and executing them. Specifically, by executing the data totaling program 200, the sensing data measured by the acceleration sensor of the user's arm is totaled, a total value for each unit time (for example, one minute) is calculated, and the recording device 1100 is stored in the aggregate data table 250. Also, by executing the sleep period extraction program 300, the calculated total value for each unit time is analyzed, and all sleep periods are detected and stored in the sleep period data table 350 of the recording device 1100.
  • the respiration frequency estimation program 400 the respiration frequency per unit time (for example, one minute) is estimated from the detected sleep period, and stored in the respiration frequency data table 450 of the recording device 1100. Also, by executing the menstrual cycle calculation program 500, the menstrual cycle is estimated from the detected respiration frequency and stored in the menstrual cycle data table 550 of the recording device 1100.
  • the PC 103 executes the data aggregation program 200, the sleep period extraction program 300, the respiration frequency estimation program 400, and the menstrual cycle estimation program 500 based on a periodic cycle or communication with the wristwatch-type sensor node 1. Further, an example is shown in which display data is presented on the display device 1031 due to the operation of the input device 1032, the activation of the PC 103, or the end of the execution of the menstrual cycle estimation program 500.
  • FIG. 5 is a flowchart showing the overall flow of data processing performed in the system to which the present invention is applied.
  • step S 1 the base station 102 transfers the sensing data transmitted from the bracelet type sensor node 1 to the PC 103 and accumulates the sensing data in the sensing data table 1150 of the PC 103. Further, an identifier for identifying the sensor node given to the sensing data and time information indicating the time when the sensing data is acquired are also stored in the sensing data table. Further, the PC 103 executes the data totaling program 200, calculates the exercise frequency per unit time from the sensing data accumulated in the recording device 1100, and stores it in the totaling data table 250 of the recording device 1100.
  • the data aggregation program 200 may be executed every predetermined period (for example, 5 minutes), may be executed due to the start or end of communication with the bracelet type sensor node 1, or the input device 1032 The operation may be performed as a cause.
  • step S ⁇ b> 2 the PC 103 executes the sleep period extraction program 300 to detect an area where the user is estimated to be in a sleep state from the aggregate data stored in the aggregate data table 250, and all sleep areas Are stored in the sleep analysis data table 350 as a set. Further, the main sleep which is the maximum sleep period of the day is detected, and the main sleep tag is added to the stored sleep region which is the main sleep, and stored in the sleep analysis data table 350.
  • the sleep period extraction program 300 may be executed every predetermined cycle (for example, 5 minutes), or may be executed due to the end of the data aggregation program 200.
  • the PC 103 executes the respiration frequency estimation program 400 to acquire sensing data within the period from the sensing data table 1150 for each sleep period stored in the sleep data table 350, and unit time (for example, , 1 second), a frequency spectrum for a predetermined period (for example, 5 minutes) is calculated, a frequency component having the maximum power is detected as a respiratory frequency candidate, and the maximum power value is sufficiently protruding. (60 * frequency) is determined as the respiratory frequency for the minute of the date and time, and if not sufficiently protruding, a tag is added as the date and time of the respiratory date not detected and stored in the respiratory frequency data table 450.
  • the respiration frequency estimation program 400 may be executed every predetermined cycle (for example, 5 minutes), or may be executed due to the end of the sleep period extraction program 300.
  • step S4 the PC 103 executes the menstrual cycle estimation program 400, calculates the average value of the respiratory frequency during the main sleep for all the days for which the main sleep is calculated, and the time series of the average respiratory rate for each day.
  • Periodic data corresponding to the menstrual cycle is extracted by applying a bandpass filter having a passband period of 20 days to 40 days to the data, and the menstrual period is detected by detecting a zero crossing point from a positive value to a negative value of the periodic data.
  • the start date candidate is calculated, and the ovulation day candidate is calculated by detecting the zero crossing point from the negative value to the positive value of the cycle data, and stored in the menstrual cycle data table 550.
  • the menstrual cycle estimation program 500 may be executed every predetermined cycle (for example, one day), or may be executed due to the end of the respiration frequency estimation program 400.
  • step S4 the respiration frequency stored in the respiration frequency estimation program 400 is presented on the display unit 1031 of the client computer 103.
  • FIG. 6 is a flowchart showing an example of processing performed by the data totaling program 200 of the PC 103.
  • sensing data corresponding to the sensor identifier held by the user is read from the sensing data table 1150.
  • the identifier of the sensor owned by the user may be acquired from, for example, the wristwatch type sensor node 1 communicating with the base station 102, or may be an identifier designated by the user with the input device 102, It may be an arbitrary identifier selected from the user sensor correspondence table.
  • the amount of sensing data to be read here is a predetermined period (for example, 5 minutes) that is an aggregation period of sensing data, or everything after the last aggregation time already stored by the execution of the past data aggregation program 200, etc. You only have to set it.
  • a total value is calculated for each predetermined time interval (for example, 1 minute) for the acceleration data of the read sensing data.
  • the number of zero crosses indicating the frequency of exercise of the wearer (user) of the bracelet type sensor node 1 within a predetermined time interval is used as the total value.
  • the sensing data detected by the bracelet type sensor node 1 includes X, Y, and Z-axis acceleration data
  • the scalar amount of the X, Y, and Z-axis acceleration ⁇ (X ⁇ 2 + Y ⁇ 2 + Z ⁇ 2) is calculated (step S12), and the obtained scalar quantity is filtered (bandpass filter) to extract only a predetermined frequency band (for example, 0.1 Hz to 5 Hz) and remove noise components (step). S13).
  • a value at which the obtained scalar amount passes a predetermined threshold for example, 0.05 G
  • the frequency at which the number of zero crossings appears within a predetermined time interval is calculated.
  • the appearance frequency is output as an exercise frequency at a predetermined time interval (1 minute) (step S14).
  • the result of calculating the exercise frequency is data obtained by sorting the exercise frequency for each unit time in time series as shown in FIG.
  • the motion frequency may be other methods, such as counting the number of times that the acceleration value in each direction of X, Y, and Z vibrates positively and negatively (frequency) within a predetermined time in each direction. In this embodiment, since the calculation can be simplified, a method of calculating the number of zero crossings is employed.
  • a flag representing the state of data within a predetermined time interval is calculated (step S15).
  • the ratio of valid data that is, data in which three-axis data of X, Y, and Z within a predetermined range exist as values other than missing values
  • the time interval is determined as missing data.
  • a threshold for example, 0.8
  • the time interval is determined as non-wearing data.
  • the flag of the data state at the time interval is determined as one of “data present”, “missing value”, and “non-attached”.
  • Execution of the data totaling program 200 obtains an exercise frequency and a data flag for each predetermined time interval, generates total data for each predetermined time interval, as shown in FIG. It accumulates together with the identifier of the wearer of the bracelet type sensor node 1 and the identifier of the bracelet type sensor node 1 (step S16).
  • FIG. 9 is an explanatory diagram showing the format of the aggregate data table 250.
  • the user ID 251 stores the identifier of the wearer of the bracelet type sensor node 1 (user of the menstrual cycle detection system), and the bracelet type sensor included in the sensing data.
  • a sensor data ID 252 that stores the identifier of the node 1
  • a measurement date and time 253 that stores a start time (measurement date and time) of a predetermined time interval
  • an exercise frequency 254 that stores an exercise frequency calculated by the execution of the data aggregation program 200
  • One entry is made up of a flag 255 that stores a status flag of data obtained by executing the data totaling program 200.
  • the user identifier may be referred to from a table (not shown) set in advance based on the identifier of the bracelet type sensor node 1.
  • FIG. 10 is a flowchart showing an example of processing performed by the sleep period extraction program 300 of the PC 103.
  • the unit time aggregate data aggregated by the execution of the data aggregation program 200 is read from the aggregate data table 250 (step S21).
  • the amount of the total data read here may be set to all after the end time of the last sleep period already stored by the execution of the past sleep analysis program 300, for example.
  • a period group estimated to be in a sleep state is detected from the read aggregated data.
  • the frequency of exercise during sleep is extremely low, the human body does exercise such as turning over during sleep, so the frequency of exercise does not become zero.
  • Several methods for determining sleep are known. For example, the Cole method (Non-patent Document 2) may be applied. The start time and end time of each period detected by such a method are held as a sleep period candidate group in a temporary storage (not shown) or the like.
  • the exercise frequency is close to zero and may be determined as sleep. For example, when a wristwatch-type sensor node has a non-wearing determination unit, sleep is not possible. It is also possible to prevent this by treating the data indicating that the flag 258 is not attached when reading the unit time totalization data table 250 as equivalent to data with high exercise frequency.
  • step S23 the adjacent sleep period candidates are combined.
  • the sleep period candidates are divided at the time of the temporary wakeup.
  • the next sleep area candidate starts within a predetermined time (for example, 30 minutes) after the end of the sleep period candidate, the two sleep period candidates are combined and handled as one large sleep period.
  • a predetermined time for example, 30 minutes
  • step S24 those that are incompatible as sleep period candidates are excluded.
  • a sleep period whose duration is a predetermined time (for example, 10 minutes) or less is excluded. If the end time of the last sleep period candidate in the candidate group is within a predetermined time (for example, 30 minutes) from the latest measurement time of the total data read by execution of the sleep analysis program 300, the next sleep analysis program This is also excluded because there is a possibility that it can be combined with a sleep period newly listed as a candidate in 300 executions.
  • the sleep period candidate group processed as described above is determined as the sleep period group.
  • step S25 main sleep is extracted from the sleep period group determined in step S24, the sleep type is determined as "main sleep", and the other sleep periods are determined as "nap”.
  • the calendar date to which each sleep period belongs is calculated. This means that if the end time of the sleep region is a predetermined time, for example, from 0 o'clock to 20 o'clock, it belongs to the same day, and if it is from 20 o'clock to 24 o'clock, it belongs to the next day. This standard is because it can be considered that sleep that ends before 20 o'clock in the ordinary life is included in the nap.
  • the sleep region that starts at 17:00 on July 23 and ends at 19:30 belongs to July 23, for example, the sleep region that starts at 16:30 on July 23 and ends at 20:30 It belongs to July 24th. From the earliest calendar date calculated in this way to the latest calendar date, the longest belonging sleep area in each calendar day is derived, and these are determined as “main sleep” on that calendar day. The type of sleep other than the “main sleep” calculated above is determined as “snapping”.
  • step S26 the confirmed sleep period group is accumulated in the sleep analysis data table 350 of the recording device 1100 as shown in FIG.
  • an identifier and sleep ID that are unique in the sleep analysis data table 350 are assigned to each sleep period. This may be selected, for example, by using a value obtained by adding 1 to the sleep ID assigned last.
  • FIG. 11 is an explanatory diagram showing the format of the sleep period data table 350.
  • the user ID 351 stores the identifier of the wearer of the bracelet type sensor node 1
  • the sleep ID 352 stores the sleep identifier
  • the start time of the sleep period One entry is configured from the start date and time 353 to be stored, the end time 354 to store the end time of the sleep period, and the sleep type 355 to store the type of sleep (whether it is main sleep or nap).
  • FIG. 12 is a flowchart showing an example of processing performed by the respiration estimation program 400 of the PC 103.
  • step S31 the sleep period extracted by the sleep period extraction program 300 is read from the sleep period data table 350.
  • the amount of the sleep period read here may be set to, for example, all after the date and time of the last respiratory data already stored by executing the past respiratory estimation program 400.
  • step S32 to step S37 each sleep period acquired here is individually processed.
  • step S32 for one period of the sleep period acquired in step S31, all sensor data corresponding to the identifier of the wristwatch sensor node 1 worn by the user included in the period are stored in the sensing data table. 1150 is read.
  • step S33 sensing data for a predetermined period (for example, 5 minutes) are cut out and scalarized for each unit time (for example, 1 minute) from the sensor data acquired in step S32.
  • the scalar amount of the acceleration data of the X, Y, and Z axes of the sensing data detected by the bracelet type sensor node 1 ⁇ (X ⁇ 2 + Y ⁇ 2 + Z ⁇ 2) is calculated.
  • step S34 the obtained scalar quantity is filtered (bandpass filter) to extract only a predetermined frequency band (for example, 0.01 Hz or more and 1 Hz or less) to remove noise components.
  • the pass band of the filter used at this time may be anything as long as it includes an average respiration frequency band during human sleep. Generally, human respiration is 0.16 Hz to 0.5 Hz. It may be 0.01 Hz or more and 1 Hz or less.
  • step S35 a frequency spectrum is obtained for the scalar quantity filtered in the previous step S34.
  • FFT Fast Fourier Transform
  • step S36 the frequency having the maximum intensity is acquired as the main frequency from the frequency spectrum obtained in the previous step S35.
  • step S37 the validity of the main frequency obtained in the previous step S36 as a respiratory frequency is verified.
  • the main frequency is within a predetermined frequency range (for example, 0.016 Hz or more and 0.33 Hz or less). Otherwise, it may be excluded as being too early or too late for breathing. good.
  • a predetermined frequency range for example, 0.016 Hz or more and 0.33 Hz or less.
  • the possibility of noise may be high and may be excluded.
  • Equation 1 the degree to which the intensity of the main frequency (peak) protrudes from other frequencies is evaluated. If the degree of protrusion is small, the possibility of noise is high and may be excluded. As a means for evaluating the degree of protrusion, for example, the following peak score (Equation 1) may be used. (Equation 1)
  • this peak score is below a threshold (for example, 8.0), it is evaluated that the degree of protrusion is insufficient.
  • step S38 the respiratory frequency per minute is calculated by multiplying the main frequency calculated in the previous step S6 by 60.
  • step S39 the respiration frequency for each minute at each date and time within each sleep period calculated above is accumulated in the respiration estimation data table 450 of the recording device 1100.
  • the respiration frequency holds a value indicating non-detection (for example, “null”), and otherwise, step S38. Holds the respiration frequency for each minute calculated in.
  • the example of detecting the main frequency component using FFT Fast Fourier Transform
  • any method that can detect the intensity of the frequency component included in the acceleration data may be used.
  • Autocorrelation may be obtained.
  • the validity as the respiration frequency in step S37 may include, for example, a case where the correlation coefficient at the selected ⁇ is greater than or equal to a threshold value.
  • one main frequency component is detected, its validity is verified, and when it is valid, the example is adopted as the respiratory frequency.
  • a plurality of protruding frequency components are detected. May be.
  • the validity of each is verified. If there are a plurality of valid frequency components, the most appropriate frequency component may be selected, for example. You may select the frequency component nearest to the respiration frequency detected at the date and time before the said date and time.
  • the average respiratory frequency (breathing trend) for each elapsed time after the start of bedtime is calculated by averaging the respiratory rate calculated for the wearer so far. The elapsed time after bedtime start may be calculated, and the frequency component closest to the respiratory trend may be selected.
  • sensing data for a predetermined period for example, 5 minutes
  • unit time for example, 1 minute
  • the respiration frequency is detected, and accumulated in the respiration estimation data table 450 is shown.
  • the sensing data for a predetermined period for example, 5 minutes
  • the respiration frequency is detected, and the unit for each predetermined unit time (for example, 1 minute).
  • sensing data may be filled with an average value of the period or treated as a missing value for a period during which the body movement is performed.
  • Sensing data may be filled with an average value of the period or treated as a missing value for a period during which the body movement is performed.
  • a means for detecting body movements other than the movements derived from breathing for example, when the movement variance of a scalar value (a variance value for a predetermined period (for example, 5 seconds) from each time point) is a predetermined threshold value or more, For example, motion may be detected. By removing such a large body motion, it is possible to prevent the high-frequency component generated by the motion with a strong body motion from obstructing the detection of the main frequency of respiration.
  • the main frequency is calculated by converting the sensing data into a scalar.
  • the X axis, the Y axis, and the Z axis are handled separately, and have the highest frequency component at each time point.
  • FIG. 13 is an explanatory diagram showing the format of the respiration estimation data table 450.
  • the user ID 451 for storing the identifier of the wearer of the bracelet type sensor node 1 and the date and time 452 for storing the date and time of a predetermined time interval are estimated. If the respiration frequency or a valid respiration is not detected, a respiration frequency 453 for storing a value (for example, null) indicating non-detection is held.
  • FIG. 14 is a flowchart showing an example of processing performed by the menstrual cycle estimation program 500 of the PC 103.
  • step S41 the main sleep period calculated by the sleep period extraction program 300 is acquired for each day included in the calculation target period, and the respiration frequency calculated in each main sleep period is acquired.
  • the period to be calculated may be set to, for example, all after the last day that has already been processed by the execution of the past menstrual cycle estimation program 500.
  • the representative respiration frequency for that day is calculated.
  • the representative respiratory frequency may be, for example, an average value of respiratory frequencies calculated during the main sleep of the day.
  • a predetermined time for example, 30 minutes
  • a predetermined time for example, 1 hour
  • the average value and standard deviation instead of using all values when calculating the average value, for example, calculate the average value and standard deviation, and use only the respiratory frequency of (average value-standard deviation) or more and (average value + standard deviation) or less.
  • the influence of noise may be reduced by calculating an average value.
  • the day when respiratory detection was not performed normally is, for example, the day when the respiratory rate is detected below a predetermined number (for example, 10) or the average respiratory rate is abnormal compared to the average respiratory rate for the past week. If it is too high. The latter is because a sudden poor physical condition such as a cold is suspected and has nothing to do with the menstrual cycle to be detected.
  • the respiratory representative value for each day is calculated as described above.
  • step S43 the representative respiratory frequencies calculated for each day are arranged in time series, and bandpass filter processing is performed.
  • the filter used here is preferably a bandpass filter that passes a period of at least 23 days and not more than 33 days as a passband. Thus, periodic data representing the menstrual cycle is obtained.
  • step S44 and step S45 the date when the calculated periodic data changes from a positive value to a negative value is calculated as the estimated menstrual start date, and the date when the calculated periodic data changes from a negative value to a positive value is calculated as the estimated ovulation date. . Further, the number of days elapsed from the previous menstrual start date to the current menstrual start date is calculated as a menstrual cycle.
  • step S39 the estimated menstrual start date, the estimated ovulation date, and the menstrual cycle calculated above are accumulated in the menstrual cycle data table 550 of the recording device 1100.
  • the estimated date of menstruation closest to each of the start dates is searched for, and the actual menstruation starts
  • the individual difference between the estimated value and the actually measured value Can be minimized.
  • this method may have a slightly larger error depending on the season. Therefore, for example, by subtracting a predetermined correction value determined for each month from the estimated value, the seasonal difference from the actual measurement value can be minimized.
  • the user's actual menstrual start date may be input, but there are other significant changes in daytime behavior, for example, near the estimated date
  • the day can also be used as the actual start date of menstruation.
  • women will have to stand more often in the bathroom.
  • the bracelet type acceleration sensor assumed in the embodiment can detect walking behavior, for example, the behavior is performed in the order of ⁇ a walking period of 3 minutes or less, a non-walking period of 10 minutes or less, a walking period of 3 minutes or less>.
  • the menstrual start date and ovulation date may be predicted for a predetermined number of times in the future (for example, three times), and stored in the menstrual cycle data table 550 in step S46.
  • FIG. 15 is an explanatory diagram showing the format of the menstrual cycle data table 550.
  • the user ID 551 stores the identifier of the wearer of the bracelet type sensor node 1, the estimated menstrual start date 552, the estimated ovulation date 553, and the like.
  • the menstrual cycle 554 for accumulating the elapsed days from the previous menstrual start date is held.
  • the ovulation date should be stored in the same record as the date of menstruation detected after the ovulation date.
  • FIG. 16 is a screen image of the menstrual cycle display screen 1200 displayed on the display unit 1031 of the PC 103 which is a client computer.
  • the display of the menstrual cycle display screen 1200 by the PC 103 may be caused by accepting a display request from a user via the input device 1032 or may be caused by the end of execution of the menstrual cycle estimation program 500.
  • the sleep period extraction program 300 may detect that the wearer has woken up.
  • a browser may be adopted, or an application that runs alone may directly display the menstrual cycle display screen 1200.
  • the menstrual cycle display screen 1200 is an example of a screen that presents information about the detected menstrual cycle to the user.
  • This includes a menstrual cycle graph 1201 that displays menstrual cycle data for a period to be displayed, a scheduled menstrual start date label 1207 that indicates the next scheduled start date of menstruation, a scheduled ovulation date label 1208 that indicates the next scheduled date of ovulation, , A pregnancy possibility label 1209 indicating the possibility of pregnancy as of today, and an advice label 1210 for displaying advice based on the phase in today's menstrual cycle.
  • the period to be displayed may be, for example, a period specified by a date period selection control (not shown) or a predetermined number of days before and after this day (for example, one month).
  • the menstrual cycle graph 1201 is a graph that displays menstrual cycle data of a period to be displayed, menstruation start date and ovulation date estimated within the period.
  • the period data calculated by the menstrual period detection program 500 is drawn on the graph, the past data is drawn with a solid line (1202), and the future part is displayed with a broken line (1203), so that the current position in the menstrual cycle is easily displayed. I can do it.
  • by drawing the estimated menstrual start date and ovulation date within the period with vertical lines (1204, 1205) it is possible to easily look back on past menstrual periods and grasp future plans.
  • the actual menstrual start date may be input by dragging the vertical line 1204 with an input device such as a mouse.
  • the actual menstrual start date is stored in a data table (not shown), and may be used for estimation of the estimated menstrual date in the processing of the menstrual detection program 500, for example.
  • the menstruation start schedule label 1207 displays the calendar date estimated next based on today's date
  • the ovulation scheduled date label 1208 shows the next estimated menstruation start based on today's date.
  • the pregnancy possibility label 1209 calculates the possibility of pregnancy from the current phase in the menstrual cycle and shows it in easy-to-understand expressions such as “low” and “high”, making it easy to grasp the ease of pregnancy today. To do.
  • the possibility of pregnancy is displayed as “high” if it is within a predetermined number of days (for example, 3 days) from any ovulation estimation date accumulated in the menstrual cycle data table 550, and “low” otherwise. Is displayed.
  • the advice label 1210 displays advice based on today's phase in the menstrual cycle. For example, since it is generally known that after two months from the start of menstruation and until the day of ovulation, the mood is clear and exercise and study progress, the next scheduled date of menstruation is later than the next scheduled date of ovulation. In the case of, by displaying a message such as "I feel refreshed today and is suitable for exercise and studying", I encourage those actions at a time suitable for exercise and studying I can do it.
  • FIG. 17 is a screen image of the long-term cycle display screen 1300 displayed on the display unit 1031 of the PC 103 which is a client computer.
  • the display of the long-term cycle display screen 1300 by the PC 103 may be caused by accepting a display request from the user via the input device 1032 or may be caused by the end of execution of the menstrual cycle estimation program 500.
  • the sleep period extraction program 300 may detect that the wearer has woken up.
  • a browser may be adopted, or an application that runs alone may directly display the long-term cycle display screen 1300.
  • the long-term cycle display screen 1300 is an example of a screen that presents long-term information about a detected menstrual cycle (the number of days from the start of menstruation to the start of the next menstruation) to the user.
  • This is composed of a menstrual cycle graph 1301 for displaying menstrual cycle data of a period to be displayed, and a regularity graph 1302 for displaying regularity of the menstrual cycle.
  • the period to be displayed may be, for example, a period specified by a date period selection control (not shown) or the entire period in which user data exists.
  • the menstrual cycle graph 1301 is a graph that displays menstrual cycle data for a period to be displayed. By arranging the menstrual cycle data 554 accumulated in the menstrual cycle data table 550 in time series, it is possible to visually grasp the tendency of the person compared with the average menstrual cycle of 28 days.
  • the regularity graph 1302 is a graph that displays the local regularity of the menstrual cycle data in the period to be displayed. By plotting the moving variance value of the data displayed on the menstrual cycle graph 1301, it is possible to visually grasp the time when the menstrual cycle is disturbed and the regular time. For example, it is known that the menstrual cycle is disturbed before menopause, and various climacteric disorders occur accordingly. One of the pains of menopause is that it is difficult to improve the symptoms. By looking at the regularity graph 1302 and confirming that the disorder of regularity gradually converges, it is useful for having hope for the convergence of menopause.
  • the three-axis acceleration sensor of the bracelet type sensor node 1 attached to the user's arm is used to measure the menstrual cycle of the user (human body) as the menstrual cycle estimation system.
  • Any sensor that can detect the activity state of the human body in a non-invasive manner may be used.
  • an angular velocity sensor attached to an arm may be used, and a biaxial or non-axial acceleration sensor may be used instead of a triaxial acceleration sensor.
  • the system for estimating the respiratory frequency during sleep is exemplified, but it is also effective to estimate the respiratory frequency during awakening rather than during sleep.
  • the breathing frequency component may be reflected in the acceleration sensor even when arms are folded during the meeting or the breath is rising immediately after exercise, but this can also be detected by the same means as described above.
  • the number of breathing during awakening varies greatly depending on the behavior and time of day. For example, in a user's life, the same behavior is always performed at the same time every day (for example, commuting, sitting on the train at the same time every day, etc.) ) And the processing of the menstrual cycle detection program 500 is performed for that period, so that the menstrual cycle can be detected even if data during sleep cannot be measured.
  • an acceleration sensor attached to the arm is used as an example of a sensor in which the state of respiration is reflected in a fragmentary manner.
  • the invention can effectively measure the respiratory frequency and detect the menstrual cycle.
  • the present invention can be applied because the movement of breathing is not reflected on the video depending on the sleeping posture.
  • an acceleration sensor or a face-to-face sensor has been mounted on a name tag badge worn by an employee in a company. Although the breathing movement is not reflected in the acceleration data from the attached badge, the breathing movement may be reflected depending on the sitting posture. Can be estimated.
  • the user of the present invention does not need to be a human being, and the present invention can be effectively applied to a living organism having periodic estrus periods such as dogs and cats.
  • the estrus period comes in a half-year cycle, so it is desirable to use a band pass having a half year as a pass band in the menstrual cycle detection program 500.
  • the gender of the sensor user can be estimated.
  • the menstrual cycle is detected by applying the present invention to long-term sensor data, it is estimated that the woman is a premenopausal woman. Therefore, whether men or women can be accurately estimated in the age before menopause, depending on whether or not menstrual cycle detection is performed normally.
  • pregnancy can be easily determined by using the present invention.
  • the advice panel 1210 of the menstrual cycle display screen 1200 may indicate “There is a possibility of pregnancy. By displaying a comment such as “Please,” you can promptly confirm pregnancy.
  • the menstrual cycle can be estimated accurately and easily from output data of a sensor in which respiratory motion is not reflected in the data in many cases, such as an acceleration sensor worn on an arm.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Physiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Abstract

　簡便に，安価に，信頼性高く，低侵襲で月経の検知を行うことを課題とする。上記課題を解決する一例として、体の動きを示すセンシングデータを取得する体動センサと，センシングデータから解析対象期間を算出する対象期間算出部と，解析対象期間内の所定の期間毎にセンシングデータを集計し，主周期成分を検出する主周期成分検出部と，対象期間内に算出された複数の主周期成分の群から主周期成分代表値を算出する代表値算出部と，算出された主周期成分代表値から生体リズムの周期を検出する生体リズム検出部と，を備えることを特徴とする生体リズム検出システムがある。

Description

生体リズム推定装置
　本発明は，対象者の生体リズムを長期的に測定する測定方法，該対象者の生体リズムを出力する出力装置，出力装置の制御方法に関するものである。
　女性の月経周期を検知することは，幅広い層の女性の健康管理やＱＯＬ向上のために高いニーズがある。例えば子どもを欲す人にとっても避妊のためにも，排卵日の予測は有用な情報であり，特に月経周期が乱れがちな人にとっては月経開始日の予測は外出のスケジュールや準備のための指針にもなる。更に月経周期の段階によって適切な化粧や肌のケア方法が違うため，月経周期を正確に知る事が出来れば化粧の美容効果を自ら最適化することも可能となる。また多くの女性を悩ませる更年期障害は，兆候として生理不順が訪れることが知られており，生理不順の検知によってホルモン治療等の早期介入によって症状を軽減することも考えられる。
　このように，簡便な月経周期の検知が実現すれば幅広い層が恩恵を受けることが出来る。しかし，誰にでも利用できるためには，検知手段が安価で，利用者に要求される判断や行動が少なく，信頼性が高く，侵襲性が軽微である事が求められる。
　以上の背景に基づき，簡便に月経周期を検知する手段として例えば特許文献１のような睡眠中の脈拍測定に基づく手法が提案されている。これは，月経の開始と共に睡眠中の平均脈拍数が上昇するという観察から，マットレス型の脈波センサにより睡眠中の脈拍を測定し，前回の月経以降の平均呼吸数よりも２日間以上睡眠中の平均脈拍数が高かった日を月経開始日候補として提示する手法である。
　また特許文献２で指摘されるように，脈拍以外では，睡眠中の呼吸数も月経周期に連動することが知られている。
特開２００６－９４９６９号公報 特表第２００９－５３２０７２号公報 特許第２００７－１４３６４２号公報
Ｔａｎａｋａ　Ｔ，　ｅｔ　ａｌ．，　"Ｌｉｆｅ　Ｍｉｃｒｏｓｃｏｐｅ：Ｃｏｎｔｉｎｕｏｕｓ　ｄａｉｌｙ－ａｃｔｉｖｉｔｙ　ｒｅｃｏｒｄｉｎｇ　ｓｙｓｔｅｍ　ｗｉｔｈ　ｔｉｎｙ　ｗｉｒｅｌｅｓｓ　ｓｅｎｓｏｒ"，　５ｔｈ　Ｉｎｔｅｒｎａｔｉｏｎａｌ　Ｃｏｎｆｅｒｅｎｅ　ｏｎ　Ｎｅｔｗｏｒｋｅｄ　Ｓｅｎｓｉｎｇ　Ｓｙｓｔｅｍｓ，　２００８ Ｃｏｌｅ　ＲＪ，Ｋｒｉｐｋｅ　ＤＦ，　Ｇｒｕｅｎ　Ｗ，　Ｍｕｌｌａｎｅｙ　ＤＪ，　Ｇｉｌｌｉｎ　ＪＣ．　Ａｕｔｏｍａｔｉｃ　ｓｌｅｅｐ／ｗａｋｅ　ｉｄｅｎｔｉｆｉｃａｔｉｏｎ　ｆｒｏｍ　ｗｒｉｓｔ　ａｃｔｉｖｉｔｙ．　Ｓｌｅｅｐ　１９９２；　１５：　４９１－４６９．
　これまでは基礎体温の測定のような，負荷の大きい検知手段しか普及していなかった。
　前記特許文献１に代表される脈拍測定に基づく技術は，安価性，信頼性，低侵襲性に欠けるため万人が用いる事は出来なかった。これは，現在存在する以下のような脈拍の測定手段はそれぞれ簡便性や侵襲性に課題を抱えるためである：
　マットレス型センサ：　人が寝るマットレスの内部に液体が満たされており，圧力センサにより人の体動や脈波などの微細な振動を抽出するセンサである。大掛かりな装置で高額であり，設置されているベッドで寝なかった場合（他人の家や，テレビを見ながらソファで寝てしまった場合など）や，マットレス上に複数人が居る場合は測定が出来ないという問題があった。月経周期の正確な検知のためには毎日の状態が欠けることなく測定出来る事は重要であるため，マットレス型センサを効果的に月経検知のために利用出来る対象者は限られている。
　胸部心電計：　主にランニング用などで胸部に巻く心電計が市場に出回っているが，運動中に使うものであるため圧迫感が強く，蒸れるため不快であり，長期的に毎日睡眠中に装着することは困難である。
　リストバンド型心拍計：　静脈の赤外線反射量を測定する方式が主流だが，心電計と比較して電気消費量が大きいため充電の手間が大きく，強く腕に接触していないと測定できないために侵襲度が高い。
　特許文献２に対しても，呼吸数の簡便な測定手段は知られていなく，また睡眠中の呼吸数が具体的にどのように振る舞うかは従来知られていないため，呼吸数を用いて月経周期を検知することは従来可能ではなかった。
　以上のように，誰もが簡便に，安価に，信頼性高く，低侵襲で月経の検知を行うというニーズは脈拍測定では満たされていない。
　ところで，日常生活において簡便に睡眠の状態を測定する手段として，手首に装着する腕輪型の加速度センサが知られている。腕輪型の加速度センサは非常に安価であり，身体への負担が少なく，また近年では活動量計や睡眠・覚醒のリズムが測定できる睡眠計として広く普及しつつある。しかし加速度計は従来あくまで体動，つまり歩行や寝返りなどを測定するものであったため，腕輪型の加速度センサを用いて呼吸や月経周期の検知を行うことはこれまで出来なかった。
　上記課題を解決する一例として，体の動きを示すセンシングデータを取得する体動センサと，センシングデータから解析対象期間を算出する対象期間算出部と，解析対象期間内の所定の期間毎にセンシングデータを集計し，主周期成分を検出する主周期成分検出部と，対象期間内に算出された複数の主周期成分の群から主周期成分代表値を算出する代表値算出部と，算出された主周期成分代表値から生体リズムの周期を検出する生体リズム検出部と，を備えることを特徴とする生体リズム検出システムがある。
　本発明によれば，手首に装着する加速度センサから得られるデータからも，精度良く月経周期を検知することができる。
本発明を適用する月経周期検知システムの構成の一例を示す全体図の例。 本実施例の生活可視化システムを構成する腕輪型センサノード１の一例を示す図の例。 腕輪型センサノードの基板に取り付けられた電子回路のブロック図の例。 本発明を適用する月経周期検知システムの各構成要素を示すブロック図の例。 データ処理の全体的な流れを示すフローチャート図の例。 ＰＣのデータ集計プログラムで行われる処理の一例を示すフローチャート図。 スカラー量が所定の閾値を通過する値をゼロクロス回数として算出する例を示した図。 単位時間毎の運動頻度を時系列的にソートしたデータの例を示した図。 集計データテーブルのフォーマットの例を示す図。 ＰＣの睡眠期間抽出プログラムで行われる処理の一例を示すフローチャート図。 睡眠期間データテーブルのフォーマットの例を示す説明図。 ＰＣの呼吸推定プログラムで行われる処理の一例を示すフローチャート図。 呼吸推定データテーブルのフォーマットの例を示す説明図。 ＰＣの月経周期推定プログラムで行われる処理の一例を示すフローチャート。 月経周期データテーブルのフォーマットの例を示す説明図。 ＰＣの表示部に表示される月経周期表示画面の画面イメージ例を示した図。 睡眠データにおける加速度波形のスカラー値の周波数スペクトルの例を示した図。 腕輪型加速度センサを装着した手首の位置と体幹の呼吸運動の影響との関係の例を示した図。 一晩中の一分間毎の呼吸頻度データの例を示した図。 一晩ずつについて呼気センサとリストバンド型加速度センサを同時測定し，一分毎の実測呼吸頻度と推定呼吸頻度を算出した散布図の例。 ある女性が腕輪型加速度センサを長期間睡眠中に装着した結果得られた，日毎の平均呼吸数データを点で示し，記録された月経の開始日を縦線にて示した例の図。 本実施例にて推定した月経周期と実測の月経周期の例を示した図。
　以下，本発明の実施の形態を，図面に基づき説明する。
　本発明者らは，非特許文献１に開示される腕輪型の加速度センサの研究を行っており，腕輪型の加速度センサで測定した加速度データの低周波成分に，呼吸運動に由来する体幹の周期的な動作が反映されている場合があり，睡眠中の呼吸数を腕輪型加速度センサから推定し得ることを発見した。更に，睡眠中の平均的な呼吸数は排卵直後に上昇し，月経開始の少し前に下降するという具体的な傾向を発見したため，呼吸数の長期変動を測定することで月経周期を検知出来る事を示した。これら２つの発見により，我々は腕輪型の加速度センサにより月経周期を簡便に，安価に，信頼性高く，低侵襲で検知する技術を実現した。以下にその示す。
　　（１）腕輪型加速度センサデータに基づく呼吸数推定
　図１７（Ａ）では実際の７分間弱の睡眠データにおける加速度波形のスカラー値の周波数スペクトルが図示されており，この例では０．２２Ｈｚ（一分間の呼吸頻度１３．２回に相当）に強いピークがある場合が示されている。この現象は，図１８（Ａ）に示すように，腕輪型加速度センサ１を装着した手首が，体幹の呼吸運動の影響を受けるような位置にある場合に生じている。そのため，例えば図１８（Ｂ）のように手首が体幹から離れてしまっている場合は図１７（Ｂ）のように，呼吸の周期成分は加速度データの中に現れない。
　特許文献１ではマットレス型センサから脈拍や呼吸数を推定する手段が開示されている。マットレス型のセンサでは体動期間中以外は体の姿勢に関わらず呼吸運動はセンサデータに反映されるため，当該文献で開示されるように，体動が検知された箇所以外のデータから主となる周波数成分を見つけるだけでも呼吸頻度を精度良く推定する事が出来る。しかし，腕から測定したデータでは，呼吸の動作がデータに反映されている場合のほうが稀である。そのため，多くの場合はセンサの値に乗るノイズ（温度変化や，センサの計測ノイズ等）の影響が優位となり，単純に主となる周波数成分を見つけるだけでは，呼吸とは関わりのない意味のない値が算出されてしまう。
　この事を図１９にて実例で示す。図１９（Ａ）は鼻口の呼気センサにより実測した，一晩中の実際の一分間毎の呼吸頻度データであり，図１９（Ｂ）は呼気センサと同期間手首に装着した腕輪型加速度センサのデータに対して，上記特許文献１で開示される手法を適用させた場合に算出される，一分間毎の推定呼吸頻度である。実測データと同じ値が算出される箇所もあるものの，全く異なる値も多く入り混じっており，この推定結果からは，対象者の睡眠中の平均呼吸数を正確に知る事は困難である。図２０（Ａ）に，２１人の一晩ずつについて呼気センサとリストバンド型加速度センサを同時測定し，一分毎の実測呼吸頻度と推定呼吸頻度を算出した散布図を示す。この散布図では，推定呼吸頻度と実測呼吸頻度は互いに大きくばらついている事が分かる。
　この問題に対しては，例えば主周波数成分のパワー値が周辺の周波数領域のパワー値と比較して充分に突出していることを評価する指標を定義し，これが閾値を超えた（即ち，充分に強いピークであった）場合にのみ，検知したピークを採用し，それ以外の場合は呼吸不検知とすることにより，呼吸運動が反映されない場合が存在する手首の加速度データからも，精度良く呼吸頻度を推定する事が出来る。この方式によると，図１９（Ｃ）に例示するように，欠損値を多く含むものの，実測値を精度良く追う呼吸頻度の推定を行う事が出来る。図２０（Ｂ）では，２１人のデータにおいても，提案方式で推定した呼吸頻度が実測呼吸頻度をよく近似出来ている事を示している。
　　（２）呼吸数に基づく月経周期検知
　我々は，腕輪型加速度センサで測定した長期的な平均呼吸数には月経周期と一致する周期性があることを発見した。図２１（Ａ）はある女性が腕輪型加速度センサを長期間睡眠中に装着した結果得られた，日毎の平均呼吸数データを点で示し，記録された月経の開始日を縦線にて示している。日毎の呼吸数は大きくばらついているが，２０日以上４０日以下の周期をバンドパスとするバンドパスフィルタを通すことで，図２１（Ｃ）のように周期が月経開始日に沿うデータを抽出することができる。正常な月経周期は２５～３８日の範囲であるため，この範囲を含むフィルタを用いればこの現象を再現出来る。
　従来特許文献１などで開示されていた脈拍データに基づく月経周期検知手法を呼吸データに適用するだけでは，精度良く月経周期を推定することが出来ない。まず，睡眠中の呼吸数は寝床環境の室内温や，本人の心身の状態に強く依存するため脈拍と比較してばらつきが大きい。そのため平均より２日間以上超えた日を月経開始日とするような方法では呼吸では誤検知が多くなってしまう。また脈拍には季節性変動がないが，呼吸数は寒い時期に少なくなり，夏には高くなるという傾向があることを我々は発見した。そのため，前記文献のように先月の平均値との比較では，呼吸数のベースライン自体が変わっているため正確な検知は出来ない。最後に，脈拍は月経の開始と共に上昇する傾向を示すが，呼吸はそのようには振る舞わない。図２１（Ｃ）で示すように，呼吸数は月経の開始の数日手前には既にピークを迎えており，最大の落ち幅の時に月経が開始する。我々が初めて具体的に確認した当該事実に基づかなくては，呼吸数に基づき月経周期の検知を行うことは出来ない。
　図２２に，以上の手法を用いて推定した月経周期と実測の月経周期を示す。縦軸は月経周期長を示しており，前回の月経開始日から何日間が経過してから次の月経が始まったかを表わす。米印で示したデータは上記手法により推定した月経周期で，丸で示したデータは本人の記録に基づく月経周期を表す。月によってはズレがあるものの，平均の絶対誤差は１．５７日に納まっている。
　本発明を適用する第１の実施例として，人（利用者）に装着したセンサによって加速度データを測定し，睡眠中の呼吸頻度を算出し，月経開始日を検知し，利用者に提示する月経周期検知システムを示す。
　図１は，本発明を適用する月経周期検知システムの構成の一例を示す全体図である。本実施例の呼吸頻度推定システムは，当該システムの利用者の動作（または状態）を検出するセンサとして，加速度センサを備えた腕輪型センサノード１を用い，生体情報として腕の加速度を検出する例を示す。腕輪型センサノード１は，利用者の腕に装着されて加速度を検出し，所定の周期で検出した加速度データをセンシングデータとして，アンテナ１０１を介して基地局１０２へ無線送信し，クライアント計算機（ＰＣ）１０３へ送信する。また，有線通信が可能な場合はＵＳＢ接続などを介して直接ＰＣ１０３へ送信する。無線，有線を問わず，センシングしたデータを逐次的に送る方法のほか，センサノードに蓄積したデータを一定周期毎に送信する方法でもよい。
　図１において，ＰＣ１０３は一つまたは複数の腕輪型センサノード１と通信を行い，各腕輪型センサノード１から利用者の動きに応じたセンシングデータを受信し，受信したセンシングデータを解析し，表示データを出力する。出力される表示データは利用者が操作するクライアント計算機（ＰＣ）１０３で閲覧することができる。
　図２は，本実施例の生活可視化システムを構成する腕輪型センサノード１の一例を示す図で，図２（Ａ）は腕輪型センサノード１の正面から見た概略図で，図２（Ｂ）は腕輪型センサノード１を側方から見た断面図である。この腕輪型センサノード１は主に利用者の動きを測定する。
　腕輪型センサノード１は，センサや制御装置を格納するケース１１と，ケース１１を人体の腕に装着するバンド１２を備える。
　ケース１１の内部には，図２（Ｂ）のように図示しないマイクロコンピュータやセンサ６等を備えた基板１０が格納される。そして，人体（生体）の動きを測定するセンサ６としては，図中Ｘ－Ｙ－Ｚの３軸の加速度をそれぞれ測定する加速度センサを採用した例を示す。なお，腕輪型センサノード１には図示しない感圧センサや静電容量センサを備え，利用者が腕時計型センサノード１を装着しているか否かを装着状態として出力しても良い。
　図３は，腕輪型センサノード１の基板１０に取り付けられた電子回路のブロック図を示す。図３において，基板１０には，基地局１０２と通信を行うアンテナ５を備えた無線通信部（ＲＦ）２と，ＰＣ１０３と有線接続するＵＳＢ通信部３９と，加速度センサであるセンサ６及び無線通信部２を制御するマイクロコンピュータ３と，マイクロコンピュータ３を間欠的に起動するためのタイマとして機能するリアルタイムクロック（ＲＴＣ）４と，各部に電力を供給する電池７と，センサ６への電力の供給を制御するスイッチ８が配置される。また，スイッチ８とセンサ６の間には，バイパスコンデンサＣ１が接続されてノイズの除去や，充放電の速度を低減して無駄な電力消費を防ぐ。
　マイクロコンピュータ３は，演算処理を実行するＣＰＵ３４と，ＣＰＵ３４で実行するプログラムなどを格納するＲＯＭ３３と，データなどを格納するＲＡＭ３２と，ＲＴＣ４からの信号（タイマ割り込み）に基づいてＣＰＵ３４に割り込みをかける割り込み制御部３５と，センサ６から出力されたアナログ信号をデジタル信号に変換するＡ／Ｄコンバータ３１と，無線通信部２との間でシリアル信号にて信号の送受を行うシリアルコミュニケーションインターフェース（ＳＣＩ）３６と，無線通信部２及びＵＳＢ通信部３９とスイッチ８を制御するパラレルインターフェース（ＰＩＯ）３７と，マイクロコンピュータ３内の上記各部へクロックを供給する発振部（ＯＳＣ）３０とを含んで構成される。そして，マイクロコンピュータ３内の上記各部はシステムバス３８を介して接続される。ＲＴＣ４は，マイクロコンピュータ３の割り込み制御部３５に予め設定されている所定の周期で割り込み信号（タイマ割り込み）を出力し，また，ＳＣＩ３６へ基準クロックを出力する。ＰＩＯ３７はＣＰＵ３４からの指令に応じてスイッチ８のＯＮ／ＯＦＦを制御し，センサ６への電力供給を制御する。
　腕輪型センサノード１は，所定の周期（例えば，１秒等）でマイクロコンピュータ３を起動して，センサ６からセンシングデータを取得し，取得したセンシングデータに腕輪型センサノード１を特定する識別子とタイムスタンプを付与して基地局１０２へ送信する。
　図４は，図１に示した本発明を適用する月経周期検知システムの各構成要素を示すブロック図である。腕輪型センサノード１が送信したセンシングデータは，基地局１０２を介してクライアント（ＰＣ）１０３の記録装置１１００のセンシングデータテーブル１１５０に蓄積される。もしくは図示しない有線通信を介して直接ＰＣ１０３と通信してもよい。ＰＣ１０３は，各種情報を表示する表示装置１０３１と利用者の操作によって様々な情報の入力を可能とする入力装置１０３２を具備する。
　表示装置１０３１は液晶ディスプレイやＣＲＴディスプレイ等の表示端末のほか，プリンタや画像ファイル出力でも良い。入力装置１０３２はキーボード，マウス等の入力用機器である。また，表示装置１０３１と入力装置１０３２はタッチパネル式ディスプレイのような，両方の機能を備える単体の機器でも良い。またＰＣ１０３は，プロセッサ１０７と，メモリ１０８及び記録装置１１００を備える。記録装置１１００は，後述する各種プログラム，各種データテーブルを記録するものであり，例えば，ハードディスクドライブやＣＤーＲＯＭドライブ，フラッシュメモリなどである。なお，複数の記録装置に各種プログラム，各種データテーブルを分割して記録するようにしてもよい。
　プロセッサ１０７は，記録装置１１００に記録されている各種プログラムをメモリ１０８に読み出して実行することにより各種機能を実現する。具体的には，データ集計プログラム２００を実行することにより，利用者の腕の加速度センサで測定されたセンシングデータを集計し，単位時間（例えば，１分間）毎の集計値を算出し，記録装置１１００の集計データテーブル２５０に格納する。また，睡眠期間抽出プログラム３００を実行することにより，算出した単位時間毎の集計値を解析し，全ての睡眠期間を検出し，記録装置１１００の睡眠期間データテーブル３５０に格納する。また呼吸頻度推定プログラム４００を実行することにより，検出した睡眠期間から単位時間（例えば，一分間）毎の呼吸頻度を推定し，記録装置１１００の呼吸頻度データテーブル４５０に格納する。また，月経周期算出プログラム５００を実行することにより，検出した呼吸頻度から月経周期を推定し，記録装置１１００の月経周期データテーブル５５０に格納する。
　なお，以下では，ＰＣ１０３が，データ集計プログラム２００と睡眠期間抽出プログラム３００と呼吸頻度推定プログラム４００と月経周期推定プログラム５００を定期的な周期，もしくは腕時計型センサノード１との通信を起因として実行し，更に入力装置１０３２の操作，もしくはＰＣ１０３の起動，もしくは月経周期推定プログラム５００の実行の終了を起因として表示装置１０３１に表示データを提示する例を示す。
　図５は，本発明を適用するシステムで行われるデータ処理の全体的な流れを示すフローチャート図である。
　まず，ステップＳ１では，腕輪型センサノード１が送信したセンシングデータを基地局１０２がＰＣ１０３へ転送し，ＰＣ１０３のセンシングデータテーブル１１５０にセンシングデータを蓄積する。また，センシングデータに付与されている，センサノードを特定するための識別子と当該センシングデータを取得した時刻を示す時刻情報もセンシングデータテーブルに蓄積する。さらに，ＰＣ１０３がデータ集計プログラム２００を実行して，記録装置１１００に蓄積されたセンシングデータから単位時間毎の運動頻度を算出し，記録装置１１００の集計データテーブル２５０に格納する。なお，データ集計プログラム２００を所定の周期（例えば，５分間）毎に実行しても良いし，腕輪型センサノード１との通信開始や終了を起因として実行しても良いし，入力装置１０３２の操作を起因として実行しても良い。
　次に，ステップＳ２では，ＰＣ１０３は睡眠期間抽出プログラム３００を実行して，集計データテーブル２５０に格納された集計データから利用者が睡眠状態にあると推定される領域を検出し，全ての睡眠領域の開始時刻と終了時刻を組みとして睡眠分析データテーブル３５０に格納する。また，一日の最大の睡眠期間である主睡眠を検出し，上記格納した睡眠領域のうち主睡眠であるものに対して主睡眠タグを付加して睡眠分析データテーブル３５０に格納する。なお睡眠期間抽出プログラム３００を所定の周期（例えば，５分間）毎に実行しても良いし，データ集計プログラム２００の終了を起因として実行しても良い。
　次に，ステップＳ３では，ＰＣ１０３は呼吸頻度推定プログラム４００を実行して，睡眠データテーブル３５０に格納された各睡眠期間について，期間内のセンシングデータをセンシングデータテーブル１１５０から取得し，単位時間（例えば，１秒）毎の，所定の期間（例えば，５分）の周波数スペクトルを算出し，最大パワーを持つ周波数成分を呼吸周波数候補として検出し，該最大パワー値が充分に突出している場合には（６０＊周波数）を当該日時の分間呼吸頻度として確定し，充分に突出していない場合は当該日時を呼吸未検出としてタグを付加し，呼吸頻度データテーブル４５０に格納する。なお呼吸頻度推定プログラム４００を所定の周期（例えば，５分間）毎に実行しても良いし，睡眠期間抽出プログラム３００の終了を起因として実行しても良い。
　次に，ステップＳ４では，ＰＣ１０３は月経周期推定プログラム４００を実行して，主睡眠が算出された全ての日について主睡眠中の呼吸頻度の平均値を算出し，日毎の平均呼吸数の時系列データにパスバンド２０日以上４０日以下の周期を持つバンドパスフィルタを施すことで月経周期に対応する周期データを抽出し，周期データの正から負の値へゼロクロスする点を検出することで月経開始日候補を算出し，周期データの負から正の値へゼロクロスする点を検出することで排卵日候補を算出し，月経周期データテーブル５５０に格納する。なお月経周期推定プログラム５００を所定の周期（例えば，１日）毎に実行しても良いし，呼吸頻度推定プログラム４００の終了を起因として実行しても良い。
　次に，ステップＳ４では，呼吸頻度推定プログラム４００に格納された呼吸頻度をクライアント計算機１０３の表示部１０３１に提示する。
　図６は，ＰＣ１０３のデータ集計プログラム２００で行われる処理の一例を示すフローチャートである。
　まずステップＳ１１では，利用者の保有するセンサの識別子に対応するセンシングデータをセンシングデータテーブル１１５０から読み込む。ここで利用者の保有するセンサの識別子は，例えば基地局１０２と通信している腕時計型センサノード１から取得しても良いし，入力装置１０２で利用者により指定された識別子でもよいし，図示しない利用者センサ対応テーブルから選択した任意の識別子でも良い。ここで読み込むセンシングデータの量は，センシングデータの集計周期である所定の周期（例えば，５分間），あるいは過去のデータ集計プログラム２００の実行によって既に格納されている最後の集計時刻以降全て，等に設定すればよい。
　次に，ステップＳ１２～Ｓ１４では，読み込んだセンシングデータの加速度データについて所定の時間間隔（例えば，１分）毎の集計値を算出する。本実施形態では，所定の時間間隔内での腕輪型センサノード１の装着者（利用者）の運動の頻度を示すゼロクロス回数を集計値として用いる。
　腕輪型センサノード１が検出したセンシングデータにはＸ，Ｙ，Ｚ軸の加速度データが含まれているので，Ｘ，Ｙ，Ｚの３軸の加速度のスカラー量＝√（Ｘ＾２＋Ｙ＾２＋Ｚ＾２）を算出し（ステップＳ１２），求めたスカラー量をフィルタ（バンドパスフィルタ）処理することで所定の周波数帯域（例えば，０．１Ｈｚ以上５Ｈｚ以下）のみを抽出しノイズ成分を除去する（ステップＳ１３）。そして，図７に示すように求めたスカラー量が所定の閾値（例えば，０．０５Ｇ）を通過する値をゼロクロス回数として算出し，ゼロクロス回数が所定時間間隔内に出現する頻度を算出し，この出現頻度を，所定の時間間隔（１分間）の運動頻度として出力する（ステップＳ１４）。この運動頻度の算出結果は図８で示すように単位時間毎の運動頻度を時系列的にソートしたデータとなる。なお，運動頻度は，Ｘ，Ｙ，Ｚの各方向の加速度の値が正と負に振動した回数（振動数）を各方向の所定時間内に数えて合計するなど他の方法でもよいが，本実施例では，計算を簡略化することができるため，ゼロクロス回数を算出する方法を採用している。
　さらに，所定時間間隔内のデータの状態を表すフラグを算出する（ステップＳ１５）。まず，所定時間間隔内に存在するデータのうち，有効なデータ（即ち，所定の範囲内のＸ，Ｙ，Ｚの３軸のデータが欠損値以外の値として存在するもの）の割合が所定の閾値（例えば０．８）より少ない場合，該時間間隔を欠損データとして確定する。また腕時計型センサノード１が装着状態を出力する場合，該時間間隔が欠損でなく所定時間間隔内に存在するデータのうち装着中である割合が所定の閾値（例えば０．８）より少ない場合，該時間間隔を非装着データとして確定する。以上により，該時間間隔のデータの状態のフラグを「データ有り」「欠損値」「非装着」のうち一つに決定する。
　データ集計プログラム２００の実行は，所定の時間間隔毎に，運動頻度，データフラグを求め，図９に示すように，所定の時間間隔毎の集計データを生成し，記録装置１１００の集計データテーブル２５０へ腕輪型センサノード１の装着者の識別子と，腕輪型センサノード１の識別子と共に蓄積する（ステップＳ１６）。
　図９は，集計データテーブル２５０のフォーマットを示す説明図で，腕輪型センサノード１の装着者（月経周期検知システムの利用者）の識別子を格納するユーザＩＤ２５１と，センシングデータに含まれる腕輪型センサノード１の識別子を格納するセンサデータＩＤ２５２と，所定の時間間隔の開始時刻（測定日時）を格納する測定日時２５３と，データ集計プログラム２００の実行により演算した運動頻度を格納する運動頻度２５４と，データ集計プログラム２００の実行により求めたデータの状態フラグを格納するフラグ２５５からひとつのエントリを構成する。なお，利用者の識別子は，腕輪型センサノード１の識別子に基づいて予め設定した図示しないテーブルから参照すればよい。
　図１０は，ＰＣ１０３の睡眠期間抽出プログラム３００で行われる処理の一例を示すフローチャートである。まず，データ集計プログラム２００の実行により集計した単位時間集計データを集計データテーブル２５０から読み込む（ステップＳ２１）。ここで読み込む集計データの量は，例えば過去の睡眠分析プログラム３００の実行によって既に格納されている最後の睡眠期間の終了時刻以降全て，等に設定すればよい。
　次に，ステップＳ２２では，読み込んだ集計データから睡眠状態であると推定される期間群を検出する。睡眠中の運動頻度は極めて低いが，睡眠中でも人体は寝返りなどの運動を行うため，運動頻度はゼロにはならない。睡眠を判定する手法はいくつか知られており，例えば，Ｃｏｌｅ法（非特許文献２）などを適用すればよい。このような手法により検出された各期間の開始時刻と終了時刻を睡眠期間の候補群として，図示しない一時ストレージ等に保持する。また，利用者がセンサノードを装着していない場合は運動頻度がゼロに近くなり，睡眠として判定されてしまう場合もあるが，例えば腕時計型センサノードが非装着判定手段を備えている場合は睡眠として判定せず，さらに単位時間集計データテーブル２５０を読み込む際にフラグ２５８が非装着を示すデータは運動頻度が高いデータと同等に扱うことによっても，これを防ぐことが可能である。
　次に，ステップＳ２３では，睡眠期間の候補の中で近接するもの同士を結合する。睡眠検出の手法によっては，例えば目覚まし時計を止めるために一時的に起床し，再び睡眠に入った場合でも一時的に起きた時刻で睡眠期間候補が区切られてしまう。しかし，生理活動としての睡眠は確かにそこで区切れているが，生活行動としての睡眠を考える場合そこで区切れることは望ましくない。そのため，睡眠期間候補の終了後，所定時間（例えば３０分）以内に次の睡眠領域候補が開始している場合，二つの睡眠期間候補を結合し，一つの大きな睡眠期間として扱う。このようにして睡眠期間候補群の中で結合できるものを探索し，結合する。
　次に，ステップＳ２４では，睡眠期間の候補として不適合であるものを排除する。一例としてあげる方法においては，まず継続時間が所定時間（例えば１０分）以下の睡眠期間は排除する。また，候補群の中の最後の睡眠期間候補の終了時刻がもし睡眠分析プログラム３００の実行により読み込んだ集計データの最新測定時刻から所定時間（例えば３０分）以内である場合，次回の睡眠分析プログラム３００の実行で新しく候補に挙がる睡眠期間と結合できる可能性があるため，これも排除する。以上のように処理した睡眠期間候補群を睡眠期間群として確定する。
　次に，ステップＳ２５では，ステップＳ２４で確定した睡眠期間群の中から，主睡眠を抽出し，睡眠の種類を「主睡眠」として確定し，それ以外の睡眠期間を「うたた寝」として確定する。まず各睡眠期間の所属するカレンダー日を算出する。これは，睡眠領域の終了時刻が所定時刻，例えば０時から２０時までであるならば同日，２０時から２４時までであるならば次の日に所属するものとする。この基準は，一般人の生活において２０時手前に終わる睡眠は昼寝に含まれると考える事も出来るからである。そして，例えば７月２３日の１７時に開始し１９時半に終了する睡眠領域は７月２３日に属し，例えば７月２３日の１６時半に開始し２０時半に終了する睡眠領域は，７月２４日に属する。こうして算出したカレンダー日の最古のカレンダー日から最新のカレンダー日まで，各カレンダー日において最長の所属睡眠領域を導き，これらをそのカレンダー日の「主睡眠」として確定する。以上で算出した「主睡眠」以外の睡眠の種類を「うたた寝」として確定する。
　最後に，ステップＳ２６では，確定した睡眠期間群を図１１に示すように，記録装置１１００の睡眠分析データテーブル３５０に蓄積する。この時，各睡眠期間に睡眠分析データテーブル３５０内で一意となるような識別子，睡眠ＩＤを割り当てる。これは例えば最後に割り当てられた睡眠ＩＤに１を足した値を利用する，等の選び方でよい。
　図１１は，睡眠期間データテーブル３５０のフォーマットを示す説明図で，腕輪型センサノード１の装着者の識別子を格納するユーザＩＤ３５１と，睡眠の識別子を格納する睡眠ＩＤ３５２と，睡眠期間の開始時刻を格納する開始日時３５３と，睡眠期間の終了時刻を格納する終了時刻３５４と，睡眠の種類（主睡眠であるか，うたた寝であるか）を格納する睡眠の種類３５５からひとつのエントリを構成する。
　図１２は，ＰＣ１０３の呼吸推定プログラム４００で行われる処理の一例を示すフローチャートである。
　ステップＳ３１では，睡眠期間抽出プログラム３００で抽出された睡眠期間を睡眠期間データテーブル３５０から読み込む。ここで読み込む睡眠期間の量は，例えば過去の呼吸推定プログラム４００の実行によって既に格納されている最後の呼吸データの日時以降全て，等に設定すればよい。
　次に，ステップＳ３２～ステップＳ３７では，ここで取得した各睡眠期間について個別に処理を行う。
　ステップＳ３２では，前記ステップＳ３１で取得した睡眠期間のうち一つの期間について，その期間内に含まれる，当該利用者が装着する腕時計型センサノード１の識別子に対応する全てのセンサデータをセンシングデータテーブル１１５０から読み込む。
　次に，ステップＳ３３では前記ステップＳ３２で取得したセンサデータについて，単位時間（例えば，１分）毎に，所定の期間（例えば，５分）のセンシングデータを切り出し，スカラー化する。腕輪型センサノード１が検出したセンシングデータのＸ，Ｙ，Ｚ軸の加速度データのスカラー量＝√（Ｘ＾２＋Ｙ＾２＋Ｚ＾２）を算出する。
　次に，ステップＳ３４では，求めたスカラー量をフィルタ（バンドパスフィルタ）処理することで所定の周波数帯域（例えば，０．０１Ｈｚ以上１Ｈｚ以下）のみを抽出しノイズ成分を除去する。この時用いるフィルタのパスバンドは人の睡眠時における平均的な呼吸の周波数帯域が含まれていればなんでもよく，一般的には人の呼吸は０．１６Ｈｚ以上０．５Ｈｚ以下であるため，例えば０．０１Ｈｚ以上１Ｈｚ以下でも良い。
　次に，ステップＳ３５では，前ステップＳ３４でフィルタ処理したスカラー量について周波数スペクトルを求める。該スカラー量についてＦＦＴ（高速フーリエ変換）を行うことで，図２１（Ａ）に例示するように，各周波数における強度を算出する。この時，各周波数の強度を，前後周波数成分の強度も含めて平均化することで，スペクトルを平滑化しても良い。
　次に，ステップＳ３６では，前ステップＳ３５で求めた周波数スペクトルから最大強度を持つ周波数を，主周波数として取得する。
　次に，ステップＳ３７では，前ステップＳ３６で求めた主周波数の，呼吸周波数としての妥当性を検証する。
　例えば，主周波数が所定の周波数領域内（例えば，０．０１６Ｈｚ以上０．３３Ｈｚ以下）であるかどうかを検証し，それ以外である場合は呼吸としては早すぎる，もしくは遅すぎるとして除外しても良い。
　また，別の例としては，主周波数の強度が所定の閾値以下であった場合には，ノイズの可能性が高いため除外しても良い。
　また，別の例としては，主周波数（ピーク）の強度が他の周波数よりも突出する度合いを評価し，突出度合いが少ない場合はノイズの可能性が高いため除外しても良い。突出する度合いを評価する手段としては，例えば下記ピークスコア（数１）を用いても良い。
（数１）
Figure JPOXMLDOC01-appb-I000001
　このピークスコアが閾値（例えば，８．０）を下回る場合は，突出度合いが足りないと評価する。
　以上で例示した条件で除外された場合は，当該日時の呼吸は「未検出」，除外されなかった場合は「検出成功」と判断する。
　次に，ステップＳ３８では，前ステップＳ６で算出した主周波数に６０をかける事で，分毎の呼吸頻度を算出する。
　最後に，ステップＳ３９では，以上で算出した，各睡眠期間内の，各日時における分毎の呼吸頻度を，記録装置１１００の呼吸推定データテーブル４５０に蓄積する。この時呼吸頻度は，ステップＳ３７で当該日時の呼吸が「未検出」であると判断された場合は，未検出を表す値（例えば，「ｎｕｌｌ」）を保持し，それ以外の場合はステップＳ３８で算出した分毎の呼吸頻度を保持する。
　また，以上の説明ではＦＦＴ（高速フーリエ変換）を用いて主となる周波数成分を検出する例を示したが，加速度データ内に含まれる周波数成分の強度が検出出来る方法であれば何でも良く，例えば自己相関を求めても良い。自己相関を用いる場合の主となる周波数成分とは，例えばτ＝０の山以降の次に現れる山の頂点としても良い。また，ステップＳ３７における呼吸周波数としての妥当性は，例えば選択されたτにおける相関係数が閾値以上である場合，などを含んでも良い。
　また，以上の説明では主となる一つの周波数成分を検出し，この妥当性を検証し，妥当である場合は呼吸頻度として採用する例を示したが，例えば複数の突出した周波数成分を検出しても良い。その場合は，例えば上記ステップＳ３７の説明のように，それぞれの妥当性を検証し，妥当であった周波数成分が複数ある場合は，例えば最も妥当であった周波数成分を選択しても良いし，当該日時の一つ手前の日時で検出された呼吸頻度に最も近い周波数成分を選択しても良い。また別の手段として，これまでに当該装着者に関して算出された睡眠期間の呼吸数を平均化することで，就寝開始後の経過時間毎の平均呼吸頻度（呼吸トレンド）を算出し，当該日時の就寝開始後経過時間を算出し，呼吸トレンドに最も近い周波数成分を選択しても良い。
　また，以上の説明では単位時間（例えば，１分）毎に，所定の期間（例えば，５分）のセンシングデータを切り出し，呼吸頻度を検出し，呼吸推定データテーブル４５０に蓄積する例を示したが，例えばより細かい間隔（例えば，１秒）毎に，所定の期間（例えば，５分）のセンシングデータを切り出し，呼吸頻度を検出し，所定の単位時間（例えば，１分）毎に，単位時間内に検出された全ての呼吸頻度を平均化することで，突発的なノイズに強くしても良い。
　また，以上の説明では単位時間毎に，所定の期間の全センシングデータを切りだして用いる例を示したが，例えば呼吸由来の運動以外の体動（例えば，寝返り）を別途検出し，このような体動が行われている期間についてはセンシングデータを当該期間の平均値で埋めてしまったり，欠損値として扱っても良い。なお，呼吸由来の運動以外の体動を検出する手段として，例えばスカラー値の移動分散（各時点から遡って所定の期間（例えば５秒間）の分散値）が所定の閾値以上である場合に体動を検出する，などでも良い。このように大きい体動を除くことで，体動の強い動作が生む高周波成分は呼吸の主周波数の検出を妨げることを未然に防ぐ事が出来る。
　また，以上の説明ではセンシングデータをスカラー化して主周波数を算出する例を示したが，例えばＸ軸，Ｙ軸，Ｚ軸をそれぞれ別に扱い，各時点において突出度の一番高い周波数成分を持つ軸を選択し，その軸の主周期を当該時点における主周波数として採用することで，特定の軸にしか呼吸運動が表れていない姿勢にいる場合でも精度よく呼吸を検出することが出来る。
　図１３は，呼吸推定データテーブル４５０のフォーマットを示す説明図で，腕輪型センサノード１の装着者の識別子を格納するユーザＩＤ４５１と，所定の時間間隔の日時を格納する日時４５２と，推定された呼吸頻度，もしくは妥当な呼吸が検出されなかった場合には未検出を表す値（例えば，ｎｕｌｌ）を格納する呼吸頻度４５３を保持する。
　図１４は，ＰＣ１０３の月経周期推定プログラム５００で行われる処理の一例を示すフローチャートである。
　ステップＳ４１では，計算対象となる期間内に含まれる各日について睡眠期間抽出プログラム３００で算出された主睡眠期間を取得し，各主睡眠期間内に算出された呼吸頻度を取得する。計算対象となる期間とは，例えば過去の月経周期推定プログラム５００の実行によって既に処理された最後の日以降以降全て，等に設定すればよい。
　次に，ステップＳ４２では，前記ステップＳ４１で呼吸頻度を取得した各日についてそれぞれ，その日の代表呼吸頻度を算出する。代表呼吸頻度とは，例えばその日の主睡眠中に算出された呼吸頻度の平均値でも良い。また，主睡眠の開始（就寝）後所定の時間（例えば３０分）から，終了（起床）前所定の時間（例えば，１時間）までの間の平均値を用いることで，飲酒した日の寝入り時の呼吸数への影響などの呼吸頻度算出の信頼性を下げる要因を排除しても良い。また平均値を算出する際全ての値を用いるのではなく，例えば平均値と標準偏差を算出し，（平均値―標準偏差）以上，（平均値＋標準偏差）以下の呼吸頻度のみを用いて平均値を算出することでノイズの影響を減らしても良い。また，呼吸検知が正常に行われなかった日は欠損日として排除しても良い。呼吸検知が正常に行われなかった日とは，例えば呼吸数が所定数（例えば，１０件）以下しか検出されなかった日や，平均呼吸数が過去一週間の平均呼吸数と比較して異常に高い場合等である。後者は，風邪等の突発的な体調不良が疑われ，本来検知したい月経周期とは無関係だからである。以上により各日の呼吸代表値を算出する。
　次に，ステップＳ４３では，前記算出した各日の代表呼吸頻度を時系列的に並べ，バンドパスフィルタ処理を行う。ここで用いるフィルタとは，少なくとも２３日以上３３日以下の周期をパスバンドとして通すバンドパスフィルタであることが望ましい。以上により，月経周期を表す周期的データを得る。
　次に，ステップＳ４４及びステップＳ４５では，前記算出した周期的データが正の値から負の値に転じる日を月経開始推定日，負の値から正の値に転じる日を排卵推定日として算出する。また前回の月経開始日から今回の月経開始日までの経過日数を月経周期として算出する。
　最後に，ステップＳ３９では，以上で算出した月経開始推定日，排卵推定日，及び月経周期を記録装置１１００の月経周期データテーブル５５０に蓄積する。
　以上の説明では，算出した周期データのゼロクロス日を月経開始日および排卵日として推定する例を示したが，実際の月経開始日と推定した月経開始日の間には数日の差がある場合がある。これは個人差と季節性変動の複合的な要因であり，それぞれを補正することでより精度よく月経開始日を推定できる場合がある。例えば図示しない手段（例えば，利用者の入力）によって過去の数回の実際の月経開始日のデータが得られる場合，それらの開始日にそれぞれ最も近い月経開始推定日を探索し，実際の月経開始日と月経開始推定日の間の日数の平均を利用者誤差として算出し，この利用者誤差を前記ステップＳ４４における月経開始推定日及び排卵推定日から差し引くことで，推定値と実測値の個人差を最小化することが出来る。また，この手法は季節によって誤差が多少大きくなる場合がある。そのため，例えば月ごとに定めた所定の補正値を推定値から差し引くことで，実測値との季節差を最小化することが出来る。
　数日程度の推定値のズレを補正するためには，利用者の実際の月経開始日の入力を受けても良いが，その他にも例えば推定日近くで日中の行動に大きな変化のあった日を実際の月経開始日として用いる事も出来る。例えば月経期間中には女性はトイレに立つ回数が増える。当該実施形態で想定する腕輪型の加速度センサであれば歩行行動を検知できるため，例えば＜３分以下の歩行期間，１０分以下の無歩行期間，３分以下の歩行期間＞の順番で行動が行われた箇所をトイレ期間として認識し，日々のトイレ期間の回数を集計し，トイレ回数が以上に高い日（例えば，過去二週間の（平均＋標準偏差）よりも高い日）を検知し，この異常日を月経開始日として用いても良い。
　また以上の説明では，全てのゼロクロス日を月経開始日もしくは排卵日として採用する例を示した。しかし，例えば生理不順によって一回分の月経周期が飛ばされた場合でも，バンドパスフィルタの特性によっては小さい山が検出され，月経開始日がご検知されてしまう可能性がある。そこで，ゼロクロス日を認識する際に，それ以前の最新の極大値の値が，過去所定の日数（例えば，半年）分の月経周期の極大値の平均と標準偏差と比較して明らかに小さい場合は月経未検知として排除しても良い。明らかに小さいとは，例えば（平均＋標準偏差／２）よりも小さい場合を含んでもいい。
　また，ステップＳ４３で算出する周期データを外挿することで，まだ実際には訪れていない未来の月経についても月経開始日及び排卵日を予測することが可能である。そのため，ステップＳ４４～Ｓ４５では未来の所定回数分（例えば，３回），月経開始日及び排卵日も予測し，ステップＳ４６で月経周期データテーブル５５０に蓄積しても良い。
　図１５は，月経周期データテーブル５５０のフォーマットを示す説明図で，腕輪型センサノード１の装着者の識別子を格納するユーザＩＤ５５１と，推定された月経開始日５５２と，推定された排卵日５５３と，前回の月経開始日からの経過日数を蓄積する月経周期５５４を保持する。排卵日は，当該排卵日の次に検出された月経開始日と同じレコードに蓄積されることが望ましい。
　図１６は，クライアント計算機であるＰＣ１０３の表示部１０３１に表示される月経周期表示画面１２００の画面イメージである。ＰＣ１０３が月経周期表示画面１２００を表示させるのは，入力装置１０３２を介した利用者からの表示要求を受け付けたことを起因としても良いし，月経周期推定プログラム５００の実行終了を起因としても良いし，例えば無線等の手段で腕輪型センサノード１のセンシングデータをリアルタイムで取得する事が可能である場合は，装着者が起床した事を睡眠期間抽出プログラム３００が検知した事を起因としても良い。なお，ＰＣ１０３で稼働するアプリケーションとしては，ブラウザを採用しても良いし，単独で稼働するアプリケーションが直接月経周期表示画面１２００を表示しても良い。
　月経周期表示画面１２００は，検知された月経周期に関する情報を利用者に提示する画面の一例である。これは，表示対象となる期間の月経周期データを表示する月経周期グラフ１２０１と，次の月経開始予定日を示す月経開始予定日ラベル１２０７と，次の排卵予定日を示す排卵予定日ラベル１２０８と，今日時点での妊娠可能性を示す妊娠可能性ラベル１２０９と，今日の月経周期における位相に基づきアドバイスを表示するアドバイスラベル１２１０から構成される。表示対象となる期間は，例えば図示しない日付期間選択コントロールによって指定された期間でも良いし，本日を含めて前後所定日数（例えば，一カ月）でも良い。
　月経周期グラフ１２０１は，表示対象となる期間の月経周期データや，当該期間内に推定された月経開始日，排卵日を表示するグラフである。月経検知プログラム５００で算出された周期データをグラフ上に，過去データは実線で描画し（１２０２），未来分については破線で表示する（１２０３）ことで，月経周期における本日の位置が分かりやすく表示出来る。また期間内に推定されている月経開始日と排卵日を縦線（１２０４，１２０５）で描画することで，過去の月経の振り返りや，未来の予定を把握しやすく出来る。また排卵前３日から排卵後２日間は妊娠の可能性が高い事が知られており，この期間を「妊娠可能日」として表示する（１２０６）ことで，妊娠を欲する人が妊娠しやすい期間を把握しやすく出来る。過去の月経開始推定日については，例えば縦線１２０４をマウス等の入力装置によりドラッグすることで実際の月経開始日を入力することを可能にしても良い。この場合，実際の月経開始日は図示しないデータテーブルに格納され，例えば月経検知プログラム５００の処理において月経推定日の推定に用いても良い。
　月経開始予定ラベル１２０７には，本日を基準にして次に推定されている月経開始日のカレンダー日を表示し，排卵予定日ラベル１２０８には，本日を基準にして次に推定されている月経開始日のカレンダー日を表示することで，近未来の体の状態を把握しやすく出来る。
　妊娠可能性ラベル１２０９には，月経周期における今日の位相から妊娠の可能性を算出し，「低い」「高い」などの分かりやすい表現で示すことで，本日の妊娠のしやすさを把握しやすくする。妊娠可能性は，例えば月経周期データテーブル５５０に蓄積された何れかの排卵推定日から今日は所定の日数（例えば３日）以内だった場合に「高い」と表示し，それ以外は「低い」と表示する。
　アドバイスラベル１２１０には，月経周期における本日の位相に基づきアドバイスを表示する。例えば月経開始以降，排卵日までの約二週間は気分がすっきりし，運動や勉強等が捗ることが一般的に知られているため，次の月経開始予定日が次の排卵予定日よりも後の場合には「今日は気分がすっきりしているため運動や勉強を行うのに適しています」等のメッセージを表示することで，運動や勉強を行うのに適した時期にそれらの行動を促す事が出来る。
　図１７は，クライアント計算機であるＰＣ１０３の表示部１０３１に表示される長期周期表示画面１３００の画面イメージである。ＰＣ１０３が長期周期表示画面１３００を表示させるのは，入力装置１０３２を介した利用者からの表示要求を受け付けたことを起因としても良いし，月経周期推定プログラム５００の実行終了を起因としても良いし，例えば無線等の手段で腕輪型センサノード１のセンシングデータをリアルタイムで取得する事が可能である場合は，装着者が起床した事を睡眠期間抽出プログラム３００が検知した事を起因としても良い。なお，ＰＣ１０３で稼働するアプリケーションとしては，ブラウザを採用しても良いし，単独で稼働するアプリケーションが直接長期周期表示画面１３００を表示しても良い。
　長期周期表示画面１３００は，検知された月経周期（月経開始日から次の月経開始までの日数）に関する長期的な情報を利用者に提示する画面の一例である。これは，表示対象となる期間の月経周期データを表示する月経周期グラフ１３０１と，月経周期の規則正しさを表示する規則性グラフ１３０２から構成される。表示対象となる期間は，例えば図示しない日付期間選択コントロールによって指定された期間でも良いし，利用者のデータの存在する全期間，などでも良い。
　月経周期グラフ１３０１は，表示対象となる期間の月経周期データを表示するグラフである。月経周期データテーブル５５０に蓄積された月経周期データ５５４を時系列的に並べることで，人の平均的な月経周期である２８日と比較した自分の傾向を視覚的に把握しやすく出来る。
　規則性グラフ１３０２は，表示対象となる期間の月経周期データの局所的な規則正しさを表示するグラフである。前記月経周期グラフ１３０１に表示したデータの移動分散値をプロットすることで，月経周期が乱れている時期と規則正しい時期を視覚的に把握しやすく出来る。例えば女性の閉経前には月経周期が乱れ，それに伴い様々な更年期障害が生じることが知られている。更年期障害の苦しみの一つは，症状の改善の目処が立ちにくいことにある。規則性グラフ１３０２を見る事で，規則性の乱れが徐々に収束していく様子を確認することで，更年期障害の収束に希望を持つこと事に役立てられる。
　なお，上記実施形態では月経周期推定システムとして利用者（人体）の月経周期を測定するために利用者の腕に装着される腕輪型センサノード１の３軸加速度センサを用いた例を示したが，非侵襲に人体の活動状態を検知可能なセンサであればよく，例えば，腕に装着した角速度センサでも良いし，３軸の加速度センサではなく，２軸や無軸の加速度センサでも良い。
　また，上記実施形態では睡眠中の呼吸頻度を推定するシステムを例示したが，睡眠中ではなく覚醒中の呼吸頻度を推定することにも効果がある。例えば会議中に腕を組んだり，運動直後に息が上がったりしている状態でも，加速度センサに呼吸の周波数成分が反映されることがあるが，これについても上記と同じ手段により検知が行える。覚醒中の呼吸数は行動内容や時間帯によって大きく変動するものであるが，例えば利用者の生活において毎日必ず同じ時間に同じ行動を行う箇所（例えば，通勤で毎日同じ時間に電車で座る，など）を検知し，その期間について月経周期検知プログラム５００の処理を行うことで，睡眠中のデータが測れなくとも月経周期の検知を行う事が出来る。
　また本発明は，呼吸の状態が断片的に反映されるセンサの例として腕に装着した加速度センサを用いた例を示したが，呼吸の状態が断片的にしか反映されないセンサであれば何でも本発明は効果的に呼吸頻度を測定し，月経周期を検知する事が出来る。例えば，寝具の上に装着した赤外線カメラで対象の動きを測定する場合は，睡眠姿勢によって呼吸の動きが映像上に反映されたりされなかったりするため，本発明が適用出来る。また，例えば特許文献３に開示されるように，近年では企業内で社員が装着する名札バッジには加速度センサや対面センサが搭載されるようになってきているが，通常時は服の上に装着したバッジからでは呼吸運動が加速度データに反映されることはないものの，座り方の姿勢によっては呼吸運動が反映される事があるため，これについても本発明を適用して呼吸頻度を精度良く推定する事が出来る。
　また，本発明の利用者は人間である必要はなく，例えば犬や猫などのように発情期を周期的に持つ生物であれば本発明は効果的に適用する事が出来る。例えば犬の場合は半年周期で発情期が訪れるため，月経周期検知プログラム５００では半年をパスバンドとして持つバンドパスを用いることが望ましい。
　また，本発明を用いることで，センサの利用者の性別を推定することも出来る。長期的なセンサデータに本発明を適用することで月経周期が検知された場合は閉経前の女性である事が推定される。そのため，月経周期検知が正常に行われるか否かで，閉経前の年代であれば男性であるか女性であるかを精度良く推定する事が出来る。
　また，本発明を用いることで，簡便に妊娠の判定を行う事が出来る。月経周期検知プログラム５００において，これまで月経が検知されていた利用者について月経が検知されなくなった場合は，例えば月経周期表示画面１２００のアドバイスパネル１２１０に「妊娠の可能性があります。医師に相談してください」等のコメントを表示することで，妊娠の確認を早期に促すこと出来る。
　上記で説明した実施形態により，腕に装着した加速度センサのように，多くの場合は呼吸運動がデータ内に反映されないセンサの出力データからも，精度良く簡便に月経周期を推定する事が出来る。
１　腕輪型センサノード
６　センサ
１０　基板
１１　ケース
１２　バンド
１０３　クライアント計算機（ＰＣ）
２５０　集計データテーブル
３５０　睡眠期間データテーブル
４５０　呼吸頻度データテーブル
１１５０　センシングデータテーブル

Claims (14)

  1. 　体の動きを示すセンシングデータを取得する体動センサと，
    　前記センシングデータから解析対象期間を算出する対象期間算出部と，
    　前記解析対象期間内の所定の期間毎に前記センシングデータを集計し，主周期成分を検出する主周期成分検出部と，
    　前記対象期間内に算出された複数の主周期成分の群から主周期成分代表値を算出する代表値算出部と，
    　前記算出された前記主周期成分代表値から生体リズムの周期を検出する生体リズム検出部と，
    　を備えることを特徴とする生体リズム検出システム。
  2. 　前記体動センサは，対象者の体に装着され、
    　前記センシングデータは、前記対象者の腕の動きを示す情報を含むことを特徴とする請求項１に記載の生体リズム検出システム。
  3. 　前記体動データは加速度データ又は角速度データを含むことを特徴とする請求項２に記載の生体リズム検出システム。
  4. 　前記対象期間算出部が算出する対象期間は，異なる日の同じ時間帯に生じる静止期間を含むことを特徴とする請求項１に記載の生体リズム検出システム。
  5. 　前記主周期成分検出部は前記センシングデータに周期成分が含まれているか判定する周期情報取得成功判定部をさらに備え，前記周期情報取得成功判定部が、前記周期成分が含まれていないと判定する場合は前記主周期成分検出部は、欠損値を出力することを特徴とする請求項１に記載の生体リズム検出システム。
  6. 　前記周期情報取得成功判定部は前記主周期成分のパワー値が予め定められた第１閾値以上である場合、前記周期成分が含まれていると判定することを特徴とする請求項５に記載の生体リズム検出システム。
  7. 　前記周期情報取得成功判定部は前記主周期成分の前記パワー値が前記主周期成分を中心とした第１の範囲内の周期成分のパワー値と比較して，所定量以上大きい場合，前記周期成分が含まれていると判定することを特徴とする請求項５に記載の生体リズム検出システム。
  8. 　前記主周期成分検出部は０．１３Ｈｚ以上０．５Ｈｚ以下の値の呼吸由来の周期成分を含めて検出することを特徴とする請求項１に記載の生体リズム検出システム。
  9. 　前記代表値算出部で算出される代表値は該対象期間内に検出された主周期成分の周波数の平均値を含むことを特徴とする請求項１に記載の生体リズム検出システム。
  10. 　前記生体リズム検出部は，前記算出した主周期に対するフィルタ処理により所定の周波数帯域を抽出し，抽出されたデータの極大点あるいは極大点の次のゼロクロス点を月経開始日として格納することを特徴とする請求項１に記載の生体リズム検出システム。
  11. 　前記生体リズム検出部は，前記算出した主周期に対する周波数解析により所定の周波数帯域を抽出し，抽出されたデータの極小点あるいは極小点の次のゼロクロス点を月経開始日として格納することを特徴とする請求項１に記載の生体リズム検出システム。
  12. 　前記生体リズム検出部は，月経／排卵開始日候補の入力を受け，前記検知した月経／排卵開始日を最寄りの開始日候補と対応させてズレ幅を算出し，計算結果をズレによって補正して出力することを特徴とする請求項１記載の生体リズム検出システム。
  13. 　前記生体リズム検出部は、前記センシングデータに基づき行動特徴量を算出し、前記算出された行動特徴量が予め定められた行動閾値より大きい場合、前記行動特徴量の算出に用いられた前記センシングデータのデータ取得日を、前記開始日候補とすることを特徴とする請求項１に記載の生体リズム検出システム。
  14. 　前記検出された生体リズムの周期に基づいて、前記センシングデータが取得された人物の性別を判定する性別判定部をさらに有することを特徴とする生体リズム検出システム。
PCT/JP2012/003258 2012-05-18 2012-05-18 生体リズム推定装置 WO2013171799A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2012/003258 WO2013171799A1 (ja) 2012-05-18 2012-05-18 生体リズム推定装置
JP2014515355A JPWO2013171799A1 (ja) 2012-05-18 2012-05-18 生体リズム推定装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/003258 WO2013171799A1 (ja) 2012-05-18 2012-05-18 生体リズム推定装置

Publications (1)

Publication Number Publication Date
WO2013171799A1 true WO2013171799A1 (ja) 2013-11-21

Family

ID=49583253

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/003258 WO2013171799A1 (ja) 2012-05-18 2012-05-18 生体リズム推定装置

Country Status (2)

Country Link
JP (1) JPWO2013171799A1 (ja)
WO (1) WO2013171799A1 (ja)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2976994A3 (en) * 2014-07-21 2016-04-27 Withings System to monitor and assist individual's sleep
JP2016131604A (ja) * 2015-01-16 2016-07-25 セイコーエプソン株式会社 生体情報測定システム、生体情報測定装置および生体情報測定方法
WO2016131630A1 (en) * 2015-02-16 2016-08-25 Ava Ag System and method for determining the fertility phase of a female
JP2017515534A (ja) * 2014-04-01 2017-06-15 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 対象の心拍を用いる目立たない排卵追跡システム及び方法
WO2018011357A3 (en) * 2016-07-13 2018-02-22 Nokia Technologies Oy Method and system to notify female fertility period
JP2018086266A (ja) * 2016-11-22 2018-06-07 パラマウントベッド株式会社 端末装置、出力方法及びコンピュータプログラム
JP2021104110A (ja) * 2019-12-26 2021-07-26 パイオニア株式会社 呼吸回数算出装置、呼吸回数算出方法、及びプログラム
WO2022085742A1 (ja) * 2020-10-21 2022-04-28 京都府公立大学法人 被検者の概日リズムの質を解析するための解析装置、解析方法、及び解析プログラム
EP4225153A4 (en) * 2020-10-07 2024-03-27 Medical Res & Development Fund For Health Services Bnai Zion Medical Center METHOD FOR DETERMINING OVULATION TIME

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006094969A (ja) * 2004-09-28 2006-04-13 Sumitomo Rubber Ind Ltd 心拍数に基づく生体リズム判定方法及び測定装置、並びに月経周期判定方法及び測定装置
JP2009532072A (ja) * 2005-11-01 2009-09-10 アーリーセンス エルティディ 臨床発作患者の監視方法及びシステム
JP2010148829A (ja) * 2008-12-26 2010-07-08 Hitachi Ltd 計算機システム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006094969A (ja) * 2004-09-28 2006-04-13 Sumitomo Rubber Ind Ltd 心拍数に基づく生体リズム判定方法及び測定装置、並びに月経周期判定方法及び測定装置
JP2009532072A (ja) * 2005-11-01 2009-09-10 アーリーセンス エルティディ 臨床発作患者の監視方法及びシステム
JP2010148829A (ja) * 2008-12-26 2010-07-08 Hitachi Ltd 計算機システム

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017515534A (ja) * 2014-04-01 2017-06-15 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 対象の心拍を用いる目立たない排卵追跡システム及び方法
US9743913B2 (en) 2014-04-01 2017-08-29 Koninklijke Philips N.V. Unobtrusive ovulation tracking system and method using a subject's heart rate
EP2976994A3 (en) * 2014-07-21 2016-04-27 Withings System to monitor and assist individual's sleep
JP2016131604A (ja) * 2015-01-16 2016-07-25 セイコーエプソン株式会社 生体情報測定システム、生体情報測定装置および生体情報測定方法
US10779802B2 (en) 2015-02-16 2020-09-22 Ava Ag System and method for determining the fertility phase of a female
WO2016131630A1 (en) * 2015-02-16 2016-08-25 Ava Ag System and method for determining the fertility phase of a female
CN107278139A (zh) * 2015-02-16 2017-10-20 阿瓦股份公司 用于确定女性的可孕期的系统和方法
WO2018011357A3 (en) * 2016-07-13 2018-02-22 Nokia Technologies Oy Method and system to notify female fertility period
JP2018086266A (ja) * 2016-11-22 2018-06-07 パラマウントベッド株式会社 端末装置、出力方法及びコンピュータプログラム
JP2021104110A (ja) * 2019-12-26 2021-07-26 パイオニア株式会社 呼吸回数算出装置、呼吸回数算出方法、及びプログラム
JP7445955B2 (ja) 2019-12-26 2024-03-08 エア・ウォーター・バイオデザイン株式会社 呼吸回数算出装置、呼吸回数算出方法、及びプログラム
EP4225153A4 (en) * 2020-10-07 2024-03-27 Medical Res & Development Fund For Health Services Bnai Zion Medical Center METHOD FOR DETERMINING OVULATION TIME
WO2022085742A1 (ja) * 2020-10-21 2022-04-28 京都府公立大学法人 被検者の概日リズムの質を解析するための解析装置、解析方法、及び解析プログラム

Also Published As

Publication number Publication date
JPWO2013171799A1 (ja) 2016-01-07

Similar Documents

Publication Publication Date Title
JP5740006B2 (ja) 呼吸測定システム及びｒｅｍ睡眠判定システム
WO2013171799A1 (ja) 生体リズム推定装置
CN108852283B (zh) 基于生理信息的睡眠评分
US10524727B2 (en) Circadian phase detection system
Chen et al. Unobtrusive sleep monitoring using smartphones
Paalasmaa et al. Unobtrusive online monitoring of sleep at home
JP5775868B2 (ja) 不眠症のための行動療法を提供するシステム及びその制御方法
US9685097B2 (en) Device and method for detecting eating activities
US20180125418A1 (en) Device and method for monitoring a physiological state of a subject
US20150245777A1 (en) Detection of emotional states
EP2032020A2 (en) Lifestyle and eating advisor based on physiological and biological rhythm monitoring
CN106073714A (zh) 一种睡眠质量的识别方法和系统
US9867597B1 (en) Method and system to notify female fertility period
JP6013668B1 (ja) 基礎体温測定システムおよび基礎体温測定装置
CN116322477A (zh) 基于温度数据的疾病检测
US10426394B2 (en) Method and apparatus for monitoring urination of a subject
CN114080180A (zh) 检测和测量打鼾
JP6813837B2 (ja) 活動リズム判定方法および活動リズム判定装置
US20220248967A1 (en) Detecting and Measuring Snoring
JP2004503282A (ja) 身体活動の測定及び分析用装置
Suzuki et al. Development of a sleep monitoring system with wearable vital sensor for home use
Jiang et al. Automatic sleep monitoring system for home healthcare
JP6702559B2 (ja) 電子機器、方法及びプログラム
Verma et al. Levels of activity identification & sleep duration detection with a wrist-worn accelerometer-based device
CN217723478U (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: 12876846

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014515355

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 12876846

Country of ref document: EP

Kind code of ref document: A1