[go: up one dir, main page]

WO2020218195A1 - 生体情報処理システムおよび生体情報取得方法 - Google Patents

生体情報処理システムおよび生体情報取得方法 Download PDF

Info

Publication number
WO2020218195A1
WO2020218195A1 PCT/JP2020/016887 JP2020016887W WO2020218195A1 WO 2020218195 A1 WO2020218195 A1 WO 2020218195A1 JP 2020016887 W JP2020016887 W JP 2020016887W WO 2020218195 A1 WO2020218195 A1 WO 2020218195A1
Authority
WO
WIPO (PCT)
Prior art keywords
peak
period
pulsation
peaks
timing
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/JP2020/016887
Other languages
English (en)
French (fr)
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.)
Sharp Corp
University Public Corporation Osaka
Original Assignee
Sharp Corp
University Public Corporation Osaka
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 Sharp Corp, University Public Corporation Osaka filed Critical Sharp Corp
Publication of WO2020218195A1 publication Critical patent/WO2020218195A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K29/00—Other apparatus for animal husbandry
    • A—HUMAN NECESSITIES
    • A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01K—ANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K67/00—Rearing or breeding animals, not otherwise provided for; New or modified breeds of animals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
    • A61B5/024—Measuring pulse rate or heart rate
    • A61B5/0245—Measuring pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/08—Measuring devices for evaluating the respiratory organs
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/346—Analysis of electrocardiograms
    • A61B5/349—Detecting specific parameters of the electrocardiograph cycle
    • A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval

Definitions

  • the following disclosure relates to technology for acquiring the mental or physical state of an organism.
  • Patent Document 1 discloses a biometric information measuring device. According to Patent Document 1, it is a biological information measuring device capable of measuring the heartbeat interval of a living body and can be worn on the living body, and is described from a plurality of electrodes in contact with the living body, a signal processing means for creating an electrocardiogram signal, and an electrocardiogram signal. A heartbeat that measures the heartbeat interval from the interval between one R wave and another R wave adjacent to the one R wave or the interval between one S wave and another S wave adjacent to the one S wave in the electrocardiogram signal.
  • An interval measuring means, a triaxial acceleration measuring means, and a temperature measuring means are provided, and the heartbeat interval obtained by the heartbeat interval measuring means, the triaxial acceleration obtained by the triaxial acceleration measuring means, and the temperature measuring means
  • a wireless transmission means for simultaneously transmitting the obtained temperature wirelessly is provided.
  • Patent Document 2 discloses a heartbeat detection method and a heartbeat detection device.
  • the heartbeat detection device includes a time difference value calculation unit that calculates a time difference value of sampling data of an electrocardiogram waveform, a time difference value determination unit that determines whether or not the time difference value exceeds a threshold value, and a time difference value determination unit. Time determination to obtain the range of the first time interval before the assumed heartbeat time, the range of the second time interval including the peak of the time difference value, and the range of the third time interval after the peak of the time difference value.
  • the relationship between the unit and the minimum value holding unit that holds the minimum value Min1, Min2, Min3 of the time difference value in the range of the first, second, and third time intervals and the minimum value Min1, Min2, Min3 determines the heartbeat time.
  • the heartbeat time determination unit is provided with the time when the time difference value exceeds the threshold value or the time when the minimum value Min2 is obtained as the heartbeat time.
  • the purpose of the present disclosure is to provide a technique for more accurately specifying the pulsation timing.
  • a beat is detected by detecting a peak of a component derived from an electrocardiogram based on a plurality of electrodes, a measuring unit that measures a potential difference between the plurality of electrodes, and a measurement result of the measuring unit.
  • a biometric information processing system including a control unit for determining dynamic timing is provided. The control unit executes the following processing for the first peak having an amplitude larger than a predetermined threshold value among the plurality of detected peaks. (1) Within the first period including the time of the first peak, three or more peaks including the first peak are specified in descending order of amplitude in the positive direction and the negative direction. (2) Within the second period shorter than the first period including the time of the first peak, one of three or more peaks in the direction opposite to the first peak is formed before and after the first peak. When it is determined that they exist one by one, the first peak is specified as the beat timing.
  • FIG. 5 is an enlarged view of electrocardiographic data in the vicinity of the target peak according to the first embodiment. It is a flowchart which shows the processing procedure for calculating the 1st autonomic nerve balance of the biological information processing system 1 which concerns on 1st Embodiment. It is a figure which shows the pulsation interval table which concerns on 1st Embodiment.
  • FIG. 5 is a Poincare plot diagram of a dog according to the first embodiment in a state in which respiration is stable in a normal state. It is a Poincare plot figure in the normal state of the dog which concerns on 1st Embodiment. It is a Poincare plot figure in the resting state of the dog which concerns on 1st Embodiment. It is a flowchart which shows the processing procedure for calculating the 2nd autonomic nerve balance of the biological information processing system 1 which concerns on 1st Embodiment.
  • FIG. 5 is an enlarged view of electrocardiographic data in the vicinity of the target peak according to the second embodiment. It is a figure which shows the functional structure of the timing acquisition part 512B which concerns on the 2nd Embodiment. It is a flowchart which shows the processing procedure of the predetermined time change processing which concerns on 2nd Embodiment.
  • FIG. 5 is a drawing showing an electrocardiographic signal and a beat timing specified based on a preferred predetermined time according to a second embodiment. It is a drawing which shows the histogram of the beat timing specified based on the predetermined time which is too short which concerns on the 2nd Embodiment.
  • FIG. 5 is a drawing showing an electrocardiographic signal and a beat timing identified based on a predetermined time that is too short, according to a second embodiment.
  • FIG. 1 is a diagram showing an overall configuration of the biometric information processing system 1 according to the present embodiment.
  • FIG. 2 is a diagram showing a functional configuration of the biometric information processing system 1 according to the present embodiment.
  • the biological information processing system 1 according to the present embodiment can also be applied to an animal having a respiratory arrhythmia. In the following, a case where the state of a dog is judged on behalf of an animal having a respiratory arrhythmia will be described.
  • the biometric information processing system 1 mainly includes electrodes 401, 402, 403 for electrocardiographic acquisition attached to the chest of a dog, and signal acquisition as a state acquisition device for processing an electrocardiographic signal.
  • the device 500 includes a communication terminal 300 as a state processing device capable of communicating with the signal acquisition device 500.
  • the electrodes 401, 402, 403 for acquiring electrocardiogram are attached to the chest or the like so as to sandwich the heart, and for example, the paws of both forefoot (or forefoot and hindfoot) have hair. It may not be in a place. Further, it is desirable that the hair is in a cut state, an electrode to which gel or the like is attached, or a structure having a protruding structure and contact with the skin even if there is hair. Alternatively, a form that induces electrocardiography through a capacitive material in a non-contact manner in the presence of hair is desirable. As a result, even an organism whose epidermis is covered with hair, such as a dog, can acquire an electrocardiogram. In the present embodiment, three electrodes 401, 402, and 403 are used, but the number of electrodes may be two or more, and more electrodes may be used. ⁇ Configuration of signal acquisition device 500>
  • FIG. 3 is a flowchart showing the entire process executed by the signal acquisition device 500 according to the present embodiment.
  • the signal acquisition device 500 includes an electrocardiographic preprocessing unit 511, a pulsation timing acquisition unit 512, and a transmission unit 560.
  • the electrocardiographic preprocessing unit 511 includes a filter and an amplifier.
  • the electrocardiographic preprocessing unit 511 converts the signal sent from the electrodes 401, 402, 403 into electrocardiographic signal data as shown in FIG. 4 and passes it to the pulsation timing acquisition unit 512 (step S002). ..
  • the electrocardiographic preprocessing unit 511 includes a filter device such as a high-pass filter and a low-pass filter, an amplification device composed of an operational amplifier, an A / D conversion device that converts an electrocardiographic analog signal into a digital signal, and the like. Is included (step S004).
  • the filter device, amplification device, and the like may be implemented by software.
  • the A / D converter it is desirable to perform sampling with a period and accuracy that can discriminate the difference in the amount of fluctuation of the beat interval. That is, it is desirable that the frequency of A / D conversion be acquired at a frequency of 25 Hz or higher.
  • the electrocardiographic signal is sampled at 100 Hz. By increasing the sampling frequency, it is possible to accurately grasp the amount of fluctuation of the beat interval.
  • the pulsation timing acquisition unit 512 is realized by, for example, a control unit including a CPU (Central Processing Unit) 510 and a memory. More specifically, the CPU 510 identifies a peak that can be determined to be an R wave among the peak signals of the electrocardiogram by executing a program stored in the memory and a determination method described later (step S006).
  • the pulsation timing acquisition unit 512 adds time information as pulsation timing to a peak that can be determined to be an R wave, and specifies the time from pulsation timing to pulsation timing as a pulsation interval (step). S008).
  • peak detection may be performed by deriving a period using an autocorrelation function, using a square wave correlation trigger, or detecting a plurality of feature points.
  • the pulsation timing acquisition unit 512 continuously specifies the pulsation timing for the continuously input electrocardiographic signals.
  • the pulsation timing acquisition unit 512 stores the time indicating the pulsation timing in the storage unit 520 for each pulsation, or transmits the time to the communication terminal 300 via the transmission unit 560 (step S010).
  • the CPU 510 has a time stamp indicating a reference time in the predetermined time zone for each predetermined time zone, and a predetermined time for each beat timing included in the predetermined time zone.
  • the detailed elapsed time from the reference time of the band is stored in the storage unit 520, or transmitted to the communication terminal 300 via the transmission unit 560.
  • the length of the predetermined time zone is 1 second
  • the detailed elapsed time from the reference time of the predetermined time zone is indicated by msec.
  • the storage unit 520 is realized by, for example, an SD card, a USB memory, or the like.
  • the transmission unit 560 is realized by a communication interface including an antenna, a connector, and the like. ⁇ How to specify the beat timing>
  • the beat timing acquisition unit 512 of the signal acquisition device 500 beats from a plurality of peaks when the electrocardiographic data with less noise as shown in FIG. 5 or the electrocardiographic data with a lot of noise as shown in FIG. 6 is obtained. Identify and extract peaks that correspond to dynamic timing.
  • the electrocardiographic signal a shown in FIG. 6 an area c in which it is difficult to specify the pulsation timing is generated due to improper attachment to an animal such as a dog or movement of the animal. Show the case. Even in such a case, there is a high possibility that the pulsation timing acquisition unit 512 according to the present embodiment can specify the accurate pulsation timing b by the following configuration and processing.
  • FIG. 7 is a block diagram showing a functional configuration of the pulsation timing acquisition unit 512 according to the present embodiment.
  • the pulsation timing acquisition unit 512 according to the present embodiment is realized by a control unit including a CPU 510 and a memory. More specifically, the CPU 510 realizes a peak detection unit 5122, a peak comparison unit 5124, a pulsation timing determination unit 5125, and the like by executing a program stored in the memory, and the storage unit 520 is the peak storage unit 5123. And beat timing storage unit 5126 and the like are realized.
  • the peak detection unit 5122 selects a target peak from the electrocardiographic data and extracts a plurality of peaks within a predetermined period based on the peak.
  • the peak storage unit 5123 stores time information and the like for identifying those peaks.
  • the peak comparison unit 514 compares a plurality of extracted peaks according to a predetermined rule.
  • the pulsation timing determination unit 5125 determines a peak satisfying a predetermined condition as the pulsation timing based on the comparison result.
  • the pulsation timing storage unit 5126 accumulates time information of the pulsation timing.
  • FIG. 8 is a flowchart showing a processing procedure of the pulsation timing specifying process according to the present embodiment.
  • the pulsation timing acquisition unit 512 acquires the electrocardiographic data
  • the pulsation timing acquisition unit 512 executes the following processing to obtain the peak in the waveform in which the electrocardiographic peak is large in the negative or positive direction. Detect and recognize the timing of beating.
  • FIG. 9 is an enlarged view of electrocardiographic data relating to the vicinity of the peak to be determined whether or not it is the pulsation timing according to the present embodiment.
  • the pulsation timing acquisition unit 512 also performs the same process for the peak protruding in the positive direction. By executing, it is configured so that the pulsation timing can be detected for the peaks protruding from both sides. It should be noted that which peak the beat timing corresponds to in the minus direction or the plus direction depends on the orientation of the electrodes attached to the dog.
  • the pulsation timing acquisition unit 512 protrudes to a threshold value V1 in the negative direction, for example, -0.2V or less with respect to the electrocardiographic waveform at a certain time point, and the product of the differential values before and after that time point is negative and immediately before. It is determined whether or not the time interval with the peak of is larger than 10 ms (step S102). If YES in step S102, the pulsation timing acquisition unit 512 determines that a peak exists at that time point. Hereinafter, the pulsation timing acquisition unit 512 determines whether or not there is a high possibility of pulsation timing with respect to the peak P1.
  • the pulsation timing acquisition unit 512 protrudes to the threshold value V1 or less in the negative direction in the first period T1 including the time point of the peak P1, and has four peaks Pm1, Pm2, Pm3, Pm4 in the positive direction in descending order of amplitude.
  • the four peaks Pp1, Pp2, Pp3, and Pp4 are identified and stored in the memory in descending order of amplitude, protruding to a threshold value V2, for example, 0.2 V or more (step S104).
  • the first period T1 is preferably, for example, 200 ms to 400 ms, that is, 400 ms to 800 ms before and after the time point of the peak.
  • the first period T1 When the first period T1 is less than 400 ms and the measurement target is a dog, noise may be detected as a peak and the heart rate may be calculated to be large. On the other hand, when the first period T1 exceeds 800 ms, adjacent beats are included in the detection period, and one of the beat timings of the beat and the adjacent beat should be detected. May detect only.
  • the pulsation timing acquisition unit 512 determines whether or not the peak Pm1 targeted this time has passed a predetermined time, for example, 200 ms or more and 10000 ms or less from the peak determined to be the previous pulsation timing. (Step S106).
  • step S106 the pulsation timing acquisition unit 512 determines whether or not the signal value Pm1 of the peak is the minimum point in the first period T1 (step S108).
  • the pulsation timing acquisition unit 512 is shorter than the first period T1 and has four peaks Pp1, Pp2, Pp3, Pp4 during the second period T2 including the peak Pm1. Of these, it is determined whether or not there is one peak before and after the peak Pm1. For example, in FIG. 9, during the second period T2, there is a peak Pp1 before the peak Pm1 and a peak Pp2 after the peak Pm1. More preferably, in the second period T2, the maximum positive peak and at least one of the 2nd to 4th peaks are present, or the 2nd positive peak and at least the 3rd to 4th peaks are present. It is determined whether or not any of the peaks is present (step S110). At this time, the second period T2 is preferably within 50 ms, preferably within 45 ms, before and after the peak Pm1.
  • the third period T3 is preferably within 20 ms before and after the peak Pm1.
  • the maximum peak Pp1 on the plus side exists in the range of 20 ms to 45 ms in the time difference of the peak Pm1 either before or after the peak Pm1 targeted this time, and the peak targeted this time. It is determined whether or not at least one of the peaks of Pp2 to Pp4 exists in the range of 20 ms to 45 ms in the time difference between before and after Pm1. If this is not the case, there is a second peak Pp2 on the plus side in the range of 20 ms to 45 ms with a time difference of either before or after the peak Pp1 targeted this time, and before the peak Pm1 targeted this time. On the other side of the above, it is determined whether or not at least one of the peaks of Pp3 and Pp4 on the plus side exists in the range of the time difference of 20 ms to 45 ms.
  • step S110 the pulsation timing acquisition unit 512 determines whether or not the amplitude of the peak Pm1 is larger than the positive peaks Pp1 and Pp2 on the plus side (step S112). ..
  • the pulsation timing acquisition unit 512 determines the timing of the peak as the pulsation timing and stores the time information indicating the timing in the storage unit 520 (step S114). ..
  • the difference (ms) from the time stamp every second is recorded as the pulsation time.
  • the signal acquisition device 500 transmits information indicating the pulsation timing for a predetermined period to the communication terminal 300.
  • the pulsation timing acquisition unit 512 identifies the accurate R wave by making determinations such as step S106, step S108, step S110, and step S112. That is, it is possible to reduce the possibility of identifying other types of peaks, body movements, electromagnetic waves, and noise due to physical impact as pulsation timing.
  • FIG. 10 is a flowchart showing a process executed by the communication terminal 300 according to the present embodiment.
  • the communication terminal 300 includes a receiving unit 361, a pulsation interval acquisition unit 321 and an analysis unit 311, a graph creating unit 312, a result output unit 313, a display 330, a data storage unit 322, and a transmitting unit 362.
  • the receiving unit 361 and the transmitting unit 362 are realized by, for example, a communication interface 360 including an antenna, a connector, and the like.
  • the receiving unit 361 receives data indicating the pulsation timing from the signal acquisition device 500 (step S100).
  • the pulsation interval acquisition unit 321 is composed of various memories 320 and the like, and stores the data received from the signal acquisition device 500.
  • the CPU 310 calculates the time between beats based on the beat timing received via the communication interface 360 (step S102).
  • the CPU 310 sequentially stores the pulsation interval as a pulsation interval table (see FIG. 11) in the memory 320 (step S104).
  • the pulsation interval is calculated in units of msec (milliseconds), for example, as shown in FIG.
  • these data may be stored in the memory 320 of the communication terminal 300, or may be stored in another device accessible from the communication terminal 300.
  • the line at the time when 9 is entered in the beat times 1 and 2 indicates a time zone in which the beat is not originally detected. This indicates a time zone in which there is no pulsation by temporarily inputting a numerical value of 9, and may be another numerical value or symbol.
  • the CPU 310 when the data indicating the pulsation timing from the signal acquisition device 500 is lost for some reason, the CPU 310 performs various plots for the period during which the pulsation interval cannot be calculated, as will be described later. Instead, various plots are restarted after the relevant period has passed. More details will be described later as a method of handling missing data.
  • the analysis unit 311, the graph creation unit 312, and the result output unit 313 are realized, for example, by the CPU 310 executing the program of the memory 320.
  • the analysis unit 311 reads the pulsation interval data from the pulsation interval acquisition unit 321 in a fixed time unit, for example, 1 minute, 10 minutes, 1 hour, etc., in a time unit necessary for determining the state, and pulsates.
  • a pulsation interval table of the interval RR (n) and the next pulsation interval RR (n + 1) is created (step S106).
  • the analysis unit 311 calculates the standard deviation of the numerical strings constituting each axis after the axis conversion (step S110).
  • the analysis unit 311 may specify the axis at which the variance is maximized by a method such as principal component analysis, and calculate the standard deviation between the axis and the axis perpendicular to the axis. Further, the analysis unit 311 may calculate the standard deviation with respect to the X-axis and the Y-axis without performing the axis conversion.
  • the directions with large variance are the X-axis direction and the Y-axis direction
  • the variation state of the beat interval plotted by Poincare is evaluated by calculating the standard deviation of the X-axis and Y-axis without performing axis conversion. it can. In this case, since it is not necessary to perform axis conversion, the amount of calculation can be reduced.
  • the biometric information processing system 1 may include a server with which the communication terminal 300 can communicate.
  • the CPU 310 as the result output unit 313 accumulates in the data storage unit 322 such as the standard deviation and the relation table, or transmits to the server via the Internet or the like by using the transmission unit 362.
  • the output result of this time can be used for grasping the short-term or long-term stress state of the observation target.
  • the graph creating unit 312 sets the pulsation interval RR (n) in the range used for calculating the standard deviation from the pulsation interval table of FIG. 11 and the next.
  • the data with the beat interval RR (n + 1) of the above is acquired, and Poincare plot diagrams as shown in FIGS. 14 to 17 are created.
  • the result output unit 313 displays the created Poincare plot diagram on its own or external output device such as a display.
  • the graph creation unit 312 may create and output a Poincare plot diagram after axis conversion by using the result of step S108.
  • FIG. 14 is a Poincare plot diagram of the dog in the excited state according to the present embodiment.
  • FIG. 15 is a Poincare plot diagram of a dog according to the present embodiment in a state in which respiration is stable in a normal state.
  • FIG. 16 is a Poincare plot diagram of the dog according to the present embodiment in a normal state.
  • FIG. 17 is a Poincare plot diagram of the dog according to the present embodiment in a resting state.
  • the heart rate increases (the pulsation interval becomes shorter), the fluctuation of the pulsation interval becomes smaller, and the plot is made. It becomes a state where the points of are gathered in a certain place.
  • the heart rate is not as low as in the resting state (the spread of the plot points is not as large as in the resting state), but at the center of the distribution of the plot points.
  • the analysis unit 311 calculates the degree of variation in the Poincare plot, that is, the standard deviation of the pulsation interval, as a numerical value indicating the autonomic nerve balance. ⁇ Another form of autonomic balance values>
  • the product of these two standard deviations may be calculated as a numerical value indicating the autonomic nerve balance.
  • FIG. 18 is a flowchart showing a processing procedure of the biometric information processing system 1 according to the present embodiment. Since steps S100 to S108 are the same as those in FIG. 10, the description is not repeated here.
  • the CPU 310 as the analysis unit 311 calculates the standard deviation for each axis after the axis conversion (step S110).
  • the analysis unit 311 may specify the axis having the maximum variance and calculate the standard deviation of the axis and the axis perpendicular to the axis.
  • the analysis unit 311 calculates the product of these two standard deviations, the square root of the product, and the like as numerical values indicating the autonomic nerve balance (step S112).
  • the result output unit 313 stores the standard deviation, the product of the standard deviations, the square root of the product, the beat interval table, etc. in the data storage unit 322, or uses the transmission unit 362 to use the transmission unit 362, etc., via the Internet or the like. And send it to the server 100.
  • the output result of this time can be used for grasping the short-term or long-term stress state of the observation target.
  • the analysis unit 311 calculates the product of the standard deviations of the two axes, the square root of the product, etc., but may also calculate the product of the standard deviations of the three or more axes, the root of the product, and the like. ..
  • the CPU 310 performs the above calculation every predetermined period, for example, several minutes, and stores the calculation result in the database of the memory 320 for creating a diagnostic graph described later.
  • the CPU 510 executes a memory program to execute a peak detection unit 5122, a peak storage unit 5123, and a peak comparison unit 5124.
  • a peak detection unit 5122 In addition to the beat timing determination unit 5125 and the beat timing storage unit 5126, it is preferable to realize the beat interval analysis unit 5127, the elapsed time setting unit 5128, and the like. Since the peak detection unit 5122, the peak storage unit 5123, the peak comparison unit 5124, the pulsation timing determination unit 5125, and the pulsation timing storage unit 5126 are the same as those in the above embodiment, the description is not repeated here.
  • the pulsation interval analysis unit 5127 frequently determines whether or not the peaks determined to be pulsation timing frequently exist at intervals of a little longer than the predetermined time in step S106.
  • the elapsed time setting unit 5128 resets the predetermined period to a longer value when the peaks determined to be the pulsation timing frequently exist at intervals slightly longer than the predetermined time.
  • the pulsation timing acquisition unit 512 executes the following processing every time a predetermined period, for example, 1 minute, elapses.
  • a predetermined period for example, 1 minute
  • the pulsation timing acquisition unit 512 generates a pulsation interval histogram showing the time between pulsation timings for a predetermined period (step S132).
  • FIGS. 23 and 24 are histogram of the beat interval in a predetermined period.
  • FIG. 24 is a graph showing the waveform of the electrocardiographic signal and the pulsation interval from the pulsation timing immediately before each peak determined to be the pulsation timing in a predetermined period.
  • the pulsation in the range of 200 ms to 400 ms is shown in FIGS. 25 and 26. Many intervals are detected.
  • the mode of the pulsation interval may be in the range of 200 ms to 400 ms, or the frequency of the pulsation interval of 200 ms to 400 ms may be within the top three.
  • the pulsation timing acquisition unit 512 determines whether or not a large peak of the histogram exists in the range of 200 ms to 400 ms (step S134). If YES in step S134, the pulsation timing acquisition unit 512 lengthens the predetermined time in step S106 by 20 ms up to 400 ms (step S136).
  • the pulsation timing acquisition unit 512 may shorten the second predetermined time by 20 ms with 200 ms as the lower limit.
  • the signal acquisition device 500 acquires the beat timing and the beat interval based on the electrocardiographic signals from the electrodes 401, 402, 403, and the communication terminal 300 beats.
  • the information for judging the state of the living thing or the information of the judgment result of the state of the living thing was calculated and output from the interval.
  • the roles of all or part of one of those devices may be shared by another device or by a plurality of devices. On the contrary, one device may play the role of all or a part of the plurality of devices, or another device may play the role.
  • the server 100 may play the role of the communication terminal 300.
  • the communication terminal 300 transmits necessary information such as the pulsation timing and the pulsation interval from the signal acquisition device 500 to the server 100 via a router, a carrier network, the Internet, or the like.
  • the server 100 calculates the information for determining the state of the living thing or the information indicating the determination result of the state of the living thing, transmits the information to the communication terminal 300, and the communication terminal 300 displays the information of the final result. And output to the speaker.
  • the communication terminal 300 or the server 100 may realize the function of the signal acquisition device 500.
  • the communication terminal 300 or the server 100 may receive an electrocardiographic signal from the signal acquisition device 500, specify the pulsation timing from a plurality of peaks, or calculate the pulsation interval.
  • the beat is detected by detecting the peak of the component derived from the electrocardiogram based on the measurement results of the plurality of electrodes, the measuring unit for measuring the potential difference between the plurality of electrodes, and the measuring unit.
  • a biometric information processing system including a control unit for determining dynamic timing is provided. The control unit executes the following processing with respect to the first peak having an amplitude larger than a predetermined threshold value among the plurality of peaks detected by the detection unit. (1) Within the first period including the time of the first peak, three or more peaks including the first peak are specified in descending order of amplitude in the positive direction and the negative direction. (2) Within the second period shorter than the first period including the time of the first peak, one of three or more peaks in the direction opposite to the first peak is formed before and after the first peak. When it is determined that they exist one by one, the first peak is specified as the beat timing.
  • the control unit sets the peak in the direction opposite to the first peak among the three or more peaks to the first peak outside the third period included in the second period.
  • the first peak is determined as the beat timing.
  • the third period is 0 ms or more and 30 ms or less.
  • the second period is 30 ms or more and 100 ms or less.
  • the first period is 400 ms or more and 800 ms or less.
  • control unit selects the target first peak from the peak at the time when a predetermined time or more has passed from the previous beat timing and the determined first peak.
  • control unit selects the target first peak from the peak at the time when a time of 200 ms or more and 10000 ms or less has elapsed from the first peak determined as the previous beat timing.
  • control unit Preferably, the control unit generates a histogram regarding the interval between beat timings, and changes a predetermined time based on the histogram.
  • control unit changes the predetermined time to 200 ms or more and 400 ms or less.
  • the electrocardiographic waveform of the living body is acquired, and among the plurality of peaks detected when the electrocardiographic waveform is acquired, the following processing is performed on the first peak having an amplitude larger than a predetermined threshold value.
  • a biometric information acquisition method to be performed is provided. (1) Within the first period including the time of the first peak, three or more peaks including the first peak are specified in descending order of amplitude in the positive direction and the negative direction. (2) Within the second period shorter than the first period including the time of the first peak, one of three or more peaks in the direction opposite to the first peak is formed before and after the first peak. When it is determined that they exist one by one, the first peak is specified as the beat timing.
  • the living body is an animal having a respiratory arrhythmia.
  • Biometric information processing system 100 Server 110: CPU 300: Communication terminal 310: CPU 311: Analysis unit 312: Graph creation unit 313: Result output unit 320: Memory 321: Beat interval acquisition unit 322: Data storage unit 330: Display 360: Communication interface 361: Reception unit 362: Transmission unit 401: Electrode 402: Electrode 403: Electrode 500: Signal acquisition device 510: CPU 511: Electrocardiographic preprocessing unit 512: Pulsation timing acquisition unit 512B: Pulsation timing acquisition unit 514: Peak comparison unit 520: Storage unit 560: Transmission unit 5122: Peak detection unit 5123: Peak storage unit 5124: Peak comparison unit 5125 : Pulsation timing determination unit 5126: Pulsation timing storage unit 5127: Pulsation interval analysis unit 5128: Elapsed time setting unit

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Animal Behavior & Ethology (AREA)
  • Environmental Sciences (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Physics & Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Husbandry (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Physiology (AREA)
  • Zoology (AREA)
  • Signal Processing (AREA)
  • Pulmonology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

複数の電極間の電位差を計測する計測部の計測結果に基づいて、心電に由来する成分のピークを検出することによって、拍動タイミングを決定する制御部（５１０）を備え、前記制御部（５１０）は、検出された複数のピークのうち、所定の閾値よりも振幅の大きい第１ピークに関して次の処理を実行する。（１）第１のピークの時点を含む第１の期間内に、第１ピークを含み、かつ、正の方向および負の方向に振幅の大きい順にそれぞれ３つ以上のピークを特定する。（２）第１ピークの時点を含む第１の期間よりも短い第２の期間内に、３つ以上ピークのうちの第１ピークと反対方向のピークが、第１ピークの前後にそれぞれ１つずつ存在すると判断したときに、第１ピークを拍動タイミングと特定する。

Description

生体情報処理システムおよび生体情報取得方法
　本国際出願は、２０１９年４月２２日に日本国特許庁に出願された日本国特許出願第２０１９－０８１３７１号に基づく優先権を主張するものであり、日本国特許出願第２０１９－０８１３７１号の全内容を参照により本国際出願に援用する。
　以下の開示は、生物の精神的状態または肉体的状態を取得するための技術に関する。
　従来から、生物の精神的または肉体的な状態を取得するための技術が知られている。例えば、特開２０１６－５４８７７号公報（特許文献１）には、生体情報測定装置が開示されている。特許文献１によると、生体の心拍間隔を測定可能で生体に着用可能な生体情報測定装置であって、生体に接触する複数の電極と、心電図信号を作成する信号処理手段と、心電図信号から該心電図信号中の一のＲ波及び該一のＲ波と隣り合う他のＲ波の間隔若しくは一のＳ波及び該一のＳ波と隣り合う他のＳ波の間隔から心拍間隔を測定する心拍間隔測定手段と、三軸加速度測定手段と、温度測定手段とを備え、前記心拍間隔測定手段によって得られた心拍間隔と前記三軸加速測定手段によって得られた三軸加速度と前記温度測定手段によって得られた温度とを、同時に無線で送信する無線送信手段を備える。
　また、特開２０１７－２９６２８号公報（特許文献２）には、心拍検出方法および心拍検出装置が開示されている。特許文献２によると、心拍検出装置は、心電図波形のサンプリングデータの時間差分値を算出する時間差分値算出部と、時間差分値が閾値を超えているかどうかを判定する時間差分値判定部と、想定される心拍時刻の手前の第一の時間間隔の範囲と時間差分値のピークを含む第二の時間間隔の範囲と時間差分値のピークの後の第三の時間間隔の範囲を求める時間判定部と、第一、第二、第三の時間間隔の範囲における時間差分値の最小値Ｍｉｎ１，Ｍｉｎ２，Ｍｉｎ３を保持する最小値保持部と、最小値Ｍｉｎ１，Ｍｉｎ２，Ｍｉｎ３の関係が心拍時刻確定条件を満たすとき、時間差分値が閾値を超えた時刻または最小値Ｍｉｎ２が得られた時刻を心拍時刻とする心拍時刻決定部を備える。
特開２０１６－５４８７７号公報 特開２０１７－２９６２８号公報
　本開示の目的は、より正確に拍動タイミングを特定するための技術を提供することにある。
　本開示の一態様に従うと、複数の電極と、複数の電極間の電位差を計測する計測部と、計測部の計測結果に基づいて、心電に由来する成分のピークを検出することによって、拍動タイミングを決定する制御部とを備える生体情報処理システムが提供される。制御部は、検出された複数のピークのうち、所定の閾値よりも振幅の大きい第１ピークに関して次の処理を実行する。（１）第１のピークの時点を含む第１の期間内に、第１ピークを含み、かつ、正の方向および負の方向に振幅の大きい順にそれぞれ３つ以上のピークを特定する。（２）第１ピークの時点を含む第１の期間よりも短い第２の期間内に、３つ以上ピークのうちの第１ピークと反対方向のピークが、第１ピークの前後にそれぞれ１つずつ存在すると判断したときに、第１ピークを拍動タイミングと特定する。
　以上のように、本開示によれば、より正確に拍動タイミングを特定するための技術が提供される。
第１の実施の形態にかかる生体情報処理システム１の全体構成を示す図である。 第１の実施の形態にかかる生体情報処理システム１の機能構成を示す図である。 第１の実施の形態にかかる信号取得装置５００の処理手順を示すフローチャートである。 第１の実施の形態にかかる心電データと拍動間隔を示す図面である。 第１の実施の形態にかかる拍動タイミングの送信データを示す図である。 第１の実施の形態にかかる心電データと拍動タイミングを示す図面である。 第１の実施の形態にかかるタイミング取得部５１２の機能構成を示す図である。 第１の実施の形態にかかる拍動タイミング決定処理の処理手順を示すフローチャートである。 第１の実施の形態にかかる対象となるピークの近傍の心電データの拡大図である。 第１の実施の形態にかかる生体情報処理システム１の第１の自律神経バランスを算出するための処理手順を示すフローチャートである。 第１の実施の形態にかかる拍動間隔テーブルを示す図である。 第１の実施の形態にかかる拍動間隔Ｒ－Ｒ（ｎ）とその次の拍動間隔Ｒ－Ｒ（ｎ＋１）との拍動間隔テーブルからＹ＝Ｘ方向とそれに垂直な方向の軸への変換を示すイメージ図である。 第１の実施の形態にかかる犬の精神的状態または肉体的状態毎の、Ｙ＝Ｘ軸に関する標準偏差と、Ｙ＝Ｘと垂直な軸に関する標準偏差との目安を示す表である。 第１の実施の形態にかかる犬の興奮状態におけるポアンカレプロット図である。 第１の実施の形態にかかる犬の通常状態で呼吸が安定している状態におけるポアンカレプロット図である。 第１の実施の形態にかかる犬の通常状態におけるポアンカレプロット図である。 第１の実施の形態にかかる犬の安静状態におけるポアンカレプロット図である。 第１の実施の形態にかかる生体情報処理システム１の第２の自律神経バランスを算出するための処理手順を示すフローチャートである。 第１の実施の形態にかかる犬の精神的または肉体的状態毎の、Ｙ＝Ｘ軸に関する標準偏差と、Ｙ＝Ｘと垂直な軸に関する標準偏差と、第２の自律神経バランスとしての標準偏差の積と、標準偏差の比との目安を示す表である。 第２の実施の形態にかかる対象となるピークの近傍の心電データの拡大図である。 第２の実施の形態にかかるタイミング取得部５１２Ｂの機能構成を示す図である。 第２の実施の形態にかかる所定時間変更処理の処理手順を示すフローチャートである。 第２の実施の形態にかかる、好ましい所定時間に基づいて特定された拍動タイミングのヒストグラムを示す図面である。 第２の実施の形態にかかる、心電信号と、好ましい所定時間に基づいて特定された拍動タイミングとを示す図面である。 第２の実施の形態にかかる、短すぎる所定時間に基づいて特定された拍動タイミングのヒストグラムを示す図面である。 第２の実施の形態にかかる、心電信号と、短すぎる所定時間に基づいて特定された拍動タイミングとを示す図面である。
　以下、図面を参照しつつ、本開示の実施の形態について説明する。以下の説明では、同一の部品には同一の符号を付してある。それらの名称および機能も同じである。したがって、それらについての詳細な説明は繰り返さない。
　＜第１の実施の形態＞
　＜生体情報処理システムの全体構成＞
　まず、図１および図２を参照して、本実施の形態にかかる生体情報処理システム１の全体構成について説明する。図１は、本実施の形態にかかる生体情報処理システム１の全体構成を示す図である。図２は、本実施の形態にかかる生体情報処理システム１の機能構成を示す図である。本実施の形態にかかる生体情報処理システム１は、呼吸性の不整脈を有する動物にも適用することができる。なお、以下では、呼吸性の不整脈を有する動物を代表して犬の状態を判断する場合について説明する。
　本実施の形態にかかる生体情報処理システム１は、主に、犬の胸部に取り付けられる心電取得用の電極４０１，４０２，４０３と、心電信号を処理するための状態取得装置としての信号取得装置５００と、信号取得装置５００と通信可能な状態処理装置としての通信端末３００とを含む。
　心電取得用の電極４０１，４０２，４０３は、胸部等において、心臓部を挟むような位置に取り付けることが望ましく、例えば、両前足（または、前足と後ろ足）の肉球部など毛の生えていない場所であってもよい。また、毛を刈った状態であるか、ゲルなどが付着した電極、あるいは、突起状の構造を持ち、毛があっても皮膚と接触する構成であることが望ましい。あるいは、毛がある状態で、非接触で容量性材料を介して心電を誘導する形態が望ましい。それにより、犬等の表皮が毛に覆われた生物であっても心電を取得することが可能となる。本実施の形態においては、３個の電極４０１，４０２，４０３を使用する構成としているが、電極は、２個以上であればよく、さらに、多くの電極を使用する構成としてもよい。
　＜信号取得装置５００の構成＞
　次に、図２および図３を参照して、信号取得装置５００の構成と処理とについて説明する。なお、図３は、本実施の形態にかかる信号取得装置５００が実行する全体処理を示すフローチャートである。信号取得装置５００は、心電前処理部５１１と拍動タイミング取得部５１２と送信部５６０を含む。
　心電前処理部５１１は、フィルタや増幅器を含む。心電前処理部５１１は、電極４０１，４０２，４０３から送られている信号を、図４に示すような心電信号データに変換して、拍動タイミング取得部５１２に受け渡す（ステップＳ００２）。
　より詳細には、心電前処理部５１１には、ハイパスフィルタ、ローパスフィルタなどのフィルタ装置、オペアンプなどから構成される増幅装置、心電のアナログ信号をデジタル信号に変換するＡ／Ｄ変換装置等が含まれる（ステップＳ００４）。尚、フィルタ装置、増幅装置などは、ソフトウェアにより実装される形態であってもよい。また、Ａ／Ｄ変換装置においては、拍動間隔のゆらぎ量の差異が判別できる周期と精度でのサンプリングを行うことが望ましい。すなわち、Ａ／Ｄ変換の周波数が、２５Ｈｚ以上の周波数で取得することが望ましい。例えば、本実施の形態においては、１００Ｈｚでの心電信号のサンプリングを行っている。サンプリングの周波数を高めることにより、拍動間隔の揺らぎ量を正確に把握することが可能となる。
　拍動タイミング取得部５１２は、たとえば、ＣＰＵ（Central Processing Unit）５１０やメモリを含む制御部によって実現される。より詳細には、ＣＰＵ５１０が、メモリに格納されているプログラムを実行することによって、後述する判断方法により、心電のピーク信号のうち、Ｒ波と判断できるピークを特定する（ステップＳ００６）。拍動タイミング取得部５１２は、Ｒ波と判断できるピークに対して拍動タイミングとしての時刻情報を付与したり、拍動タイミングから拍動タイミングまでの時間を拍動間隔として特定したりする（ステップＳ００８）。なお、拍動間隔の算出方法としては、自己相関関数を用いた周期の導出や、矩形波相関トリガを用いる方法、特徴点を複数検出する方法などによるピーク検出で行ってもよい。
　本実施の形態においては、拍動タイミング取得部５１２は、連続して入力される心電信号に対して連続して拍動タイミングの特定を実行する。拍動タイミング取得部５１２は、拍動毎に、拍動タイミングを示す時刻を、記憶部５２０に記憶したり、送信部５６０を介して通信端末３００に送信したりする（ステップＳ０１０）。例えば、ＣＰＵ５１０は、図５に示すように、所定の時間帯毎に、当該所定の時間帯の基準時刻を示すタイムスタンプと、当該所定の時間帯に含まれる拍動タイミング毎の、所定の時間帯の基準時刻からの詳細な経過時間とを記憶部５２０に記憶させたり、送信部５６０を介して通信端末３００に送信したりする。本実施の形態においては、所定の時間帯の長さは１秒間であって、所定の時間帯の基準時刻からの詳細な経過時間はｍｓｅｃで示される。
　なお、記憶部５２０は、一例として、ＳＤカードやＵＳＢメモリなどによって実現される。送信部５６０は、アンテナやコネクタなどを含む通信インターフェイスによって実現される。
　＜拍動タイミングの特定方法＞
　次に、信号取得装置５００における拍動タイミングを特定するための構成について説明する。信号取得装置５００の拍動タイミング取得部５１２は、図５に示すようなノイズが少ない心電データや図６に示すようなノイズの多い心電データが得られた場合に、複数のピークから拍動タイミングに対応するピークを特定したり抽出したりする。なお、図６に示す心電信号ａに関しては、犬などの動物への取り付けが上手くいかなかったり、動物が動いてしまったりすることなどによって、拍動タイミングが特定しにくいエリアｃが生じている場合を示す。このような場合などにおいても、本実施の形態にかかる拍動タイミング取得部５１２は、下記の構成や処理によって、正確な拍動タイミングｂを特定できる可能性が高い。
　図７は、本実施の形態にかかる拍動タイミング取得部５１２の機能構成を示すブロック図である。図７および図２に示すように、本実施の形態にかかる拍動タイミング取得部５１２は、ＣＰＵ５１０やメモリを含む制御部によって実現される。より詳細には、ＣＰＵ５１０が、メモリに格納されているプログラムを実行することによって、ピーク検出部５１２２やピーク比較部５１２４や拍動タイミング決定部５１２５などを実現し、記憶部５２０がピーク記憶部５１２３や拍動タイミング記憶部５１２６などを実現する。
　ピーク検出部５１２２は、心電データから対象となるピークを選択したり、当該ピークを基準とした所定の期間内の複数のピークを抽出したりする。ピーク記憶部５１２３は、それらのピークを特定するための時間情報などを記憶する。ピーク比較部５１４は、抽出された複数のピークを、所定のルールに従って比較する。拍動タイミング決定部５１２５は、比較結果に基づいて、所定の条件を満たすピークを拍動タイミングとして決定する。拍動タイミング記憶部５１２６は、拍動タイミングの時間情報を蓄積する。
　次に、図８と図９とを参照して、拍動タイミング取得部５１２の処理について説明する。なお、図８は、本実施の形態にかかる拍動タイミング特定処理の処理手順を示すフローチャートである。本実施の形態にかかる拍動タイミング取得部５１２は、心電データを取得すると以下の処理を実行することによって、マイナス方向またはプラス方向に心電ピークが大きく出ている波形に対して当該ピークを検出して拍動のタイミングを認識する。
　なお、フィルタを通過した後のＡＤ変換後の心電信号は、図９に示すように歪を生じた形状になっていることが多い。図９は、本実施の形態にかかる、拍動タイミングであるか否かの判断の対象となるピークの近辺に関する心電データの拡大図である。
　以下では、マイナス方向に突出するピークに関して、拍動タイミングであるか否かを判断する場合の処理について説明するが、拍動タイミング取得部５１２は、同様の処理をプラス方向に突出するピークに関しても実行することによって、両方に突出したピークに関して拍動タイミングを検知できるように構成されている。なお、拍動タイミングがマイナス方向とプラス方向のどちらのピークに該当するかは、犬に取り付けられた電極の向きに依存する。
　拍動タイミング取得部５１２は、ある時点における心電波形に関して、マイナス方向の閾値Ｖ１、例えば－０．２Ｖ、以下まで突出しており、その時点の前後の微分値の積がマイナスであって、直前のピークとの時間間隔が１０ｍｓよりも大きいか否かを判断する（ステップＳ１０２）。ステップＳ１０２にてＹＥＳである場合、拍動タイミング取得部５１２は、当該時点にピークが存在すると判断する。以下、拍動タイミング取得部５１２は、当該ピークＰ１に関して、拍動タイミングの可能性が高いか否かを判断する。
　拍動タイミング取得部５１２は、当該ピークＰ１の時点を含む第１の期間Ｔ１において、マイナス方向の閾値Ｖ１以下まで突出し、振幅の大きな順に４つのピークＰｍ１、Ｐｍ２、Ｐｍ３、Ｐｍ４と、プラス方向の閾値Ｖ２、例えば０．２Ｖ、以上まで突出し、振幅の大きな順に４つのピークＰｐ１、Ｐｐ２、Ｐｐ３、Ｐｐ４とを特定して、メモリに記憶する（ステップＳ１０４）。ここで、第１の期間Ｔ１は、例えば、当該ピークの時点の前後に２００ｍｓ～４００ｍｓ、すなわち４００ｍｓ～８００ｍｓとすることが好ましい。第１の期間Ｔ１が４００ｍｓ未満の場合、測定の対象が犬の場合、ノイズをピークとして検出し、心拍数が大きく算出される恐れがある。一方で、第１の期間Ｔ１が８００ｍｓを超える場合、隣接する拍動を検出期間に含んでしまい、当該拍動と隣接する拍動のそれぞれの拍動タイミングを検出するべきところを、どちらか一方のみを検出する恐れがある。
　拍動タイミング取得部５１２は、今回対象としているピークＰｍ１が、前回の拍動タイミングと判断されたピークから所定時間、例えば２００ｍｓなど、以上１００００ｍｓ以下経過しているものであるか否かを判断する（ステップＳ１０６）。
　ステップＳ１０６にてＹＥＳである場合、拍動タイミング取得部５１２は、当該ピークの信号値Ｐｍ１が、第１の期間Ｔ１において、最小点であるか否かを判断する（ステップＳ１０８）。
　ステップＳ１０８にてＹＥＳである場合、拍動タイミング取得部５１２は、第１の期間Ｔ１よりも短く、当該ピークＰｍ１を含む第２の期間Ｔ２の間において、４つのピークＰｐ１、Ｐｐ２、Ｐｐ３、Ｐｐ４のうち当該ピークＰｍ１の前後にそれぞれ１つずつピークが存在するか否かを判断する。例えば、図９では、第２の期間Ｔ２の間において、当該ピークＰｍ１の前のピークＰｐ１が存在し、当該ピークＰｍ１の後ろにピークＰｐ２が存在する。さらに好ましくは、第２の期間Ｔ２において、プラス側の最大のピークと２番から４番の少なくともいずれかのピークが存在する、または、プラス側の２番目のピークと３番から４番の少なくともいずれかのピークが存在する、か否かを判断する（ステップＳ１１０）。このとき、第２の期間Ｔ２は、当該ピークＰｍ１の前後５０ｍｓ以内、好ましくは４５ｍｓ以内であることが好ましい。
　より好ましくは、第２の期間Ｔ２よりも短く、当該ピークＰｍ１を含む第３の期間Ｔ３の期間には、４つのピークＰｐ１、Ｐｐ２、Ｐｐ３、Ｐｐ４が存在するか否かも判断してもよい。例えば、第３の期間Ｔ３は、当該ピークＰｍ１の前後２０ｍｓ以内であることが好ましい。
　より詳細には、今回対象としている当該ピークＰｍ１よりも前か後の一方において、当該ピークＰｍ１の時間差が２０ｍｓから４５ｍｓの範囲に、プラス側の最大のピークＰｐ１が存在し、今回対象としているピークＰｍ１よりも前か後の他方において、時間差が２０ｍｓから４５ｍｓの範囲に、Ｐｐ２～Ｐｐ４の少なくともいずれかのピークが存在するか否かを判断する。当てはまらない場合は、今回対象としているピークＰｐ１よりも前か後の一方において、時間差が２０ｍｓから４５ｍｓの範囲に、プラス側の２番目のピークＰｐ２が存在し、今回対象としているピークＰｍ１よりも前か後の他方において、時間差が２０ｍｓから４５ｍｓの範囲に、プラス側のＰｐ３かＰｐ４の少なくともいずれかのピークが存在するか否かを判断する。
　ステップＳ１１０にてＹＥＳである場合、拍動タイミング取得部５１２は、当該ピークＰｍ１の振幅が、プラス側の１番または２番のピークＰｐ１、Ｐｐ２よりも大きいか否かを判断する（ステップＳ１１２）。
　ＣＰＵ１１０は、ステップＳ１１２にてＹＥＳである場合、拍動タイミング取得部５１２は、当該ピークのタイミングを拍動タイミングと判断して、当該タイミングを示す時間情報を記憶部５２０に記憶する（ステップＳ１１４）。
　本実施の形態においては、１秒毎のタイムスタンプからの差分（ｍｓ）を拍動時刻として記録する。信号取得装置５００は、所定の期間分の拍動タイミングを示す情報を通信端末３００に送信する。
　以上の通り、本実施の形態にかかる拍動タイミング取得部５１２は、ステップＳ１０６、ステップＳ１０８、ステップＳ１１０、ステップＳ１１２などの判断を行うことによって、正確なＲ波を特定するものである。すなわち、他の種類のピークや体動、電磁波、物理的な衝撃によるノイズを、拍動タイミングであると特定してしまう可能性を低減することができるものである。
　＜通信端末３００の構成＞
　次に、図２および図１０を参照して、通信端末３００の構成と処理について説明する。なお、図１０は、本実施の形態にかかる通信端末３００が実行する処理を示すフローチャートである。通信端末３００は、受信部３６１、拍動間隔取得部３２１、解析部３１１と、グラフ作成部３１２と、結果出力部３１３と、ディスプレイ３３０と、データ記憶部３２２と、送信部３６２とを含む。
　まず、受信部３６１と送信部３６２は、例えば、アンテナやコネクタなどを含む通信インターフェイス３６０によって実現される。受信部３６１は、信号取得装置５００からの拍動タイミングを示すデータを受信する（ステップＳ１００）。
　拍動間隔取得部３２１は各種のメモリ３２０などによって構成され、信号取得装置５００から受信したデータを格納する。本実施の形態においては、ＣＰＵ３１０が、通信インターフェイス３６０を介して受信した拍動タイミングに基づいて拍動と拍動の間の時間を計算する（ステップＳ１０２）。ＣＰＵ３１０は、拍動間隔を拍動間隔テーブル（図１１参照。）として逐次メモリ３２０に蓄積していく（ステップＳ１０４）。本実施の形態においては、拍動間隔は、例えば、図１１に示すようにｍｓｅｃ（ミリセック）の単位で計算される。ただし、これらのデータは、通信端末３００のメモリ３２０に記憶されてもよいし、通信端末３００からアクセス可能な他の装置に記憶されてもよい。なお、図１１において拍動時刻１及び２に９が記入されている時刻の行は、もともと拍動が検出されていない時間帯示すものである。これは、拍動が無い時間帯を仮に９の数値を入力して示すものであり、その他の数値や記号であってもよい。
　なお、本実施の形態においては、ＣＰＵ３１０は、何らかの原因によって信号取得装置５００からの拍動タイミングを示すデータが欠損した場合、後述する通り、拍動間隔が計算できない期間に関しては各種のプロットを行わず、当該期間を経過したところから各種のプロットを再開する。より詳細には、欠損データの取り扱い方法として後述する。
　解析部３１１と、グラフ作成部３１２と、結果出力部３１３とは、例えばＣＰＵ３１０がメモリ３２０のプログラムを実行することによって実現される。解析部３１１は、一定時間単位、例えば、１分、１０分、１時間など、状態を判定するために必要な時間単位で、拍動間隔取得部３２１から拍動間隔データを読み出して、拍動間隔Ｒ－Ｒ（ｎ）とその次の拍動間隔Ｒ－Ｒ（ｎ＋１）との拍動間隔テーブルを作成する（ステップＳ１０６）。
　解析部３１１は、図１２に示すように、拍動間隔Ｒ－Ｒ（ｎ）とその次の拍動間隔Ｒ－Ｒ（ｎ＋１）との拍動間隔テーブルからＹ＝Ｘ方向とそれに垂直な方向の軸への変換を行う（ステップＳ１０８）。
　解析部３１１は、軸の変換を行った後のそれぞれの軸を構成する数値列に関する標準偏差を算出する（ステップＳ１１０）。なお、解析部３１１は、Ｙ＝Ｘ軸に関する標準偏差だけを算出してもよいし、Ｙ＝Ｘと垂直な軸に関する標準偏差だけを算出してもよいし、両方を算出してもよい。図１３は、犬の精神状態または肉体的状態毎の、Ｙ＝Ｘ軸に関する標準偏差と、Ｙ＝Ｘと垂直な軸に関する標準偏差との目安を示す表である。
　なお、解析部３１１は、主成分分析などの方法により分散が最大になる軸を特定し、当該軸と当該軸に垂直な軸に関する標準偏差を算出してもよい。さらには、解析部３１１は、軸変換を行わずに、Ｘ軸とＹ軸に関する標準偏差を算出するものであってもよい。分散の大きい方向がＸ軸方向とＹ軸方向である場合には、軸変換を行わなくとも、Ｘ軸とＹ軸の標準偏差を算出することで、ポアンカレプロットした拍動間隔のばらつき状態を評価できる。この場合、軸変換を行う必要が無いために、計算量を低減することができる。
　結果出力部３１３は、例えば、自身の、あるいは外部の、ディスプレイ３３０やスピーカなどの出力装置に、標準偏差を表示させたり、音声メッセージを出力させたりする（ステップＳ１１４）。より詳細には、結果出力部３１３は、Ｙ＝Ｘ軸に関する標準偏差だけを出力させてもよいし、Ｙ＝Ｘと垂直な軸に関する標準偏差だけを出力させてもよいし、両方を出力させてもよいし、大きい方だけを出力させてもよいし、小さい方だけを出力させてもよい。
　標準偏差を計算することにより、拍動間隔Ｒ－Ｒ（ｎ）とその次の拍動間隔Ｒ－Ｒ（ｎ＋１）とをそれぞれ軸としてポアンカレプロットした拍動間隔のばらつき状態が評価できる。
　なお、本実施の形態にかかる生体情報処理システム１は、通信端末３００が通信可能なサーバを含む形態であってもよい。その場合、結果出力部３１３としてのＣＰＵ３１０は、標準偏差や関係テーブルなどデータ記憶部３２２に蓄積したり、送信部３６２を利用することによって、インターネットなどを介してサーバに送信したりする。これによって、今回の出力結果を観察対象の短期または長期のストレス状態の把握などに利用することができる。
　本実施の形態においては、ステップＳ１０８とは別に、同時にグラフ作成部３１２は、図１１の拍動間隔テーブルから、標準偏差の計算に使用した範囲の拍動間隔Ｒ－Ｒ（ｎ）とその次の拍動間隔Ｒ－Ｒ（ｎ＋１）とのデータを取得して、図１４～図１７に示すようなポアンカレプロット図を作成する。
　そして、結果出力部３１３は、作成されたポアンカレプロット図を、自身の、または外部の、ディスプレイなどの出力装置に表示させる。なお、グラフ作成部３１２は、ステップＳ１０８の結果を利用して、軸変換後のポアンカレプロット図を作成して出力してもよい。
　ここで、ポアンカレプロット図に関して説明する。図１４は、本実施の形態にかかる犬の興奮状態におけるポアンカレプロット図である。図１５は、本実施の形態にかかる犬の通常状態で呼吸が安定している状態におけるポアンカレプロット図である。図１６は、本実施の形態にかかる犬の通常状態におけるポアンカレプロット図である。図１７は、本実施の形態にかかる犬の安静状態におけるポアンカレプロット図である。
　まず、例えば犬などの呼吸性の不整脈を有する生物の場合、図１４のような興奮状態においては、心拍数が上昇し（拍動間隔は短くなる）、拍動間隔の揺らぎは小さくなり、プロットの点が一定の場所に集まるような状態になる。
　そして、図１５のような呼吸が安定している通常の状態においては、心拍数が安静状態ほどは少なくない（プロットの点の広がりが安静状態ほど大きくない）が、プロット点の分布の中心にプロットが少ない（穴の空白）領域が存在する。このような形状になるのは、犬の心拍が呼吸の影響を大きく受けるため、拍動変動が周期的に変化することが原因と考えられる（呼吸性不整脈）。そのため、リラックスした緩やかな拍動ではないが、呼吸が安定して行われているため、空白の存在する状態になると考えられる。
　そして、図１６のような通常状態においては、拍動に揺らぎがみられ、ばらつきは大きくなる（プロット点が広がる）が、プロット点が散乱している状態となる。
　そして、図１７の安静状態においては、犬がリラックスしているために拍動の間隔が大きくなり、さらに呼吸性不整脈の影響を大きく受けるために、プロット点の広がりが大きくなると共に、円形や四角形に近い形状や、三角形に近い形状となる。そのいずれの形状においても、安静状態ではポアンカレプロットのプロット点の分布の中心部に空白部分が見られる形状となる。
　このように、本実施の形態においては、算出結果に基づいて間接的に、ポアンカレプロットのプロット点の分布の広がりの大きさや形状、中心部にプロットが多くみられるか少なくみられるかを予想することができ、その結果、生物の精神的状態または肉体的状態を予想することができる。そして、上述した通り、解析部３１１は、自律神経バランスを示す数値として、ポアンカレプロットのバラツキ具合すなわち拍動間隔の標準偏差を算出するものである。
　＜自律神経バランスの数値に関する別の形態＞
　上記の実施の形態においては、通信端末３００が、ポアンカレプロットのＹ＝Ｘの軸に沿った標準偏差またはＹ＝Ｘと垂直な軸に沿った標準偏差を出力するものであった。しかしながら、自律神経バランスを示す数値として、それら２つの標準偏差の積を算出してもよい。以下では、図１８を参照して、本実施の形態にかかる生体情報処理システム１の処理手順について説明する。
　図１８は、本実施の形態にかかる生体情報処理システム１の処理手順を示すフローチャートである。ステップＳ１００～ステップＳ１０８は、図１０のものと同様であるため、ここでは説明を繰り返さない。
　解析部３１１としてのＣＰＵ３１０は、軸の変換を行った後のそれぞれの軸に関する標準偏差を算出する（ステップＳ１１０）。なお、解析部３１１は、分散が最大になる軸を特定し、当該軸と当該軸に垂直な軸に関する標準偏差を算出してもよい。
　そして、解析部３１１は、自律神経バランスを示す数値として、それらの２つの標準偏差の積や積の平方根などを計算する（ステップＳ１１２）。
　結果出力部３１３は、例えば、通信端末３００の、または外部の、ディスプレイやスピーカなどの出力装置に、標準偏差の積や積の平方根などを表示させたり、音声メッセージを出力させたりする（ステップＳ１１４）。より詳細には、結果出力部３１３は、Ｙ＝Ｘ軸に関する標準偏差と、Ｙ＝－Ｘの軸に関する標準偏差と、両者の積や積の平方根などとを出力させてもよい。
　図１９は、犬の精神状態または肉体的状態毎の、Ｙ＝Ｘ軸に関する標準偏差と、Ｙ＝Ｘと垂直な軸に関する標準偏差と、自律神経バランスを示す数値としての標準偏差の積や積の平方根などと、標準偏差の比との目安を示す表である。
　標準偏差の積を計算することにより、拍動間隔Ｒ－Ｒ（ｎ）とその次の拍動間隔Ｒ－Ｒ（ｎ＋１）とをそれぞれ軸としてポアンカレプロットした拍動間隔の分布の広がりの大きさや形状、一様に分散している、中心に空白がある等のばらつき状態が評価できる。また、縦横比が同じで大きさのみ変化している状態や分布の広がり面積が同じで中心部のばらつき状態が異なる場合などに有効にばらつき状態を評価できる。
　この場合も、結果出力部３１３は、標準偏差や標準偏差の積や積の平方根や拍動間隔テーブルなどをデータ記憶部３２２に蓄積したり、送信部３６２を利用することによって、インターネットなどを介してサーバ１００に送信したりする。これによって、今回の出力結果を観察対象の短期または長期のストレス状態の把握などに利用することができる。
　解析部３１１は、２つの軸の標準偏差の積や積の平方根などを計算するものであるが、３つ以上の軸の標準偏差の積やその累乗根などを計算するものであってもよい。
　ＣＰＵ３１０は、所定の期間、例えば数分間、毎に上記の計算を行い、当該計算結果を後述する診断グラフ作成のためにメモリ３２０のデータベースに蓄積していく。
　＜第２の実施の形態＞
　心電信号の取得や拍動間隔の取得においては、様々な要因により、図２０に示すように、拍動タイミング以外の時点で、拍動タイミングのピークと同様に大きいピークが検出されてしまう場合がある。このような場合は、両者を拍動タイミングであると認識してしまう可能性が高い。このような不具合が生じる可能性をさらに低減するために、ステップＳ１０６の所定時間やステップＳ１０４の第１の所定期間が、自動的に調整されることが好ましい。
　例えば、本実施の形態にかかる拍動タイミング取得部５１２Ｂに関しては、図２１に示すように、ＣＰＵ５１０は、メモリのプログラムを実行することによって、ピーク検出部５１２２やピーク記憶部５１２３やピーク比較部５１２４や拍動タイミング決定部５１２５や拍動タイミング記憶部５１２６などに加えて、拍動間隔解析部５１２７や、経過時間設定部５１２８などを実現することが好ましい。なお、ピーク検出部５１２２やピーク記憶部５１２３やピーク比較部５１２４や拍動タイミング決定部５１２５や拍動タイミング記憶部５１２６は、上記の実施の形態と同様であるため、ここでは説明を繰り返さない。
　拍動間隔解析部５１２７は、拍動タイミングと判断したピーク同士がステップＳ１０６の所定時間よりも少しだけ長い程度の間隔をあけて存在することが頻繁に生じているか否かを判断する。経過時間設定部５１２８は、拍動タイミングと判断したピーク同士が所定時間よりも少しだけ長い程度の間隔に存在することが頻繁に生じている場合に、所定期間を長めに設定し直す。
　具体的には、拍動タイミング取得部５１２は、所定の期間、例えば１分など、が経過するたびに以下のような処理を実行する。図２２を参照して、拍動タイミング取得部５１２は、所定の期間分の、拍動タイミング間の時間を示す拍動間隔のヒストグラムを生成する（ステップＳ１３２）。
　例えば、拍動タイミングの近傍において、拍動タイミングのピークと同様に大きいピークが検出されてしまうことが少ない状況に関しては、図２３や図２４に示すように、２００ｍｓ～４００ｍｓの範囲の拍動間隔が検出されることが少ない。なお、図２３は、所定の期間における拍動間隔のヒストグラムである。図２４は、所定の期間における、心電信号の波形と、拍動タイミングと判断されたピーク毎の直前の拍動タイミングからの拍動間隔と、を示すグラフである。
　逆に、拍動タイミングの近傍において、拍動タイミングのピークと同様に大きいピークが検出されてしまうことが多い状況に関しては、図２５や図２６に示すように、２００ｍｓ～４００ｍｓの範囲の拍動間隔が多く検出される。例えば、拍動間隔の最頻値が２００ｍｓ～４００ｍｓの範囲にあったり、２００ｍｓ～４００ｍｓの拍動間隔の頻度が上位３番目以内であったりする。
　このようにして、拍動タイミング取得部５１２は、２００ｍｓ～４００ｍｓの範囲にヒストグラムの大きなピークが存在するか否かを判断する（ステップＳ１３４）。ステップＳ１３４にてＹＥＳである場合、拍動タイミング取得部５１２は、４００ｍｓを上限として、ステップＳ１０６の所定時間を２０ｍｓだけ長くする（ステップＳ１３６）。
　なお、ステップＳ１３４にてＹＥＳである場合、拍動タイミング取得部５１２は、２００ｍｓを下限として、第２の所定時間を２０ｍｓだけ短くしてもよい。
　＜第３の実施の形態＞
　上記の実施の形態にかかる生体情報処理システム１は、電極４０１，４０２，４０３からの心電信号に基づいて信号取得装置５００が拍動タイミングや拍動間隔を取得し、通信端末３００が拍動間隔から生物の状態を判断するための情報または生物の状態の判定結果の情報を算出して出力するものであった。しかしながら、それらの１つの装置の全部または一部の役割が、別の装置によって担われてもよいし、複数の装置によって分担されてもよい。逆に、それら複数の装置の全部または一部の役割を、１つの装置が担ってもよいし、別の装置が担ってもよい。
　例えば、通信端末３００の役割をサーバ１００が担ってもよい。例えば、通信端末３００が信号取得装置５００からの拍動タイミングや拍動間隔などの必要な情報をルータやキャリア網やインターネットなどを介してサーバ１００に送信する。そして、サーバ１００が生物の状態を判断するための情報または生物の状態の判定結果を示す情報を算出し、当該情報を通信端末３００に送信し、通信端末３００が最終的な結果の情報をディスプレイやスピーカに出力する。
　あるいは、信号取得装置５００の機能を通信端末３００やサーバ１００が実現しても良い。例えば、通信端末３００やサーバ１００が、信号取得装置５００から心電信号を受信して、複数のピークから拍動タイミングを特定したり、拍動間隔を計算したりしてもよい。
　＜まとめ＞
　上記の実施の形態においては、複数の電極と、複数の電極間の電位差を計測する計測部と、計測部の計測結果に基づいて、心電に由来する成分のピークを検出することによって、拍動タイミングを決定する制御部とを備える生体情報処理システムが提供される。制御部は、検出部で検出された複数のピークのうち、所定の閾値よりも振幅の大きい第１ピークに関して次の処理を実行する。（１）第１のピークの時点を含む第１の期間内に、第１ピークを含み、かつ、正の方向および負の方向に振幅の大きい順にそれぞれ３つ以上のピークを特定する。（２）第１ピークの時点を含む第１の期間よりも短い第２の期間内に、３つ以上ピークのうちの第１ピークと反対方向のピークが、第１ピークの前後にそれぞれ１つずつ存在すると判断したときに、第１ピークを拍動タイミングと特定する。
　好ましくは、制御部は、（２）の処理において、３つ以上ピークのうちの第１ピークと反対方向のピークが、第２の期間に含まれる第３の期間外に、第１のピークの前後にそれぞれ１つずつ存在することを特定した時に、第１ピークを拍動タイミングと決定する。第３の期間は０ｍｓ以上３０ｍｓ以下である。第２の期間は３０ｍｓ以上１００ｍｓ以下である。
　好ましくは、第１の期間は、４００ｍｓ以上８００ｍｓ以下である。
　好ましくは、制御部は、１つ前の拍動タイミングと決定された第１のピークから所定の時間以上が経過した時点のピークから対象とする第１のピークを選択する。
　好ましくは、制御部は、１つ前の拍動タイミングと決定された第１のピークから２００ｍｓ以上１００００ｍｓ以下の時間が経過した時点のピークから対象とする第１ピークを選択する。
　好ましくは、制御部は、拍動タイミングの間隔に関するヒストグラムを生成し、当該ヒストグラムに基づいて、所定の時間を変更する。
　好ましくは、制御部は、所定の時間を２００ｍｓ以上４００ｍｓ以下に変更する。
　上記の実施の形態においては、生体の心電波形を取得し、心電波形を取得するときに検出された複数のピークのうち、所定の閾値よりも振幅の大きい第１ピークに関して次の処理を実行する生体情報取得方法が提供される。
（１）第１のピークの時点を含む第１の期間内に、第１ピークを含み、かつ、正の方向および負の方向に振幅の大きい順にそれぞれ３つ以上のピークを特定し、
（２）第１ピークの時点を含む第１の期間よりも短い第２の期間内に、３つ以上ピークのうちの第１ピークと反対方向のピークが、第１ピークの前後にそれぞれ１つずつ存在すると判断したときに、第１ピークを拍動タイミングと特定する。
　好ましくは、生体は、呼吸性の不整脈を有する動物である。
　今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本開示の範囲は、上記した説明ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれることが意図される。
１　　　　：生体情報処理システム
１００　　：サーバ
１１０　　：ＣＰＵ
３００　　：通信端末
３１０　　：ＣＰＵ
３１１　　：解析部
３１２　　：グラフ作成部
３１３　　：結果出力部
３２０　　：メモリ
３２１　　：拍動間隔取得部
３２２　　：データ記憶部
３３０　　：ディスプレイ
３６０　　：通信インターフェイス
３６１　　：受信部
３６２　　：送信部
４０１　　：電極
４０２　　：電極
４０３　　：電極
５００　　：信号取得装置
５１０　　：ＣＰＵ
５１１　　：心電前処理部
５１２　　：拍動タイミング取得部
５１２Ｂ　：拍動タイミング取得部
５１４　　：ピーク比較部
５２０　　：記憶部
５６０　　：送信部
５１２２　：ピーク検出部
５１２３　：ピーク記憶部
５１２４　：ピーク比較部
５１２５　：拍動タイミング決定部
５１２６　：拍動タイミング記憶部
５１２７　：拍動間隔解析部
５１２８　：経過時間設定部

Claims (9)

  1. 　複数の電極と、
    　前記複数の電極間の電位差を計測する計測部と、
    　前記計測部の計測結果に基づいて、心電に由来する成分のピークを検出することによって、拍動タイミングを決定する制御部とを備え、
    　前記制御部は、
    　検出された複数のピークのうち、所定の閾値よりも振幅の大きい第１ピークに関して次の処理を実行する、
    （１）前記第１のピークの時点を含む第１の期間内に、前記第１ピークを含み、かつ、正の方向および負の方向に振幅の大きい順にそれぞれ３つ以上のピークを特定し、
    （２）前記第１ピークの時点を含む前記第１の期間よりも短い第２の期間内に、前記３つ以上ピークのうちの前記第１ピークと反対方向のピークが、前記第１ピークの前後にそれぞれ１つずつ存在すると判断したときに、前記第１ピークを前記拍動タイミングと特定する、生体情報処理システム。
  2. 　前記制御部は、前記（２）の処理において、
    　前記３つ以上ピークのうちの前記第１ピークと反対方向のピークが、前記第２の期間に含まれる第３の期間外に、前記第１のピークの前後にそれぞれ１つずつ存在することを特定した時に、前記第１ピークを前記拍動タイミングと決定し、
    　前記第３の期間は０ｍｓ以上３０ｍｓ以下であって、
    　前記第２の期間は３０ｍｓ以上１００ｍｓ以下である、請求項１に記載の生体情報処理システム。
  3. 　前記第１の期間は、４００ｍｓ以上８００ｍｓ以下である、請求項１または２に記載の生体情報処理システム。
  4. 　前記制御部は、１つ前の拍動タイミングと決定された第１のピークから所定の時間以上が経過した時点のピークから対象とする前記第１のピークを選択する、請求項１から３のいずれか１項に記載の生体情報処理システム。
  5. 　前記制御部は、１つ前の拍動タイミングと決定された第１のピークから２００ｍｓ以上１００００ｍｓ以下の時間が経過した時点のピークから対象とする前記第１ピークを選択する、請求項４に記載の生体情報処理システム。
  6. 　前記制御部は、拍動タイミングの間隔に関するヒストグラムを生成し、当該ヒストグラムに基づいて、前記所定の時間を変更する、請求項４または５に記載の生体情報処理システム。
  7. 　前記制御部は、前記所定の時間を２００ｍｓ以上４００ｍｓ以下に変更する、請求項６に記載の生体情報システム。
  8. 　生体の心電波形を取得し、
    　前記心電波形を取得するときに検出された複数のピークのうち、所定の閾値よりも振幅の大きい第１ピークに関して次の処理を実行する
    （１）前記第１のピークの時点を含む第１の期間内に、前記第１ピークを含み、かつ、正の方向および負の方向に振幅の大きい順にそれぞれ３つ以上のピークを特定し、
    （２）前記第１ピークの時点を含む前記第１の期間よりも短い第２の期間内に、前記３つ以上ピークのうちの前記第１ピークと反対方向のピークが、前記第１ピークの前後にそれぞれ１つずつ存在すると判断したときに、前記第１ピークを前記拍動タイミングと特定する、生体情報取得方法。
  9. 　前記生体は、呼吸性の不整脈を有する動物である、ことを特徴とする請求項８に記載の生体情報取得方法。
PCT/JP2020/016887 2019-04-22 2020-04-17 生体情報処理システムおよび生体情報取得方法 Ceased WO2020218195A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019081371A JP2022096000A (ja) 2019-04-22 2019-04-22 生体情報処理システムおよび生体情報取得方法
JP2019-081371 2019-04-22

Publications (1)

Publication Number Publication Date
WO2020218195A1 true WO2020218195A1 (ja) 2020-10-29

Family

ID=72942032

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/016887 Ceased WO2020218195A1 (ja) 2019-04-22 2020-04-17 生体情報処理システムおよび生体情報取得方法

Country Status (2)

Country Link
JP (1) JP2022096000A (ja)
WO (1) WO2020218195A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001346771A (ja) * 2000-06-09 2001-12-18 Canon Inc Ｒ波認識方法及びｒ−ｒ間隔測定方法及び心拍数測定方法及びｒ−ｒ間隔測定装置及び心拍数測定装置
JP2003000561A (ja) * 2001-06-18 2003-01-07 Canon Inc Ｒ波認識方法及びｒ−ｒ間隔測定方法及び心拍数測定方法及びｒ−ｒ間隔測定装置及び心拍数測定装置
JP2009112625A (ja) * 2007-11-08 2009-05-28 Almedio Inc 心拍測定装置
WO2016208152A1 (ja) * 2015-06-25 2016-12-29 パナソニックＩｐマネジメント株式会社 心電位検出装置および心電位検出方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001346771A (ja) * 2000-06-09 2001-12-18 Canon Inc Ｒ波認識方法及びｒ−ｒ間隔測定方法及び心拍数測定方法及びｒ−ｒ間隔測定装置及び心拍数測定装置
JP2003000561A (ja) * 2001-06-18 2003-01-07 Canon Inc Ｒ波認識方法及びｒ−ｒ間隔測定方法及び心拍数測定方法及びｒ−ｒ間隔測定装置及び心拍数測定装置
JP2009112625A (ja) * 2007-11-08 2009-05-28 Almedio Inc 心拍測定装置
WO2016208152A1 (ja) * 2015-06-25 2016-12-29 パナソニックＩｐマネジメント株式会社 心電位検出装置および心電位検出方法

Also Published As

Publication number Publication date
JP2022096000A (ja) 2022-06-29

Similar Documents

Publication Publication Date Title
JP7303854B2 (ja) 評価装置
US20180256048A1 (en) Pulse wave transit time measurement device and living body state estimation device
EP3795069B1 (en) Apparatus and method for estimating bio-information
JP6647415B2 (ja) 状態取得コンピュータ、状態の取得方法、および情報処理システム
KR20190120684A (ko) 생체신호 측정 상태 모니터링 장치 및 방법과, 생체정보 측정 장치 및 방법
WO2016072940A1 (en) Multi-channel ballistocardiography with cepstrum smoothing and quality-based dynamic channel selection
JP6766710B2 (ja) 血圧測定装置、方法及びプログラム
KR20160107390A (ko) 심전도 신호 측정장치
US20180055373A1 (en) Monitoring device to identify candidates for autonomic neuromodulation therapy
JP2021129693A (ja) 心臓疾患重症度判定システム
WO2018055996A1 (ja) コンピュータ、呼吸数の取得方法、および情報処理システム
KR20160039506A (ko) 일상생활 중 심부 체온 일주기 리듬 추정 방법
KR20020018541A (ko) 신체 및 감정 상태 인식 장치 및 방법
WO2020218195A1 (ja) 生体情報処理システムおよび生体情報取得方法
US20230165507A1 (en) Biological signal processing device, watching system, and watching method
CN108712883A (zh) 生物体信息测定装置、个人识别装置、个人识别方法以及个人识别程序
JP2021137332A (ja) 生体情報処理装置、生体情報処理システムおよび生体情報処理方法
JPWO2019198691A1 (ja) 情報処理装置、およびウェアラブル端末
JPWO2019193998A1 (ja) 状態取得装置、ネットワークシステム、および状態処理装置
EP3930570B1 (en) Measurement unit and monitoring system for monitoring indicator of parkinson's disease in person
JP2021137250A (ja) 生体情報処理装置、生体情報処理システムおよび生体情報処理方法
JP2022074177A (ja) 生体情報測定用ハーネスおよび生体情報処理システム
WO2022091195A1 (ja) Ｒｒｉ計測装置、ｒｒｉ計測方法およびｒｒｉ計測プログラム
KR101519960B1 (ko) 12채널 심전도 측정을 위한 일체형 다중전극 및 이를 이용하여 12채널 심전도 검출방법
US20200323478A1 (en) Information processing device, state acquisition program, server, and information processing method

Legal Events

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

Ref document number: 20796343

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20796343

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: JP