WO2008035611A1 - Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données - Google Patents

Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données Download PDF

Info

Publication number
WO2008035611A1
WO2008035611A1 PCT/JP2007/067832 JP2007067832W WO2008035611A1 WO 2008035611 A1 WO2008035611 A1 WO 2008035611A1 JP 2007067832 W JP2007067832 W JP 2007067832W WO 2008035611 A1 WO2008035611 A1 WO 2008035611A1
Authority
WO
WIPO (PCT)
Prior art keywords
measurement signal
data
basic
signal
data processing
Prior art date
Application number
PCT/JP2007/067832
Other languages
English (en)
French (fr)
Inventor
Atsushi Okumoto
Original Assignee
Mitsubishi Chemical Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Mitsubishi Chemical Corporation filed Critical Mitsubishi Chemical Corporation
Publication of WO2008035611A1 publication Critical patent/WO2008035611A1/ja

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/02Instruments for indicating weather conditions by measuring two or more variables, e.g. humidity, pressure, temperature, cloud cover or wind speed
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20ICT 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

Definitions

  • Data processing apparatus data processing method, and data processing program
  • the present invention relates to a data processing device, a data processing method, and a data processing program for analyzing various states of an object such as a state of a living body such as a machine, animals, plants, and microorganisms, or a natural phenomenon such as weather or an earthquake. It is about.
  • parameters that vary according to the state of an object are measured as time-series measurement signals.
  • Patent Document 1 collects a subject's pulse wave waveform as a measurement signal and extracts the order contained in the measurement signal (that is, the deterministic structure that governs the fluctuation of the measurement signal).
  • a technique for diagnosing the condition of a subject by performing an operation to perform the above is disclosed. In this technology, by calculating the chaos attractor and Lyapunov exponent of the measurement signal, it is possible to logically extract the order inherent in the measurement signal, and to enable objective diagnosis to the subject. .
  • Patent Document 2 discloses a technique for diagnosing the health condition of a subject by analyzing data of a biological signal detected from the subject.
  • Data analysis methods in this technology include chaos analysis, detrend fluctuation analysis (DFA), frequency conversion, wavelet analysis, and multifractal analysis.
  • This preprocessing includes, for example, noise reduction processing (noise removal processing) for removing noise mixed in the measurement signal, filter processing (filtering) for extracting a predetermined frequency component from the measurement signal, orthogonal transformation processing, Fourier transform Conversion processing, envelope processing, etc.
  • Patent Document 1 Japanese Patent Publication No. 6-9546
  • Patent Document 2 Japanese Patent Laid-Open No. 2001-299766
  • the information that is removed in the process of pre-processing may include important information that affects the subsequent arithmetic processing and the arithmetic results in addition to simple noise and unnecessary information.
  • the characteristic omission (information omission) associated with pre-processing may reduce the calculation accuracy in subsequent calculation processing, resulting in inaccurate calculation results.
  • chaos with irregular behavior exists in normal measurement signals that measure animals and plants and natural phenomena.
  • Chaos is an unpredictable complex behavior caused by the nonlinear deterministic system inherent in the measurement object, and it is a noise (irregular behavior that does not depend on a deterministic system, Although it is conceptually different from unnecessary information other than the information that becomes (), it is extremely difficult to distinguish between these strictly!
  • the present invention has been devised in view of such problems, and an object of the present invention is to provide a data processing device, a data processing method, and a program that can improve data processing speed and data processing accuracy with a simple configuration. And Also provided is a data processing device, a data processing method, and a data processing program that can obtain an accurate calculation result in a short time with a simple configuration in calculation processing for extracting a non-linear structure in a measurement signal. With the goal.
  • a data processing device of the present invention includes measurement signal detection means (signal detection device) that detects, as a measurement signal, a parameter that varies according to the state of the object.
  • Signal processing means for generating a basic measurement signal by performing signal processing as preprocessing for grasping fluctuations of the parameters with respect to the measurement signal detected by the measurement signal detection means (basic measurement A basic data extraction unit (basic data extraction unit) that extracts, as basic data, a measurement signal that characterizes the variation of the parameter based on the basic measurement signal generated by the signal processing unit.
  • An extraction range setting unit for setting a predetermined region defined by the basic data extracted by the basic data extraction unit as the extraction range of the measurement signal, and the measurement signal detection hand Feature data for extracting, as feature data, a measurement signal characterizing the state of the object from the measurement signals included in the extraction range set by the extraction range setting means from among the measurement signals detected in step It is characterized by having an extraction means (feature data extraction unit)
  • the preprocessing in the signal processing means is easy to find variations in the parameters.
  • the data processing device of the present invention according to claim 2 is characterized in that, in the configuration according to claim 1, the measurement signal detecting means detects an animal vital sign as the measurement signal.
  • the vital sign here means a physical quantity as a vital sign detected from the body of an animal (including a person).
  • body movements and respiratory rate associated with exercise such as walking, heart rate, body temperature, skin surface temperature, skin potential, pulse wave (pulse rate), brain wave, blood flow, saliva, and other body fluid components, respiratory air and blood Oxygen saturation, blood glucose level, electrocardiogram, electrical conductivity, body weight (pressure on the seating surface), number and cycle of blinks, sweating, and other electromagnetic wave intensity and chemical substance concentration .
  • the signal processing means performs the signal processing on the measurement signal using a linear analysis technique. It is characterized by giving.
  • the data processing device of the present invention is characterized in that, in the configuration described in any one of claims;! To 3, the signal processing means is the measurement signal detected by the measurement signal detection means. It has a filter processing means (basic measurement signal generation unit) for filtering a predetermined frequency component set in advance from the signal.
  • the filter processing means basic measurement signal generation unit
  • the data processing device of the present invention according to claim 5 is the configuration described in any one of claims;! To 4, wherein the measurement signal detection means is substantially periodic as the measurement signal.
  • the signal processing unit generates a wave obtained by smoothing the measurement signal as the basic measurement signal, and the basic data extraction unit detects a fluctuation peak in the basic measurement signal.
  • the detection time is extracted as the basic data.
  • “to smooth the measurement signal” means to make the waveform of the measurement signal data a smooth wave (waveform).
  • the extraction range setting means sets a time near the detection time of the peak as the extraction range. Is special.
  • the feature data extraction means extracts the detection time of the fluctuation peak in the measurement signal included in the extraction range as the feature data. It is characterized by doing.
  • the data processing device of the present invention is based on the feature data extracted by the feature data extraction means in the configuration described in any one of claims;! To 7.
  • the method further comprises arithmetic processing means (second data processing unit) for extracting a nonlinear structure in the measurement signal and analyzing the state of the object.
  • a data processing method wherein a measurement signal detecting step for detecting a parameter that varies according to the state of an object as a measurement signal, and the measurement signal detected in the measurement signal detection step A signal processing step for generating a basic measurement signal by performing signal processing as a pre-processing for grasping fluctuations of the parameter, and a processing based on the basic measurement signal generated in the signal processing step.
  • the data processing method of the present invention described in claim 10 detects an animal vital sign as the measurement signal in the measurement signal detection step. Is special.
  • the data processing method of the present invention according to claim 11 is characterized in that, in addition to the configuration according to claim 9 or 10, the signal processing step performs the signal processing on the measurement signal using a linear analysis technique.
  • the data processing method of the present invention according to claim 12 is the data processing method according to any one of claims 9 to 11;
  • a predetermined frequency component set in advance is filtered from the measurement signal detected in the measurement signal detection step.
  • the data processing method of the present invention described in claim 13 is, in addition to the configuration described in any one of claims 9 to 12, fluctuating substantially periodically as the measurement signal in the measurement signal processing step.
  • the signal processing step a wave obtained by smoothing the measurement signal is generated as the basic measurement signal, and in the basic data extraction step, the detection time of the fluctuation peak in the basic measurement signal is calculated.
  • the special feature is to extract it as basic data.
  • the data processing method of the present invention described in claim 14 sets the time near the peak detection time as the extraction range in the extraction range setting step.
  • the data processing method of the present invention according to claim 15 is characterized in that, in addition to the configuration according to claim 14, in the feature data extraction step, the detection time of the fluctuation peak in the measurement signal included in the extraction range is represented by the feature data. It is characterized by extracting as.
  • the data processing method of the present invention is based on the feature data extracted in the feature data extraction step in addition to any one of claims 9 to 15; It is characterized by further comprising an arithmetic processing step for extracting a nonlinear structure from the measurement signal and analyzing the state of the object.
  • the data processing program of the present invention is a data processing program for causing a computer to function as a signal processing means, basic data extraction means, extraction range setting means, and feature data extraction means.
  • the means performs signal processing as a preprocessing for grasping the fluctuation of the parameter detected according to the state of the object detected as the measurement signal, and generates a basic measurement signal.
  • the basic data extraction means extracts the measurement signal characterizing the variation of the parameter as basic data based on the basic measurement signal generated by the signal processing means, and the extraction range setting means A predetermined area defined by the basic data extracted by the basic data extracting means is set as an extraction range of the measurement signal, and the feature data extracting means The Among the measurement signals detected by the measurement signal detection means, the measurement signal characterizing the state of the object is extracted as feature data from the measurement signals included in the extraction range set by the extraction range setting means. It is characterized by doing.
  • the measurement signal subjected to signal processing as preprocessing in setting the extraction range of the measurement signal In the extraction of more specific feature data, the measurement signal that characterizes the state of the object can be extracted extremely accurately because it is extracted from the measurement signal included in the extraction range.
  • signal processing as preprocessing is easy.
  • processing can be completed in a short time with a simple configuration.
  • the data processing device and the data processing method of the present invention since the peak of the basic measurement signal generated by smoothing the measurement signal is extracted, the basic data can be easily obtained. Can be detected.
  • the peak of the basic measurement signal is used. Measurement signals included in a distant range can be excluded from the feature data extraction targets. In other words, the correlation with the measurement signal that characterizes the state of the object is strong V, and it is possible to easily extract the information S.
  • the measurement signal included in the measurement signal and characterizing the state of the object can be accurately extracted.
  • the data processing apparatus and the data processing method of the present invention (claims 8 and 16), the nonlinear structure in the measurement signal can be extracted based on the accurately extracted feature data, and the data analysis is highly reliable. It can be performed.
  • FIG. 1 is a block diagram showing an overall configuration of a data processing apparatus according to an embodiment of the present invention.
  • FIG. 2 is a graph for explaining the contents of data processing in this data processing device, (a) is a time-series graph of measurement signals related to data processing in the basic measurement signal generator and basic data extractor, (b) ) Is a time series graph of measurement signals related to data processing in the feature data extraction unit.
  • FIG. 3 is a flowchart showing control details in the data processing apparatus.
  • FIG. 4 is a schematic diagram showing a configuration example of the data processing apparatus using a computer. Explanation of symbols
  • this data processing device is a device that detects a parameter (for example, acceleration) corresponding to a human walking state as a detection signal, performs data processing on the detected signal, and outputs the detected signal.
  • a parameter for example, acceleration
  • the data processing apparatus includes a signal detection device (measurement signal detection means) 1, a first data processing unit 9, and a second data processing unit 10.
  • the first data processing unit 9 performs arithmetic processing to make it easy to grasp the characteristics of the signals detected by the signal detection device 1, while the second data processing unit 10 performs substantial data processing.
  • the walking state is analyzed as follows.
  • the first data processing unit 9 and the second data processing unit 10 are functional parts that are arithmetically processed inside the computer, and each function is configured as an individual program.
  • the first data processing section 9 in this embodiment functions as a signal processing means, basic data extraction means, extraction range setting means, and feature data extraction means.
  • the second data processing unit 10 functions as an arithmetic processing means.
  • Fig. 4 shows a configuration example of the data processing apparatus using a computer.
  • the computer 12 includes the signal detection device 1 described above, a storage device (ROM, RAM, etc.) 13, a central processing unit (CPU) 14, a monitor as an output interface 15, a keyboard 16 and a mouse 17 as an input interface. It is configured.
  • the first data processing unit 9 and the second data processing unit 10 according to the data processing apparatus are stored in the storage device 13 as programs.
  • the signal detection device 1 detects various conditions (variation factors) related to the state of life such as machines, animals, plants and microorganisms, natural phenomena such as weather and earthquakes, and various object states. It is a sensor. This parameter includes not only information directly detected from the sensor but also information obtained by processing sensor detection information by calculation or the like and using the corresponding parameter value as an estimated value.
  • an acceleration sensor for detecting an acceleration signal that is a calculation target of the data processing apparatus is applied as the signal detection apparatus 1 and is mounted on the body of a target person.
  • This acceleration sensor may be one to three axes depending on the measurement object and purpose, but it detects acceleration acting in three directions: vertical direction, horizontal front-rear direction and horizontal left-right direction during walking. It is preferable to use a three-axis acceleration sensor.
  • a triaxial acceleration sensor is used, and the detected information of the acceleration in the vertical direction and the detection time information detected here are processed as a measurement signal S as shown in FIG. Input to part 9.
  • the first data processing unit 9 is a functional part that applies the processing specified in claim 1 of the present application to the measurement signal S. As shown in FIG. 1, the measurement signal storage unit 11, the basic measurement signal generation unit 2, and the basic data Data extraction unit 3, extraction range setting unit 4 and feature data extraction unit 5.
  • the processing performed here is to convert the signal by applying mathematical and electrical processing to the measurement signal S consisting of an optical signal, audio signal, electromagnetic signal, etc., in order to make it easier to grasp its characteristics. (In other words, signal processing).
  • the processing can be classified into analog signal processing and digital signal processing depending on the type of signal to be processed.
  • the first data processing unit 9 performs processing included in the category of digital signal processing.
  • the measurement signal storage unit 11 is a functional part that stores the measurement signal S input from the signal detection device 1.
  • the measurement signal S stored here is stored in the measurement signal storage unit 11 as shown in FIG. Are divided into two systems and input to the basic measurement signal generator 2 and the feature data extractor 5, respectively. That is, raw information that is not processed at all is input to each of the basic measurement signal generation unit 2 and the feature data extraction unit 5.
  • the measurement signal S input to the measurement signal storage unit 11 may be a series of time-series data detected during a predetermined time in the signal detection device 1, or may be signal detection. It may be individual measurement data detected at any time by the apparatus 1.
  • the measurement signal S as time series data is input to the basic measurement signal generation unit 2 and the feature data extraction unit 5 as it is.
  • each measurement signal S is stored in the measurement signal storage unit 11 until the total number of measurement signals S from the signal detection device 1 becomes equal to or greater than a predetermined number, and the number of data When all of these are fully prepared, the entire measurement signal S is input to the basic measurement signal generation unit 2 and the feature data extraction unit 5.
  • the latter case will be described.
  • the basic measurement signal generation unit 2 is a functional part that performs processing for extracting a necessary signal component from the time-series measurement signal S input from the signal detection device 1.
  • a digital filter process for detecting the peak of acceleration change during walking is performed, and three types of filter processes, a low-pass filter, a phase filter, and a high-pass filter, are applied to the time series data of acceleration. It's like! /
  • the specific content of the processing performed by the basic measurement signal generator 2 differs depending on the type and content of the measurement signal data and the purpose of the analysis. However, in general, linear analysis methods such as filter processing, Fourier transform processing, and envelope processing are preferably performed alone or in combination.
  • the low-pass filter is a filter that reduces (or removes) high-frequency vibration components from an input signal including various frequency components.
  • This low-pass filter suppresses noise components with higher frequencies than the frequency of acceleration fluctuations associated with a person walking (that is, the reciprocal of the time interval from when one foot touches the ground until the other foot touches the ground) It is like that.
  • the vibrations at fine time intervals seen in the measurement signal S indicated by the solid line in Fig. 2 (a) are high-frequency noise components.
  • the phase filter is a filter that delays the phase change of the input signal.
  • the phase of the input signal is delayed by a quarter period of the acceleration fluctuation.
  • the high-pass filter is a filter that reduces (removes!
  • the drift component of the signal detected by the signal detection device 1 is suppressed.
  • the time series data of the acceleration fluctuations accompanying the walking of the person from the measurement signal S is clarified.
  • a signal as indicated by the broken line is obtained. It is taken out.
  • the characteristic of the walking time interval is extracted from the raw information detected by the signal detection device 1 as a wave having a corresponding period.
  • the time series signal thus extracted is hereinafter referred to as a basic measurement signal S.
  • this basic measurement signal S is a signal extracted through a phase filter.
  • the basic data extraction unit 3 includes the basic measurement signal S generated by the basic measurement signal generation unit 2.
  • information characterizing acceleration fluctuations is extracted.
  • information that characterizes acceleration fluctuations is estimated from the basic measurement signal S obtained by filtering.
  • the peak position of the input measurement signal S is extracted.
  • the basic data extraction unit 3 changes the signal value from negative to positive before and after the 0 points of the basic measurement signal S.
  • the measurement signal S corresponding to the zero point is extracted. More specifically, the measurement signal S at the same time as the time 0 is extracted as basic data D in Fig. 2 (a).
  • a method of adjusting the peak position of the measurement signal S and the basic data D by adjusting the content of the filter processing is the following method.
  • the basic data D is simply extracted based on the basic measurement signal S.
  • the configuration includes an extraction range setting unit 4 and a feature data extraction unit 5 described below.
  • the extraction range setting unit 4 determines the time near the time when the basic data D is detected as the extraction range A.
  • the neighborhood time is before or after the time when the basic data D is detected.
  • the extraction range A has a predetermined time t before the time when the basic data D is detected and a predetermined time after the time when the basic data D is detected.
  • the extraction range A set here will be explained subsequently. It is input to the feature data extraction unit 5.
  • the predetermined times t and t are times that can be arbitrarily set. For example, basic
  • It may be determined as the ratio of the measurement signal S to the wavelength, or it may be a preset value.
  • each predetermined time t and t force are estimated to be generated by the above-described filter processing.
  • the predetermined times t and t may be relatively short.
  • 1 t should be set relatively long.
  • the feature data extraction unit 5 extracts, as feature data D, a signal that characterizes the walking state of a person, that is, an actual peak position, from the measurement signal S included in the extraction range A set by the extraction range setting unit 4. To do.
  • the calculation process here is shown in Fig. 2 (b).
  • the extraction range setting unit 4 and the feature data extraction unit 5 are based on the basic data D.
  • the extraction range setting unit 4 sets the extraction range A on the assumption that there is information indicating the future characteristics. Furthermore, the feature data extraction unit 5 reduces only the measurement signal S included in the extraction range A to be calculated, thereby reducing calculation labor and calculation time, and extracting the feature data D from the measurement signal S. The accuracy is ensured.
  • the second data processing unit 10 is a functional part for performing substantial data processing on the data processed in the first data processing unit 9, and as shown in FIG.
  • the analyzing unit 7 and the determining unit 8 are provided.
  • This second data processing unit 10 The feature of peak interval time of the acceleration data of either the left or right foot during walking is analyzed.
  • a non-linear analysis method for observing fluctuations in peak interval time is used.
  • the term "fluctuation” here refers to a slight waveform shift (spatial, temporal change or movement that is partially irregular) that is observed when a certain wave changes every moment. pointing.
  • vital signs such as respiratory rate, heart rate, and brain waves
  • the peak interval and period of those waves are It is known to exhibit complex fluctuations that are not constant.
  • a number of methods have been proposed to analyze the structure that seems to dominate the behavior among the complex variations that appear irregular.
  • the second data processing unit 10 analyzes the nonlinear structure that would exist behind the fluctuation by observing the degree of fluctuation of the peak interval time by using such a method. is there.
  • Specific analysis techniques include known analysis techniques such as spectrum analysis (FFT analysis), fractal analysis (multi-fractal analysis, detrend fluctuation analysis, etc.), chaos analysis, and wavelet analysis. Then, detrend fluctuation analysis, which is one of the fractal analysis methods, is used.
  • FFT analysis spectrum analysis
  • fractal analysis multi-fractal analysis
  • detrend fluctuation analysis etc.
  • chaos analysis chaos analysis
  • wavelet analysis wavelet analysis
  • the detrend analysis method is a statistical analysis method that evaluates the complexity of the wave to be analyzed using a value called the scaling index.
  • the acceleration data of either the left or right foot is aligned in the data alignment unit 6, the scaling index of the acceleration data is calculated in the analysis unit 7, and the evaluation is performed in the determination unit 8. It has become.
  • the data alignment unit 6 arranges the feature data D extracted by the feature data extraction unit 5 for convenience of calculation in the analysis unit 7.
  • feature data D extracted by the feature data extraction unit 5 for convenience of calculation in the analysis unit 7.
  • the feature data D is alternately sorted and divided in this data alignment unit 6 in order to observe the walking state in more detail.
  • the extracted feature data D are arranged in the order of their detection times
  • the odd-numbered feature data D group is
  • the interval (the time interval at which one foot contacts the ground) can be grasped.
  • the analysis unit 7 individually specifies a scaling instruction for each data group of the even-numbered and odd-numbered feature data D groups input from the data alignment unit 6 of the second data processing unit 10.
  • one feature data group D is assigned n based on the detection time.
  • the least square error (variance) F between the interval time and the trend is calculated, and the slope ⁇ of the logarithmic plot of each of the division number n and the variance F is calculated as a scaling index.
  • the trend means the trend of data transition in each section. For example, the data in each section is approximated to a straight line.
  • the magnitude of the self-similarity of the degree of variation in the actual walking interval time when the observation interval is changed is calculated as the scaling index ⁇ .
  • the determination unit 8 determines the walking state based on the scaling index ⁇ calculated by the analysis unit 7.
  • l / f fluctuation means fluctuation among the above-mentioned fluctuations such that the magnitude of the fluctuation component (power spectrum) is 1 / f with respect to the frequency f. ing.
  • 1 / f fluctuations whose power spectrum is inversely proportional to the frequency f can be observed.
  • 1 / f fluctuations whose power spectrum is inversely proportional to the frequency f can be observed.
  • 1 / f fluctuation is used as an index for judging the health condition of the observation target.
  • step A10 for explaining the control contents in the data processing apparatus using the flowchart shown in FIG. 3, the acceleration detection information and the detection time information are detected as the measurement signal S by the acceleration sensor as the signal detection apparatus 1.
  • the detected measurement signal S is input to the measurement signal storage unit 11 of the first data processing unit 9 and stored.
  • step ⁇ 20 it is determined whether or not the total number of measurement signals S stored in the measurement signal storage unit 11 is equal to or greater than a predetermined number set in advance. That is, in this step, it is determined whether or not the number of data to be signal processed is sufficient. If the total number of measurement signals S is greater than or equal to the predetermined number, the process proceeds to step ⁇ ⁇ ⁇ ⁇ 30. If the total number of measurement signals S is less than the predetermined number, the process returns to step A10. As a result, step ⁇ 10 20 is repeatedly executed until the number of data is sufficient.
  • the basic measurement signal generation unit 2 performs three types of filter processing on the time series data of the measurement signal S: a low-pass filter, a phase filter, and a high-pass filter. Through these filter processes, high acceleration centered on acceleration fluctuations caused by human walking. The frequency and low frequency vibration components are reduced, and the feature of walking time interval is extracted as a wave having a corresponding period.
  • Basic measurement signal S force shown by the broken line in Fig. 2 (a) This wave. Note that the magnitude of the basic measurement signal S is reduced by phase filtering.
  • Time force of 0 Corresponds to the time of the peak position of the time series measurement signal S.
  • the basic data extraction unit 3 performs a basic measurement from the basic measurement signal S.
  • This data D is extracted.
  • basic measurement in which the signal value changes from negative to positive before and after.
  • the measurement signal S corresponding to the zero point of the constant signal S is extracted as basic data D.
  • the extraction range setting unit 4 detects the basic data D.
  • Extraction range A is basic data D
  • the predetermined time t before and the predetermined time t after that are
  • Measured time range Each predetermined time t and t is estimated to be caused by filtering.
  • the peak of the actual measurement signal S is located in the extraction range A. That is, as shown in Fig. 2 (b), the extraction range A can absorb some errors that occur between the actual measurement signal S peak position and the basic data D.
  • the feature data extraction unit 5 extracts the maximum value of the measurement signal S included in the extraction range A as feature data D, and the process proceeds to step A70.
  • the feature data D is the maximum value of the actual measurement signal S, and it represents the walking state.
  • step A70 the data alignment unit 6 divides the feature data D extracted in the previous step, and the time and acceleration when the right foot touches the ground in the person's walking state.
  • the measurement signal S input from the signal detection device 1 is processed in each of the two signal processes.
  • One is signal processing for setting the extraction range ⁇ ⁇ ⁇ ⁇ in the basic measurement signal generator 2 and basic data extraction unit 3, and the other is signal processing for extraction of feature data D in the feature data extraction unit 5.
  • the extraction range A that is set has a width that can tolerate an error.
  • the effects of errors can be offset.
  • feature data D is extracted directly without preprocessing, so accurate feature data D
  • the feature data D is extracted from the original measurement signal without “missing” while referring to the result of the preprocessing, so that the reliability of the data processing can be improved.
  • the feature data extraction unit 5 extracts the maximum value of the measurement signals S included in the extraction range A as the feature data D. This configuration features a lot of noise
  • the filter processing as preprocessing in the basic measurement signal generation unit 2 is general signal processing such as a low-pass filter, a phase filter, and a high-pass filter, which is easy to implement and can be processed in a short time. it can.
  • the calculation processing in the feature data extraction unit 5 does not require a complicated calculation and can obtain a result quickly.
  • the processing content in the first data processing unit 9 of the present embodiment is configured by signal processing using a linear analysis method! /, For example, compared with signal processing using a nonlinear analysis method. Has the advantage of being simple.
  • the data processing apparatus simply extracts the basic data D based on the basic measurement signal S.
  • the configuration includes an extraction range setting unit 4 and a feature data extraction unit 5 that can be simply extracted.
  • an extraction range setting unit 4 As a conventional technique for reducing the error generated by the filter processing, there is a power of finely adjusting the delay amount of the phase change in the phase filter.
  • accurate data is extracted. Therefore, there is no need for such fine adjustment, and there is no need to confirm the basic data D after the actual preprocessing. Therefore, measure
  • the feature data D group input to the second data processing unit 10 has noise. Extracted from raw information that has not been removed
  • each of the first data processing unit 9 and the second data processing unit 10 may be configured as a one-chip microcomputer incorporating a ROM, RAM, CPU, or the like, or as an electronic circuit such as a digital circuit or an analog circuit. It may be formed.
  • the signal processing process from detection of the measurement signal S to output of the result can be automated, so that the small microcomputer as described above is used.
  • a small display device having the same function as the motor 15 according to the present invention or a microsensor having the same function as the signal detection device 1 is mounted. It is also possible to manufacture a small processing device with integrated input / output.
  • the force S to which the acceleration sensor for detecting the acceleration signal is applied as the signal detection means, and the signal to be calculated by the data processing apparatus include various objects. Various parameters related to the state of the can be considered.
  • body movement and respiration rate associated with exercise such as walking, heart rate, body temperature, skin surface temperature, skin potential, pulse wave (pulse rate) ), Electroencephalogram, blood flow, body fluid components such as saliva, respiratory oxygen saturation, blood glucose level, electrocardiogram, electrical conductivity, body weight (pressure on the seating surface), number and cycle of blinks, sweating Amount, electromagnetic wave intensity emitted from the body, chemical substance concentration, etc.
  • the measurement signal S detected by the signal detection device 1 is directly input to the first data processing unit 9, but the signal detection device 1 and the first data processing are configured.
  • a configuration in which the unit 9 is separated may be employed.
  • the time-series data of the measurement signal S detected by the signal detection device 1 is stored in some storage medium, and the time-series data is input to the first data processing unit 9 when calculation processing is required. It is possible.
  • the time series data may be input to the measurement signal storage unit 11, but may be input to each of the basic measurement signal generation unit 2 and the feature data extraction unit 5 without passing through the measurement signal storage unit 11. .
  • the basic measurement signal generation unit 2 is subjected to three types of filter processing: a low-pass filter, a phase filter, and a high-pass filter. May be.
  • the preprocessing in the basic measurement signal generation unit 2 refers to all irreversible (with irreversible changes) arithmetic processing for making it easy to find parameter fluctuations.
  • the specific processing content need not be the filter processing as long as it is an arithmetic processing for making it easy to find parameter fluctuations.
  • Hilbert transform processing, envelope processing, Fourier transform processing, signal processing using an averaging method, wavelet analysis processing, fractal analysis processing, etc. may be used.
  • arbitrary signal addition and subtraction, proportional processing, integration processing, differentiation processing, and the like are examples of filters, and the like.
  • the first data processing unit 9 and the second data processing unit 10 are configured using electronic circuits such as a digital circuit and an analog circuit, instead of the digital filter as described in the above embodiment.
  • an analog filter may be applied.
  • the data processing performed in the first data processing unit 9 may be analog signal processing.
  • the basic data extraction unit 3 uses 0 of the basic measurement signal S.
  • the force with which the measurement signal S corresponding to the point is extracted can be appropriately set according to the calculation target in the data processing apparatus.
  • the position and width of the extraction range A in the extraction range setting unit 4 and the feature data in the feature data extraction unit 5 The same applies to the position where the data D is taken out.
  • the detrend fluctuation analysis method is used in the analysis unit 7, but the analysis method is not limited to this. Considering the general impact of preprocessing on the analysis results, if the analysis method in the analysis unit 7 is a non-linear analysis method, a more accurate analysis is possible compared to the linear analysis method. The result can be expected.
  • Applications of the data processing apparatus, data processing method, and data processing program of the present invention are not particularly limited, and grasp the state of the object based on data obtained by measuring machines, animals and plants, and natural phenomena. Therefore, it can be suitably used as a data processing apparatus, a data processing method, and a data processing program.
  • the present invention is useful for use in extracting a non-linear structure from measured data.

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Description

明 細 書
データ処理装置,データ処理方法及びデータ処理プログラム
技術分野
[0001] 本発明は、機械,動植物や微生物等の生命体の状態や天候や地震等の自然現象 といった、種々の対象体の状態を分析するためのデータ処理装置,データ処理方法 及びデータ処理プログラムに関するものである。
背景技術
[0002] 従来、対象体の状態に応じて変動するパラメータを時系列の測定信号として測定し
、その測定信号の中から特徴的な信号を抽出することで、対象体の状態を分析する データ処理技術が多数開発されてレ、る。 例えば、特許文献 1には、被験者の脈波波形を測定信号として採取するとともに、 その測定信号の中に内包されている秩序 (すなわち、測定信号の変動を支配する決 定論的な構造)を抽出するための演算を行い、被験者の状態を診断する技術が開示 されている。この技術では、測定信号のカオスアトラクター及びリアプノフ指数を演算 することで、測定信号に内在する秩序を論理的に抽出することができ、被験者に対す る客観的な診断が可能になるとされている。
[0003] また、特許文献 2には、被験者から検出される生体信号のデータを解析することに よって被験者の健康状態を診断する技術が開示されている。この技術におけるデー タの解析手法としては、カオス解析のほか、デトレンド変動解析 (DFA)や周波数変 換,ウェーブレット解析,マルチフラクタル解析等が挙げられている。
ところで、これらのような従来の技術では一般的に、測定信号に対しその特徴を把 握しやすくするための前処理としての種々の信号処理を施した後、実質的なデータ 処理が実施されるようになっている。この前処理には、例えば測定信号に混入してい るノイズを取り除くためのノイズ低減処理(ノイズ除去処理)や、測定信号のうち所定の 周波数成分を取り出すフィルタ処理 (フィルタリング),直交変換処理,フーリエ変換 処理,エンベロープ処理等がある。
[0004] 上述の特許文献 1に記載の技術にお!/、ても、脈波センサで検出された測定信号を A/D変換器に通してデジタル信号へと変換する際に、信号の演算処理速度を向上 させることを目的として、 A/D変換器力 出力される離散データを整数型に限定す る前処理が実施されると記載されて!/、る。これらのような前処理を施すことによって、 その後の演算処理においてあまり重要ではない情報を取り除くことができ、データ処 理の速度や精度を高めることができるようになつている。
特許文献 1:特公平 6— 9546号公報
特許文献 2:特開 2001— 299766号公報
発明の開示
発明が解決しょうとする課題
[0005] しかしながら、前処理の過程で取り除かれる情報には、単なるノイズや不要な情報 だけでなぐその後の演算処理や演算結果に影響を与える重要な情報が含まれてい る場合がある。
例えば、前述の種々の前処理のうち、不可逆的な演算が含まれる前処理を実施し た場合には、その処理内容に関わらず、測定信号中に含まれる何らかの情報が失わ れることになる。つまり、前処理に伴う特徴的な情報の抜け (情報の脱落)によって、 その後の演算処理における演算精度が低下してしまい、演算結果が不正確となるお それがある。
[0006] また、測定信号中の非線形構造を抽出するための非線形解析手法を用いる場合 には、より本質的な課題が生じる。
すなわち、動植物や自然現象を測定対象とした通常の測定信号中には、不規則な 挙動を示すカオスが存在している。カオスとは、その測定対象に内在する非線形の 決定論的システムに起因する予測不能の複雑な振る舞いのことであって、ノイズ (決 定論的なシステムに依らない不規則な挙動であり、処理対象となる情報以外の不要 な情報)とは概念的に異なるものではあるものの、厳密にこれらを区別することは現状 では極めて困難とされて!/、る。
[0007] そのため、例えば測定信号中から一般的なノイズを除去する処理を行ってしまえば 、ノイズと共にカオスの情報が部分的に除去されてしまうことは免れ得ない。つまり、 測定信号に前処理を施すことによって測定信号中の非線形構造も部分的に取り除 かれてしまい、却ってその構造を把握しに《なってしまう。線形解析手法においてノ ィズとして切り捨てられていた部分の情報が重要な意味を持つ非線形解析手法では 、前処理によって解析結果に与えられる影響力 線形解析手法におけるそれとは比 較にならなレ、程大きレ、のである。
[0008] 一方、上記のような前処理を全く行わなければ、特徴的な情報の抜けは生じないも のの、不要な情報を含んだ全ての測定信号を演算処理することになるため、多大な 処理労力と処理時間とが必要となる。
本発明はこのような課題に鑑み案出されたもので、簡素な構成で、データ処理速度 及びデータ処理精度を向上させることができるデータ処理装置,データ処理方法及 びプログラムを提供することを目的とする。また、測定信号中の非線形構造を抽出す る演算処理において、簡素な構成で、短時間で正確な演算結果を得られるようにし た、データ処理装置,データ処理方法及びデータ処理プログラムを提供することを目 的とする。
課題を解決するための手段
[0009] 上記目的を達成するために、請求項 1記載の本発明のデータ処理装置は、対象体 の状態に応じて変動するパラメータを測定信号として検出する測定信号検出手段( 信号検出装置)と、該測定信号検出手段で検出された該測定信号に対し、該パラメ ータの変動を把握するための前処理としての信号処理を施して、基本測定信号を生 成する信号処理手段 (基本測定信号生成部)と、該信号処理手段で生成された該基 本測定信号に基づレ、て、該パラメータの変動を特徴付ける測定信号を基本データと して抽出する基本データ抽出手段 (基本データ抽出部)と、該基本データ抽出手段 で抽出された該基本データによって規定される所定の領域を該測定信号の抽出範 囲として設定する抽出範囲設定手段 (抽出範囲設定部)と、該測定信号検出手段で 検出された該測定信号のうち、該抽出範囲設定手段で設定された該抽出範囲に含 まれる測定信号の中から、該対象体の状態を特徴付ける測定信号を特徴データとし て抽出する特徴データ抽出手段(特徴データ抽出部)とを備えたことを特徴としてい
[0010] なお、該信号処理手段における該前処理とは、該パラメータの変動を見つけやすく するための不可逆的な(非可逆変化を伴う)演算処理全般のことを指しており、例え ば該パラメータ中に含まれる不必要な情報(ノイズ)を除去するためのフィルタ処理や ヒルベルト変換処理,エンベロープ処理,フーリエ変換処理,加算平均の手法を用い た信号処理,ウェーブレット解析処理,フラクタル解析処理等を含む。また、任意の信 号加算や減算,比例処理,積分処理,微分処理等も含む。
[0011] また、請求項 2記載の本発明のデータ処理装置は、請求項 1記載の構成において 、該測定信号検出手段が、動物のバイタルサインを該測定信号として検出することを 特徴としている。
なお、ここでいうバイタルサインとは、動物(人を含む)の身体から検出される生命徴 候としての物理量のことを意味している。例えば、歩行等の運動に伴う体動や呼吸数 ,心拍数,体温,皮膚表面温度,皮膚電位,脈波 (脈拍数),脳波,血流量,唾液な どの体液成分,呼吸気中や血中の酸素飽和度,血糖値,心電,電気伝導度,体重( 着座面への圧力),まばたきの数や周期,発汗量,その他身体から発せられる電磁 波の強度や化学物質濃度等が挙げられる。
[0012] また、請求項 3記載の本発明のデータ処理装置は、請求項 1又は 2記載の構成に おいて、該信号処理手段が、線形解析手法を用いて該測定信号に対し該信号処理 を施すことを特徴としている。
また、請求項 4記載の本発明のデータ処理装置は、請求項;!〜 3の何れか 1項に記 載の構成において、該信号処理手段が、該測定信号検出手段で検出された該測定 信号から、予め設定された所定の周波数成分を濾波するフィルタ処理手段(基本測 定信号生成部)を有してなることを特徴として!/、る。
[0013] また、請求項 5記載の本発明のデータ処理装置は、請求項;!〜 4の何れか 1項に記 載の構成において、該測定信号検出手段が、該測定信号として略周期的に変動す るパラメータを検出し、該信号処理手段が、該測定信号を平滑化した波動を該基本 測定信号として生成するとともに、該基本データ抽出手段が、該基本測定信号にお ける変動ピークの検出時刻を該基本データとして抽出することを特徴としている。ここ で、「該測定信号を平滑化する」とは、測定信号データの波形を滑らかな波動 (波形) にすることをいう。 [0014] また、請求項 6記載の本発明のデータ処理装置は、請求項 5記載の構成において 、該抽出範囲設定手段が、該ピークの検出時刻の近傍時刻を該抽出範囲として設 定することを特 ί毁としている。
また、請求項 7記載の本発明のデータ処理装置は、請求項 6記載の構成において 、該特徴データ抽出手段が、該抽出範囲に含まれる測定信号における変動ピークの 検出時刻を該特徴データとして抽出することを特徴としている。
[0015] また、請求項 8記載の本発明のデータ処理装置は、請求項;!〜 7の何れか 1項に記 載の構成において、該特徴データ抽出手段で抽出された該特徴データに基づき、 該測定信号中の非線形構造を抽出し、該対象体の状態を解析する演算処理手段( 第二データ処理部)をさらに備えたことを特徴としている。
請求項 9記載の本発明のデータ処理方法は、対象体の状態に応じて変動するパラ メータを測定信号として検出する測定信号検出ステップと、該測定信号検出ステップ で検出された該測定信号に対し、該パラメータの変動を把握するための前処理とし ての信号処理を施して、基本測定信号を生成する信号処理ステップと、該信号処理 ステップで生成された該基本測定信号に基づレ、て、該パラメータの変動を特徴付け る測定信号を基本データとして抽出する基本データ抽出ステップと、該基本データ抽 出ステップで抽出された該基本データによって規定される所定領域を抽出範囲とし て設定する抽出範囲設定ステップと、該測定信号検出ステップで検出された該測定 信号のうち、該抽出範囲設定ステップで設定された該抽出範囲に含まれる測定信号 の中から、該対象体の状態を特徴付ける測定信号を特徴データとして抽出する特徴 データ抽出ステップとを備えたことを特徴としている。
[0016] また、請求項 10記載の本発明のデータ処理方法は、請求項 9記載の構成に加え、 該測定信号検出ステップにおレ、て、動物のバイタルサインを該測定信号として検出 することを特 ί毁としている。
また、請求項 11記載の本発明のデータ処理方法は、請求項 9又は 10記載の構成 に加え、該信号処理ステップにおいて、線形解析手法を用いて該測定信号に対し該 信号処理を施すことを特徴として!/、る。
[0017] また、請求項 12記載の本発明のデータ処理方法は、請求項 9〜; 11の何れか 1項に 記載の構成に加え、該信号処理ステップにおいて、該測定信号検出ステップで検出 された該測定信号から、予め設定された所定の周波数成分を濾波することを特徴と している。
また、請求項 13記載の本発明のデータ処理方法は、請求項 9〜; 12の何れか 1項に 記載の構成に加え、該測定信号処理ステップにおいて、該測定信号として略周期的 に変動するパラメータを検出し、該信号処理ステップにおいて、該測定信号を平滑 化した波動を該基本測定信号として生成するとともに、該基本データ抽出ステップに おいて、該基本測定信号における変動ピークの検出時刻を該基本データとして抽出 することを特 ί毁としている。
[0018] また、請求項 14記載の本発明のデータ処理方法は、請求項 13記載の構成に加え 、該抽出範囲設定ステップにおいて、該ピークの検出時刻の近傍時刻を該抽出範囲 として設定することを特 ί毁として!/ヽる。
また、請求項 15記載の本発明のデータ処理方法は、請求項 14記載の構成に加え 、該特徴データ抽出ステップにおいて、該抽出範囲に含まれる測定信号における変 動ピークの検出時刻を該特徴データとして抽出することを特徴としている。
[0019] また、請求項 16記載の本発明のデータ処理方法は、請求項 9〜; 15の何れか 1項に 記載に加え、該特徴データ抽出ステップで抽出された該特徴データに基づき、該測 定信号中の非線形構造を抽出し、該対象体の状態を解析する演算処理ステップをさ らに備えたことを特徴としている。
請求項 17記載の本発明のデータ処理プログラムは、コンピュータを、信号処理手 段,基本データ抽出手段,抽出範囲設定手段及び特徴データ抽出手段として機能 させるためのデータ処理プログラムであって、該信号処理手段が、測定信号として検 出された、対象体の状態に応じて変動するパラメータの変動を把握するための前処 理としての信号処理を該測定信号に対して施すとともに、基本測定信号を生成し、該 基本データ抽出手段が、該信号処理手段で生成された該基本測定信号に基づ!/、て 、該パラメータの変動を特徴付ける測定信号を基本データとして抽出し、該抽出範囲 設定手段が、該基本データ抽出手段で抽出された該基本データによって規定される 所定の領域を該測定信号の抽出範囲として設定し、該特徴データ抽出手段が、該 測定信号検出手段で検出された該測定信号のうち、該抽出範囲設定手段で設定さ れた該抽出範囲に含まれる測定信号の中から、該対象体の状態を特徴付ける測定 信号を特徴データとして抽出することを特徴としている。
発明の効果
[0020] 本発明のデータ処理装置,データ処理方法及びデータ処理プログラム(請求項 1 , 9, 17)によれば、測定信号の抽出範囲の設定においては前処理としての信号処理 を施した測定信号を用い、より具体的な特徴データの抽出においては、その抽出範 囲に含まれる測定信号から取り出すため、対象体の状態を特徴付ける測定信号を極 めて正確に抽出することができる。
[0021] 本発明のデータ処理装置及びデータ処理方法 (請求項 2, 10)によれば、ノイズが 多く特徴の抽出が困難とされているような動物のバイタルサインの中からでも、正確 に特徴データを抽出することができる。
本発明のデータ処理装置及びデータ処理方法 (請求項 3, 11)によれば、前処理と しての信号処理が容易である。また、簡素な構成で短時間に処理を済ませることがで きる。
[0022] 本発明のデータ処理装置及びデータ処理方法 (請求項 4, 12)によれば、容易に 信号処理を行うことができる。また、パラメータの変動を特徴付ける基本データを抽出 するのに十分な情報を素早く選別することができる。
本発明のデータ処理装置及びデータ処理方法 (請求項 5, 13)によれば、測定信 号を平滑化して生成された基本測定信号のピークを抽出するようになっているため、 容易に基本データを検出することができる。
[0023] 本発明のデータ処理装置及びデータ処理方法 (請求項 6, 14)によれば、基本測 定信号のピークの検出時刻の近傍時刻を抽出範囲として設定するため、基本測定 信号のピークから離れている範囲に含まれる測定信号を、特徴データの抽出対象か ら除外することができる。つまり、対象体の状態を特徴付ける測定信号との相関が強 V、部分の情報を容易に取り出すこと力 Sできる。
[0024] 本発明のデータ処理装置及びデータ処理方法 (請求項 7, 15)によれば、測定信 号中に含まれる、対象体の状態を特徴付ける測定信号を正確に取り出すことができ 本発明のデータ処理装置及びデータ処理方法 (請求項 8, 16)によれば、正確に 抽出された特徴データに基づいて測定信号中の非線形構造を抽出することができ、 信頼性の高いデータ解析を行うことができる。
図面の簡単な説明
[0025] [図 1]本発明の一実施形態に係るデータ処理装置の全体構成を示すブロック図であ
[図 2]本データ処理装置におけるデータ処理内容を説明するためのグラフであり、 (a )は基本測定信号生成部及び基本データ抽出部におけるデータ処理に係る測定信 号の時系列グラフ、 (b)は特徴データ抽出部におけるデータ処理に係る測定信号の 時系列グラフである。
[図 3]本データ処理装置における制御内容を示すフローチャートである。
[図 4]コンピュータを利用した本データ処理装置の構成例を示す模式図である。 符号の説明
[0026] 1 信号検出装置 (測定信号検出手段)
2 基本測定信号生成部 (信号処理手段)
3 基本データ抽出部 (基本データ抽出手段)
4 抽出範囲設定部 (抽出範囲設定手段)
5 特徴データ抽出部 (特徴データ抽出手段)
6 データ整列部
7 解析部
8 判定部
9 第一データ処理部
10 第二データ処理部(演算処理手段)
11 測定信号記憶部
12 コンピュータ
13 記憶装置
14 中央処理装置(CPU) 15 モニタ
発明を実施するための最良の形態
[0027] 以下、図面により、本発明の実施の形態について説明する力 本発明は以下の実 施の形態に限定されるものではなぐその要旨を逸脱しない範囲内であれば種々に 変更して実施することができる。
[0028] [1.構成]
[1 - 1.全体構成]
本実施形態では、人間の歩行状態を解析対象としたデータ処理装置を具体例とし て説明する。すなわち、本データ処理装置は、人間の歩行状態に対応するパラメ一 タ(例えば、加速度)を検出信号として検出し、それにデータ処理を施して出力を行う 装置である。
[0029] 図 1に示すように、本データ処理装置は、信号検出装置 (測定信号検出手段) 1 ,第 一データ処理部 9及び第二データ処理部 10を備えて構成される。第一データ処理 部 9は、信号検出装置 1で検出された信号の特徴を把握しやすくするための演算処 理を施すものであり、一方、第二データ処理部 10は、実質的なデータ処理として歩 行状態を解析するものである。
これらの第一データ処理部 9及び第二データ処理部 10は、コンピュータの内部で 演算処理される機能部位であり、各機能は個別のプログラムとして構成されている。 なお、本実施形態における第一データ処理部 9は、信号処理手段,基本データ抽出 手段,抽出範囲設定手段及び特徴データ抽出手段として機能するものである。また 、第二データ処理部 10は、演算処理手段として機能するものである。
[0030] コンピュータを利用した本データ処理装置の構成例を図 4に示す。このコンピュータ 12は、上述の信号検出装置 1 ,記憶装置 (ROM, RAM等) 13, 中央処理装置(CP U) 14,出力インタフェースとしてのモニタ 15,入力インタフェースとしてのキーボード 16及びマウス 17を備えて構成されている。ここで、本データ処理装置に係る第一デ ータ処理部 9及び第二データ処理部 10は、記憶装置 13の内部にプログラムとして記 憶されている。
[0031] 以下、本データ処理装置における信号処理内容について、図 1のブロック図を用い て概念的に説明する。
[1 - 2.信号検出装置]
信号検出装置 1は、機械,動植物や微生物等の生命体の状態や天候や地震等の 自然現象と!/、つた、種々の対象体の状態に関わる様々なパラメータ(変動要素)を検 出するセンサである。このパラメータには、センサから直接検出される情報のほか、セ ンサでの検出情報を演算等によって処理して、対応するパラメータの値を推定値とし て求めたものも含まれる。
[0032] 本実施形態ではこの信号検出装置 1として、本データ処理装置の演算対象となる 加速度信号を検出するための加速度センサが適用されており、対象となる人物の体 に装着されている。なお、この加速度センサは、測定対象や目的に合わせて一軸〜 三軸のものを任意に用いてよいが、歩行時における鉛直方向,水平前後方向及び 水平左右方向の三方向へ作用する加速度を検出するための三軸加速度センサを用 いるのが好ましい。今回の具体例では三軸加速度センサを用いており、ここで検出さ れた鉛直方向の加速度の検出情報及びその検出時刻情報が、図 1に示すように、測 定信号 Sとして第一データ処理部 9へ入力されるようになっている。
[0033] [1 3·第一データ処理部]
第一データ処理部 9は、本願請求項 1に規定された処理を測定信号 Sに施す機能 部位であり、図 1に示すように、測定信号記憶部 11 ,基本測定信号生成部 2,基本デ ータ抽出部 3,抽出範囲設定部 4及び特徴データ抽出部 5を備えて構成される。ここ で施される処理とは、光学信号,音声信号,電磁気信号等からなる測定信号 Sに対し 、その特徴を把握しやすくするために、数学的,電気的な加工を施して信号を変換 すること(換言すれば、信号処理)である。なお、該処理はその処理対象となる信号の 種類によって、アナログ信号処理とデジタル信号処理とに分類することができる。本 第一データ処理部 9は、デジタル信号処理の範疇に含まれる処理を実施するもので ある。
[0034] [1 - 3- 1.測定信号記憶部]
測定信号記憶部 11は、信号検出装置 1から入力された測定信号 Sを記憶する機能 部位である。ここに記憶された測定信号 Sは、図 1に示すように、測定信号記憶部 11 から二系統に分かれて基本測定信号生成部 2及び特徴データ抽出部 5のそれぞれ へと入力されるようになっている。つまり、基本測定信号生成部 2及び特徴データ抽 出部 5の各々に対して、何ら加工されていない生の情報が入力されることになる。
[0035] なお、この測定信号記憶部 11に入力される測定信号 Sは、信号検出装置 1におい て所定時間の間に検出された一連の時系列データであってもよいし、あるいは、信号 検出装置 1で随時検出された個別の測定データであってもよい。
前者の場合には、時系列データとしての測定信号 Sがそのまま、基本測定信号生 成部 2及び特徴データ抽出部 5へ入力される。また、後者の場合には、信号検出装 置 1からの測定信号 Sの総数が予め設定された所定数以上となるまでの間、各測定 信号 Sが測定信号記憶部 11に記憶され、データ数が十分に揃った段階でそれら全 体の測定信号 Sが基本測定信号生成部 2及び特徴データ抽出部 5へと入力されるよ うになつている。なお、本実施形態では、後者の場合について説明する。
[0036] [1 3— 2.基本測定信号生成部]
基本測定信号生成部 2は、信号検出装置 1から入力された時系列の測定信号 Sか ら必要な信号成分を取り出すための処理を実施する機能部位である。本具体例では 、歩行時の加速度変化のピークを検出するためのデジタルフィルタ処理がなされて おり、加速度の時系列データに対し、ローパスフィルタ,位相フィルタ及びハイパスフィ ルタの三種のフィルタ処理が施されるようになって!/、る。
[0037] この基本測定信号生成部 2で施される処理の具体的な内容としては、測定信号デ ータの種類や内容、あるいは、どのような目的のもとに解析するのかによつて異なるが 、一般的には、フィルタ処理,フーリエ変換処理,エンベロープ処理といった線形解 析手法を単独で、又は、それらを組み合わせて行うのが好ましい。
ローパスフィルタは、多様な周波数成分を含んだ入力信号の中から高周波数の振 動成分を低減させる(あるいは除去する)フィルタである。このローパスフィルタにより、 人物の歩行に伴う加速度変動の周波数 (すなわち、一方の足が地面についてから他 方の足が地面につくまでの時間間隔の逆数)よりも高周波数のノイズ成分が抑制され るようになっている。なお、図 2 (a)中に実線で示された測定信号 Sに見られる細かい 時間間隔の振動が、高周波数のノイズ成分である。 [0038] また、位相フィルタは、入力信号の位相変化を遅延させるフィルタである。ここでは、 入力信号の位相をその加速度変動の四半周期分だけ遅延させるようになつている。 一方、ハイパスフィルタは、人物の歩行に伴う加速度変動の周波数よりも低周波数 の振動成分を低減させる(ある!/、は除去する)フィルタである。このハイパスフィルタに より、信号検出装置 1で検出された信号のドリフト成分が抑制されるようになっている。 これらのフィルタ処理を施すことにより、測定信号 Sの中から人物の歩行に伴う加速度 変動の時系列データが明確化されるようになつている。
[0039] 例えば、図 2 (a)中に実線で示された測定信号 Sとしての加速度の時系列信号に対 して、上記の各フィルタ処理を施すと、破線で示されたような信号が取り出される。つ まり、信号検出装置 1で検出された生の情報の中から、歩行の時間間隔という特徴が 、それに対応する周期を有する波として取り出される。こうして取り出された時系列信 号のことを、以下、基本測定信号 Sと呼ぶ。
B
[0040] なお、この基本測定信号 Sは位相フィルタを介して取り出された信号であるため、こ
B
の信号の大きさが 0となる時刻(すなわち 0点の位置) 1S 元の入力信号のピーク位置
(最大値,最小値)が検出された時刻に対応している。
[0041] [1 3— 3·基本データ抽出部]
基本データ抽出部 3は、基本測定信号生成部 2で生成された基本測定信号 Sの中
B
から、加速度変動を特徴づける情報を抽出するものである。ここでは、加速度変動を 特徴づける情報として、フィルタ処理によって得られた基本測定信号 Sから推定され
B
る、入力された測定信号 Sのピーク位置を抽出するようになっている。つまり、基本デ ータ抽出部 3では、基本測定信号 Sの 0点のうち、その前後で信号値が負から正へと
B
変化するものに着目し、その 0点に対応する測定信号 Sを抽出する。より具体的には 、図 2 (a)中において 0点の時刻と同一時刻の測定信号 Sを基本データ Dとして抽出
B
するようになつている。
[0042] なお、この図 2 (a)に示すように、実際の測定信号 Sのピーク位置と基本データ Dと
B
の間には若干のズレが生じていることがわかる。これは、前述の通り、基本測定信号 生成部 2におけるフィルタ処理によって失われた情報によるものである。すなわち、フ ィルタリングを施すことで測定信号 S中に含まれる特徴的な情報に抜けが生じた結果 、基本測定信号 sの 0点の位置と実際のピーク位置とが相違してしまっているのであ
B
[0043] 一般に、入力信号に対する前処理として不可逆的な演算処理操作を行うほど、そ の処理内容に関わらず、入力信号中に含まれる本来の情報が失われる。つまり、たと えそれが入力信号の変動を見つけやすくするための処理操作であったとしても、そ の結果得られる情報には何らかの誤差が含まれてしまうことになる。もちろん、入力信 号の変動把握が容易となる限度内において可能な限り誤差の小さい前処理を追求 することも考えられるものの、その結果得られる演算の精度には限界がある。
[0044] また、このような誤差を小さくするための典型的な手法として、フィルタ処理内容を 調整することで測定信号 Sのピーク位置と基本データ Dとを一致させるとレ、う手法が
B
ある。例えば、位相フィルタにおける位相変化の遅延量を微調整して、測定信号 Sの ピーク位置と基本データ Dとの誤差を一様に狭めることが考えられる。
B
しかしながら、このような微調整は、熟練した技術者の勘に頼らざるを得ないやや不 確実な手法であり、データ処理精度が不安定となるおそれがある。また、このような微 調整は、対象とする個々のデータ群毎 (例えば、信号を検出した個体毎)にそれぞれ に応じて行う必要があるため、複数の(大量の)データ群を扱う場合、処理に時間が かかり、調整作業も繁雑なものとなる。さらに、実際に前処理を施した後の基本デー タ Dを確認してからでなければ調整ができないため、このような手法では十分なデー
B
タ処理速度が得られな!/、のである。
[0045] そこで、本発明では、単に基本測定信号 Sに基づいて基本データ Dを抽出するだ
B B
けでなぐ以下に説明する抽出範囲設定部 4及び特徴データ抽出部 5を備えた構成 としている。
[0046] [1 3— 4.抽出範囲設定部]
抽出範囲設定部 4は、基本データ Dが検出された時刻の近傍時間を抽出範囲 Aと
B
して設定するものである。近傍時間とは基本データ Dが検出された時刻の前後の時
B
間のことを意味している。ここでは図 2 (b)に示すように、抽出範囲 Aが、基本データ D が検出された時刻以前の所定時間 tと基本データ Dが検出された時刻以後の所定
B 1 B
時間 とから構成されている。なお、ここで設定された抽出範囲 Aは、続いて説明する 特徴データ抽出部 5へと入力されるようになっている。
[0047] なお、所定時間 t及び tは、それぞれ任意に設定可能な時間である。例えば、基本
1 2
測定信号 Sの波長えに対する割合として定めてもよいし、予め設定した値としてもよ
B
い。
本実施形態では、各所定時間 t及び t力 上述のフィルタ処理によって生じると推定
1 2
される時間誤差よりも大きくなるように設定されている。なお、フィルタ処理によって失 われる情報量が比較的少ない場合には所定時間 t及び tが比較的短くてもよぐ一
1 2
方、失われる情報量が比較的多いフィルタ処理を実施する場合には所定時間 t及び
1 tを比較的長く設定するとよい。
[0048] [1 3— 5·特徴データ抽出部]
特徴データ抽出部 5は、抽出範囲設定部 4で設定された抽出範囲 A内に含まれる 測定信号 Sの中から、人物の歩行状態を特徴付ける信号、すなわち、実際のピーク 位置を特徴データ Dとして抽出するものである。ここでの演算処理は、図 2 (b)に示
C
すように、抽出範囲 A内の測定信号 Sの最大値を検出することで求められている。こ れにより、歩行状態において左右何れかの足が地面についた時刻及びその時の鉛 直方向の加速度が、特徴データ Dとして正確に抽出されることになる。なお、図 2 (b)
C
中においては、特徴データ Dのグラフ上の位置が記号 +で示されている。
C
[0049] このように、抽出範囲設定部 4及び特徴データ抽出部 5は、基本データ Dに基づい
B
て再び元の測定信号 Sへ立ち返り、測定信号 Sの中から特徴データ Dを抽出するよ
C
うに機能する。つまり、基本データ Dには誤差が含まれているものの、その近傍に本
B
来の特徴を示す情報が存在するものと見なして抽出範囲設定部 4で抽出範囲 Aを設 定している。さらに、特徴データ抽出部 5では、抽出範囲 A中に含まれる測定信号 S のみを演算対象とすることで演算労力や演算時間を低減させるとともに、測定信号 S から特徴データ Dを抽出することでその精度を確保しているのである。
C
[0050] [1 4·第二データ処理部]
第二データ処理部 10は、第一データ処理部 9において処理が施されたデータに対 する実質的なデータ処理を行うための機能部位であり、図 1に示すように、データ整 列部 6,解析部 7及び判定部 8を備えて構成される。この第二データ処理部 10では、 歩行時における左右何れか一方の足の加速度データのピーク間隔時間の特徴が解 析されるようになつている。なお、本実施形態における解析の手法としては、ピーク間 隔時間の揺らぎを観察する非線形解析手法が用いられてレ、る。
[0051] ここでいう「揺らぎ」とは、ある波動が刻々と変化する際に観察される僅かな波形の ズレ (空間的,時間的変化や動きが部分的に不規則な動き)のことを指している。例 えば、歩行に伴う体動を加速度変動として検出した時系列データだけでなぐ呼吸数 や心拍数,脳波等のバイタルサインを時系列データとした場合にも、それらの波動の ピーク間隔や周期は一定ではなぐ複雑な変動を示すことが知られている。一方で、 このような不規則に見える複雑な変動の中から、その挙動を支配していると考えられ る構造を解析するための数々の手法が提案されている。第二データ処理部 10は、こ れらのような手法を利用して、ピーク間隔時間の揺らぎの度合いを観察することにより 、その変動の背後に存在するであろう非線形構造を解析するものである。
[0052] なお、具体的な解析手法としては、スペクトル解析(FFT解析) ,フラクタル解析(マ ルチフラクタル解析,デトレンド変動解析等),カオス解析及びウェーブレット解析等 の公知の解析手法が挙げられる力 ここでは、フラクタル解析法の一つであるデトレ ンド変動解析が用いられてレ、る。
デトレンド変動解析の手法は、解析対象となる波動の複雑性をスケーリング指数と 呼ばれる値で評価する統計的な解析手法である。本実施形態では、左右何れか一 方の足の加速度データがデータ整列部 6において整列され、その加速度データのス ケーリング指数が解析部 7において演算され、さらにその評価が判定部 8でなされる ようになつている。
[0053] [1 4 1 ·データ整列部]
まず、データ整列部 6は、解析部 7における演算の便宜を図るベぐ特徴データ抽 出部 5で抽出された特徴データ Dを整歹 IJさせるものである。つまり、特徴データ Dを
C C
どのように整列させるかは、解析部 7におけるデータ処理の種類等に応じて適宜設定 される。なお、解析部 7での解析手法に応じて、このデータ整列部 6における演算を 省略してもよい。
[0054] 本実施形態では、特徴データ抽出部 5で抽出された特徴データ D において、人物 の左右の足が区別されていないため、歩行状態をより詳細に観察するために、この データ整列部 6において特徴データ Dが交互に並べ替えられ、分割されるようにな
C
つている。
具体的には、抽出された特徴データ Dをその検出時刻の順に並べ、奇数番目のも
C
のと偶数番目のものとに分離する。つまり、例えば奇数番目の特徴データ D群には
C
、左右何れか一方の足が地面についた時刻及びその時の加速度のデータが羅列さ れ、偶数番目の特徴データ D群には、他方の足のデータが羅列されるようになって
C
いる。これにより、一方の特徴データ D群の時間間隔を演算することで、片足の歩行
C
間隔(片足が地面に接触する時間間隔)を把握することができるようになつている。
[0055] [1 4 2·解析部]
解析部 7は、第二データ処理部 10のデータ整列部 6から入力された偶数番目及び 奇数番目の特徴データ D群のそれぞれのデータ群について個別にスケーリング指
C
数を演算する。具体的には、片方の特徴データ D群をその検出時刻に基づいて n
C
個の区間に分割し、各区間において各特徴データ Dが検出された時間間隔 (歩行
C
間隔時間)とそのトレンドとの最小二乗誤差 (分散) Fを算出して、分割数 n及び分散 F の各々の対数プロットの勾配 αをスケーリング指数として演算する。なお、トレンドとは 、各区間内におけるデータの推移傾向を意味しており、例えば各区間内のデータを 直線に近似したものとする。
[0056] この方法では、区間の分割数 ηを変化させれば、特徴データ D群の観察スケール
C
も変化することになり、算出される分散 Fも変化することになる。一方、分割数 η及び 分散 Fの各々の対数プロットに線形関係が認められれば、分割数 ηと分散 Fとの間に は自己相似におけるスケールが存在するということになる。つまりここでは、観察する 区間を変化させた場合における、実際の歩行間隔時間のばらつきの度合いの自己 相似性の大きさをスケーリング指数 αとして演算していることになる。
[0057] [1 4 3·判定部]
判定部 8は、解析部 7で演算されたスケーリング指数 αに基づいて、歩行状態を判 定する。一般に、スケーリング指数 αの値によって、分割数 ηと分散 Fとの間の相関を 判断することができることが知られている。例えば、 0. 5< α < 1である場合には、長 距離相関が認められ、 α = 1である場合には 1/f揺らぎの相関が認められる。また、 相関はあるもののフラクタル性が認められない場合には、 α〉1となる。
[0058] なお、「l/f揺らぎ」とは、前述の揺らぎのうち、揺らぎ成分の大きさ(パワースぺタト ノレ)が周波数 fに対して 1/fとなるような揺らぎのことを意味している。例えば、小川の せせらぎ音やそよ風の風圧,木目の形状,小鳥のさえずり音といった自然界に存在 する波動をスペクトル解析すると、パワースペクトルが周波数 fに反比例する 1/f揺ら ぎを観測することができる。近年では、人間や動物などから発せられる生体信号 (バイ タルサイン)にも揺らぎが観察されることが判明しており、特に、観察対象の健康状態 を判断するための指標として 1/f揺らぎを用いることの有用性が多数報告されている
[0059] これらのような特性に基づき、判定部 8は、スケーリング指数 αが α = 1に近いほど 、良好な歩行状態であると判定するようになっている。ここでの判定結果は前述のモ ユタ 15 出力されるようになっている。
[0060] [2.フローチャート]
図 3に示すフローチャートを用いて、本データ処理装置における制御内容を説明す ステップ A10では、信号検出装置 1としての加速度センサにより、加速度の検出情 報及びその検出時刻情報が測定信号 Sとして検出される。ここで検出された測定信 号 Sは、第一データ処理部 9の測定信号記憶部 11へ入力され、記憶される。
続くステップ Α20では、測定信号記憶部 11にお!/、て、記憶された測定信号 Sの総 数が予め設定された所定数以上であるか否かが判定される。つまりこのステップでは 、信号処理すべきデータ数が十分に揃っているか否かが判定される。ここで、測定信 号 Sの総数が所定数以上である場合にはステップ Α30へ進み、測定信号 Sの総数が 所定数未満である場合にはステップ A10 戻る。これにより、データ数が十分に揃う までの間、ステップ Α10 20が繰り返し実行されることになる。
[0061] ステップ Α30では、基本測定信号生成部 2において、測定信号 Sの時系列データ に対しローパスフィルタ,位相フィルタ及びハイパスフィルタの三種のフィルタ処理が 施される。これらのフィルタ処理により、人物の歩行に伴う加速度変動を中心として高 周波数及び低周波数の振動成分が低減され、歩行の時間間隔という特徴が、それ に対応する周期を有する波として取り出される。図 2 (a)に破線で示される基本測定 信号 S力 この波である。なお位相フィルタ処理により、基本測定信号 Sの大きさが
B B
0となる時刻力 時系列の測定信号 Sのピーク位置の時刻に対応するものとなる。
[0062] 続くステップ A40では、基本データ抽出部 3において、基本測定信号 Sの中から基
B
本データ Dが抽出される。ここでは、前後で信号値が負から正へと変化する基本測
B
定信号 Sの 0点に対応する測定信号 Sが、基本データ Dとして抽出される。図 2 (a)
B B
に示すように、実際の測定信号 Sのピーク位置と基本データ Dとの間には、若干の
B
誤差が生じて!/、ること力 Sわ力、る。
[0063] さらに続くステップ A50では、抽出範囲設定部 4において、基本データ Dが検出さ
B
れた時刻の近傍時刻が抽出範囲 Aとして設定される。抽出範囲 Aは、基本データ D
B
が検出された時刻を挟んで、それ以前の所定時間 tとそれ以後の所定時間 tとを合
1 2 計した時間の幅である。各所定時間 t及び tは、フィルタ処理によって生じると推定さ
1 2
れる時間誤差よりも大きく設定されているため、抽出範囲 Aの中に実際の測定信号 S のピークが位置することになる。つまり、図 2 (b)に示すように、抽出範囲 Aは、実際の 測定信号 Sのピーク位置と基本データ Dとの間に生じている若干の誤差を吸収しうる
B
幅を備えている。
[0064] 続くステップ A60では、特徴データ抽出部 5において、抽出範囲 A内に含まれる測 定信号 Sの最大値が特徴データ Dとして抽出され、ステップ A70へと進む。図 2 (b)
C
に示すように、特徴データ Dは実際の測定信号 Sの最大値であり、歩行状態を特徴
C
付ける正確な信号となる。
ステップ A70では、データ整列部 6において、前ステップで抽出された特徴データ Dが分割され、人物の歩行状態における右足が地面についた時刻及び加速度の時
C
系列データと、左足に係る時系列データとが生成される。そして、続くステップ A80で は、第二データ処理部 10においてこれらの特徴データ Dが解析され、このフローが
C
終了する。
[0065] 第二データ処理部 10における具体的な解析フローや判定及びその出力フローに ついては説明を省略するが、右足に係る時系列特徴データ D及び左足に係る時系 列特徴データ Dの各々に対して、歩行間隔時間のゆらぎのスケーリング指数 αが演
C
算され、 α = 1である状態を基準として、歩行状態が良好であるか否かが判定される 。また、このような判定結果はモニタ 15へ出力される。なお、特徴データ Dの解析が
C
一旦終了した時点で、ステップ Α20での判定に係る測定信号 Sの総数力 Sリセットされ [0066] [3.効果]
このように、本実施形態に係るデータ処理装置によれば、信号検出装置 1から入力 された測定信号 Sが二系統の信号過程の各々で処理される。一方は基本測定信号 生成部 2及び基本データ抽出部 3における抽出範囲 Αの設定のための信号処理で あり、他方は特徴データ抽出部 5における特徴データ Dの抽出のための信号処理で
c
ある。
[0067] 前者の信号処理過程においては測定信号 Sに前処理が施されている力 その結果 設定される抽出範囲 Aは、誤差を許容しうる幅を有するものであるため、前処理に伴 う誤差の影響を相殺することができる。また、後者の信号処理過程においては前処理 を施さずに直接特徴データ Dの抽出処理がなされているため、正確な特徴データ D
C
を取り出すこと力 sでさる。
C
[0068] つまり、一般に、前処理としての信号処理を行うと、測定信号の特徴が把握しやすく するなる反面、特徴的な情報の「抜け」が生じることになる力 本データ処理装置によ れば、前処理の結果を参照ながら、「抜け」のない元の測定信号の中から特徴データ Dを抽出するため、データ処理の信頼性を向上させることができる。
C
また、特徴データ抽出部 5における特徴データ Dの抽出に際し、抽出範囲 A内に
C
含まれる測定信号 Sのみが参照され、抽出範囲 A以外の測定信号 Sが除外されるた め、歩行状態を特徴付ける測定信号 Sとの相関が強いと考えられる部分の情報のみ を容易に取り出すことができ、十分なデータ処理速度を確保することが可能となる。
[0069] さらに、特徴データ抽出部 5は、抽出範囲 Aに含まれる測定信号 Sのうちの最大値 を特徴データ Dとして抽出するようになっている。この構成により、ノイズが多く特徴
C
の抽出が困難とされているような動物のバイタルサインの中からでも、正確に特徴デ 一タを由出すること力できる。 また、基本測定信号生成部 2における前処理としてのフィルタ処理は、ローパスフィ ルタ,位相フィルタ及びハイパスフィルタといった一般的な信号処理であり、実施が容 易であるとともに、短時間に処理を済ませることができる。一方、特徴データ抽出部 5 における演算処理に関しても、複雑な演算が不要であり素早く結果を得ることができ る。特に、本実施形態の第一データ処理部 9における処理内容は、線形解析の手法 による信号処理から構成されて!/、るため、例えば非線形解析の手法を用いた信号処 理と比較して構成が簡素であるという利点がある。
[0070] また、本データ処理装置は、単に基本測定信号 Sに基づいて基本データ Dを抽
B B
出するだけでなぐ抽出範囲設定部 4及び特徴データ抽出部 5を備えた構成となって いる。つまり、フィルタ処理によって生成される誤差を小さくするための従来の手法と して、位相フィルタにおける位相変化の遅延量を微調整するというものがある力 本 発明によれば、正確なデータが抽出されるが故にこのような微調整の必要がないうえ 、実際に前処理を施した後の基本データ Dの確認も不要である。したがって、測定
B
信号 Sの検出及び測定信号 Sの第一データ処理部 9への入力から第二データ処理 部 10における結果の出力に至るまでの全信号処理過程を完全に自動化することが 可能となり、信号処理の労力を格段に低減させることが可能となる。
[0071] さらに、本実施形態のデータ処理装置では、第一データ処理部 9におけるデータ処 理の後、特徴データ D群のピーク間隔時間の揺らぎを観察する非線形解析手法が
C
用いられている。前述の通り、仮に前処理でノイズを除去してしまえば却って非線形 構造を把握しにくくなりかねないが、本データ処理装置では第二データ処理部 10へ 入力される特徴データ D群は、ノイズが除去されていない生の情報から抽出された
c
ものであるため、特徴的な情報の抜けが生じず、正確に非線形構造を把握すること が可能となる。つまり、一見ノイズのように見える部分の情報を切り捨てることなく正確 な特徴部分の情報を抽出して非線形解析を実施することができる。
[0072] このように、本データ処理装置によれば、簡素な構成で、データ処理速度及びデー タ処理精度を向上させることができる。
[0073] [4.その他]
以上、本発明の実施形態について説明したが、本発明は、上記実施形態に限定さ れず、本発明の趣旨を逸脱しない範囲で種々変形することが可能である。
例えば、上述の実施形態では、第一データ処理部 9及び第二データ処理部 10に おけるデータ処理機能がプログラムとして構成されたものを例示した力 これらの機 能を実現手段はこれに限定されない。例えば、各第一データ処理部 9及び第二デー タ処理部 10を、 ROM, RAM, CPU等を内蔵したワンチップマイコンとして構成して もよいし、あるいは、デジタル回路やアナログ回路といった電子回路として形成しても よい。
[0074] なお、前述の通り、本発明のデータ処理装置においては、測定信号 Sの検出から 結果の出力に至るまでの信号処理過程を自動化することが可能なため、上述のよう な小型のマイコンで本発明に係るデータ処理装置を構成する場合、本発明に係るモ ユタ 15と同様の機能を備えた小型表示装置や、信号検出装置 1と同様の機能を備え たマイクロセンサを搭載させて、入出力一体型の小型処理装置を製造することも可能 である。
[0075] また、上述の実施形態では、信号検出手段として加速度信号を検出するための加 速度センサが適用されている力 S、本データ処理装置の演算対象となる信号としては、 種々の対象体の状態に関わる様々なパラメータが考えられる。
まず、人間や動物の状態に関係するバイタルサインを演算対象の信号とした場合、 歩行等の運動に伴う体動や呼吸数,心拍数,体温,皮膚表面温度,皮膚電位,脈波 (脈拍数),脳波,血流量,唾液などの体液成分,呼吸気中や血中の酸素飽和度, 血糖値,心電,電気伝導度,体重 (着座面への圧力),まばたきの数や周期,発汗量 ,その他身体から発せられる電磁波の強度や化学物質濃度等が挙げられる。
[0076] また、機械の作動状態に関係する物理量を演算対象の信号とした場合、その機械 の作動出力変動や仕事率の変動,その機械の作動によってなされた作業精度の変 動等を用いることが考えられる。さらに、天候,地震,火山活動といった自然現象を観 察する場合には、気圧や気温,風速,風向,地殻変動等を演算対象の信号とするこ とが考えられる。
なお、上述の実施形態では、第一データ処理部 9で演算処理される測定信号 Sとし て、その定義域が時間領域からなる時系列データを用いている力 これの代わりに、 その定義域が二次元空間領域 (及び時間領域)からなる画像信号を用いることも考え られる。
[0077] また、上述の実施形態では、信号検出装置 1で検出された測定信号 Sが第一デー タ処理部 9 直接入力される構成となっているが、信号検出装置 1と第一データ処理 部 9とを分離した構成としてもよい。例えば、信号検出装置 1で検出された測定信号 S の時系列データを何らかの記憶媒体に保存しておき、演算処理が必要となった時点 でそれらの時系列データを第一データ処理部 9 入力することが考えられる。この場 合それらの時系列データを測定信号記憶部 11 入力してもよいが、測定信号記憶 部 11を介さずに基本測定信号生成部 2及び特徴データ抽出部 5の各々へ入力して あよい。
[0078] また、上述の実施形態では、基本測定信号生成部 2においてローパスフィルタ,位 相フィルタ及びハイパスフィルタの三種のフィルタ処理が施されて!/、る力 バンドパス フィルタやノッチフィルタを併用してもよい。なお、基本測定信号生成部 2における前 処理とは、パラメータの変動を見つけやすくするための不可逆的な(非可逆変化を伴 う)演算処理全般のことを指している。つまり、パラメータの変動を見つけやすくするた めの演算処理であれば、具体的な処理内容がフィルタ処理でなくてもよい。例えば、 ヒルベルト変換処理,エンベロープ処理,フーリエ変換処理,加算平均の手法を用い た信号処理,ウェーブレット解析処理,フラクタル解析処理等を用いることが考えられ る。また、任意の信号加算や減算,比例処理,積分処理,微分処理等も含まれる。
[0079] なお、デジタル回路やアナログ回路といった電子回路を使って第一データ処理部 9 及び第二データ処理部 10を構成する場合には、上述の実施形態に記載されたよう なデジタルフィルタの代わりに、アナログフィルタを適用すればよい。すなわち、第一 データ処理部 9において実施されるデータ処理は、アナログ信号処理であってもよい また、上述の実施形態では、基本データ抽出部 3において、基本測定信号 Sの 0
B
点に対応する測定信号 Sが抽出されるようになっている力 このような抽出対象は、本 データ処理装置における演算対象に応じて適宜設定することができる。抽出範囲設 定部 4における抽出範囲 Aの位置や幅、及び、特徴データ抽出部 5において特徴デ ータ Dを取り出す位置についても同様である。
C
[0080] また、上述の実施形態では、解析部 7においてデトレンド変動解析の手法が用いら れているが、解析方法はこれに限定されない。なお、前処理が解析結果に与える一 般的な影響の大きさを考慮すると、解析部 7における解析手法が非線形解析手法で ある場合には、線形解析手法の場合と比較してより正確な解析結果が期待できるも のと考えられる。
産業上の利用可能性
[0081] 本発明のデータ処理装置,データ処理方法及びデータ処理プログラムの用途は特 に制限されず、機械,動植物や自然現象を測定対象として得られたデータに基づい て当該対象の状態を把握するためのデータ処理装置,データ処理方法及びデータ 処理プログラムとして好適に利用することが可能である。とりわけ、本発明は測定され たデータ中の非線形構造を抽出する用途に用いて有用である。
以上、本発明を特定の態様を用いて詳細に説明したが、本発明の意図と範囲を離 れることなく様々な変更が可能であることは当業者に明らかである。
なお、本出願は、 2006年 9月 19日付で出願された日本特許出願(特願 2006— 2
52356号)に基づいており、その全体が引用により援用される。

Claims

請求の範囲
[1] 対象体の状態に応じて変動するパラメータを測定信号として検出する測定信号検 出手段と、
該測定信号検出手段で検出された該測定信号に対し、該パラメータの変動を把握 するための前処理としての信号処理を施して、基本測定信号を生成する信号処理手 段と、
該信号処理手段で生成された該基本測定信号に基づ!/、て、該パラメータの変動を 特徴付ける測定信号を基本データとして抽出する基本データ抽出手段と、
該基本データ抽出手段で抽出された該基本データによって規定される所定の領域 を該測定信号の抽出範囲として設定する抽出範囲設定手段と、
該測定信号検出手段で検出された該測定信号のうち、該抽出範囲設定手段で設 定された該抽出範囲に含まれる測定信号の中から、該対象体の状態を特徴付ける 測定信号を特徴データとして抽出する特徴データ抽出手段と
を備えたことを特徴とする、データ処理装置。
[2] 該測定信号検出手段が、動物のバイタルサインを該測定信号として検出する ことを特徴とする、請求項 1記載のデータ処理装置。
[3] 該信号処理手段が、線形解析手法を用いて該測定信号に対し該信号処理を施す ことを特徴とする、請求項 1又は 2記載のデータ処理装置。
[4] 該信号処理手段が、該測定信号検出手段で検出された該測定信号から、予め設 定された所定の周波数成分を濾波するフィルタ処理手段を有してなる
ことを特徴とする、請求項;!〜 3の何れか 1項に記載のデータ処理装置。
[5] 該測定信号検出手段が、該測定信号として略周期的に変動するパラメータを検出 し、
該信号処理手段が、該測定信号を平滑化した波動を該基本測定信号として生成 するとともに、
該基本データ抽出手段が、該基本測定信号における変動ピークの検出時刻を該 基本データとして抽出する
ことを特徴とする、請求項;!〜 4の何れか 1項に記載のデータ処理装置。
[6] 該抽出範囲設定手段が、該ピークの検出時刻の近傍時刻を該抽出範囲として設 定する
ことを特徴とする、請求項 5記載のデータ処理装置。
[7] 該特徴データ抽出手段が、該抽出範囲に含まれる測定信号における変動ピークの 検出時刻を該特徴データとして抽出する
ことを特徴とする、請求項 6記載のデータ処理装置。
[8] 該特徴データ抽出手段で抽出された該特徴データに基づき、該測定信号中の非 線形構造を抽出し、該対象体の状態を解析する演算処理手段
をさらに備えたことを特徴とする、請求項 1〜7の何れ力、 1項に記載のデータ処理装 置。
[9] 対象体の状態に応じて変動するパラメータを測定信号として検出する測定信号検 該測定信号検出ステップで検出された該測定信号に対し、該パラメータの変動を 把握するための前処理としての信号処理を施して、基本測定信号を生成する信号処 該信号処理ステップで生成された該基本測定信号に基づ!/、て、該パラメータの変 動を特徴付ける測定信号を基本データとして抽出する基本データ抽出ステップと、 該基本データ抽出ステップで抽出された該基本データによって規定される所定領 域を抽出範囲として設定する抽出範囲設定ステップと、
該測定信号検出ステップで検出された該測定信号のうち、該抽出範囲設定ステツ プで設定された該抽出範囲に含まれる測定信号の中から、該対象体の状態を特徴 付ける測定信号を特徴データとして抽出する特徴データ抽出ステップと
を備えたことを特徴とする、データ処理方法。
[10] 該測定信号検出ステップにおいて、動物のバイタルサインを該測定信号として検出 することを特徴とする、請求項 9記載のデータ処理方法。
[11] 該信号処理ステップにお!/、て、線形解析手法を用いて該測定信号に対し該信号 処理を施す
ことを特徴とする、請求項 9又は 10記載のデータ処理方法。
[12] 該信号処理ステップにおいて、該測定信号検出ステップで検出された該測定信号 から、予め設定された所定の周波数成分を濾波する
ことを特徴とする、請求項 9〜; 11の何れか 1項に記載のデータ処理方法。
[13] 該測定信号処理ステップにおいて、該測定信号として略周期的に変動するパラメ ータを検出し、
該信号処理ステップにおいて、該測定信号を平滑化した波動を該基本測定信号と して生成するとともに、
該基本データ抽出ステップにおいて、該基本測定信号における変動ピークの検出 時刻を該基本データとして抽出する
ことを特徴とする、請求項 9〜; 12の何れか 1項に記載のデータ処理方法。
[14] 該抽出範囲設定ステップにおいて、該ピークの検出時刻の近傍時刻を該抽出範囲 として設定する
ことを特徴とする、請求項 13記載のデータ処理方法。
[15] 該特徴データ抽出ステップにおいて、該抽出範囲に含まれる測定信号における変 動ピークの検出時刻を該特徴データとして抽出する
ことを特徴とする、請求項 14記載のデータ処理方法。
[16] 該特徴データ抽出ステップで抽出された該特徴データに基づき、該測定信号中の 非線形構造を抽出し、該対象体の状態を解析する演算処理ステップ
をさらに備えたことを特徴とする、請求項 9〜; 15の何れ力、 1項に記載のデータ処理方 法。
[17] コンピュータを、信号処理手段,基本データ抽出手段,抽出範囲設定手段及び特 徴データ抽出手段として機能させるためのデータ処理プログラムであって、
該信号処理手段が、測定信号として検出された、対象体の状態に応じて変動する パラメータの変動を把握するための前処理としての信号処理を該測定信号に対して 施すとともに、基本測定信号を生成し、
該基本データ抽出手段が、該信号処理手段で生成された該基本測定信号に基づ いて、該パラメータの変動を特徴付ける測定信号を基本データとして抽出し、 該抽出範囲設定手段が、該基本データ抽出手段で抽出された該基本データによ つて規定される所定の領域を該測定信号の抽出範囲として設定し、
該特徴データ抽出手段が、該測定信号検出手段で検出された該測定信号のうち、 該抽出範囲設定手段で設定された該抽出範囲に含まれる測定信号の中から、該対 象体の状態を特徴付ける測定信号を特徴データとして抽出する
ことを特徴とする、データ処理プログラム。
PCT/JP2007/067832 2006-09-19 2007-09-13 Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données WO2008035611A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2006-252356 2006-09-19
JP2006252356A JP4952162B2 (ja) 2006-09-19 2006-09-19 データ処理装置,データ処理方法及びデータ処理プログラム

Publications (1)

Publication Number Publication Date
WO2008035611A1 true WO2008035611A1 (fr) 2008-03-27

Family

ID=39200438

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2007/067832 WO2008035611A1 (fr) 2006-09-19 2007-09-13 Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données

Country Status (2)

Country Link
JP (1) JP4952162B2 (ja)
WO (1) WO2008035611A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012036135A1 (ja) * 2010-09-17 2012-03-22 三菱化学株式会社 情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2012213624A (ja) * 2011-03-30 2012-11-08 Denso It Laboratory Inc 身体能力判定装置及び身体能力判定方法
CN111582675A (zh) * 2020-04-22 2020-08-25 北京启安智慧科技有限公司 Natech事件的关键特征分析系统及方法
WO2022160498A1 (zh) * 2021-01-29 2022-08-04 深圳市科曼医疗设备有限公司 基于血液细胞分析仪的自动检测方法和装置

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5522338B2 (ja) * 2008-10-28 2014-06-18 日本電気株式会社 状況判定装置、状況判定システム、その方法及びプログラム
JP5321002B2 (ja) * 2008-11-18 2013-10-23 オムロンヘルスケア株式会社 体動バランス検出装置、体動バランス検出プログラム、体動バランス検出方法
JP5812381B2 (ja) * 2010-11-25 2015-11-11 公立大学法人首都大学東京 振動体の異常検知方法および装置
JP6233837B2 (ja) * 2013-11-19 2017-11-22 公立大学法人首都大学東京 睡眠段階判定装置、睡眠段階判定プログラムおよび睡眠段階判定方法
JP6455852B2 (ja) * 2015-02-09 2019-01-23 公立大学法人首都大学東京 信号解析システム、方法およびプログラム
JP6778369B2 (ja) * 2015-03-10 2020-11-04 東京都公立大学法人 生体信号解析システムおよびプログラム

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10260066A (ja) * 1997-01-20 1998-09-29 Akiyuuto Kk 波形検出装置およびその装置を利用した状態監視システム
JP2000166877A (ja) * 1998-12-09 2000-06-20 Mitsubishi Chemicals Corp 生体リズム検査装置及び生体リズム検査方法
JP2001299766A (ja) * 2000-04-26 2001-10-30 Sharp Corp 健康状態診断方法およびその方法を用いた健康状態診断装置
JP2007244495A (ja) * 2006-03-14 2007-09-27 Sony Corp 体動検出装置、体動検出方法および体動検出プログラム

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3564501B2 (ja) * 2001-03-22 2004-09-15 学校法人明治大学 乳幼児の音声解析システム
JP2003331265A (ja) * 2002-05-14 2003-11-21 Matsushita Electric Ind Co Ltd 目画像撮像装置及び虹彩認証装置
JP3935080B2 (ja) * 2003-01-10 2007-06-20 株式会社Ctiサイエンスシステム 伝播音の帯域別音圧解析装置、周波数解析装置及び周波数解析方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10260066A (ja) * 1997-01-20 1998-09-29 Akiyuuto Kk 波形検出装置およびその装置を利用した状態監視システム
JP2000166877A (ja) * 1998-12-09 2000-06-20 Mitsubishi Chemicals Corp 生体リズム検査装置及び生体リズム検査方法
JP2001299766A (ja) * 2000-04-26 2001-10-30 Sharp Corp 健康状態診断方法およびその方法を用いた健康状態診断装置
JP2007244495A (ja) * 2006-03-14 2007-09-27 Sony Corp 体動検出装置、体動検出方法および体動検出プログラム

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012036135A1 (ja) * 2010-09-17 2012-03-22 三菱化学株式会社 情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体
JPWO2012036135A1 (ja) * 2010-09-17 2014-02-03 三菱化学株式会社 情報処理方法、情報処理装置、出力装置、情報処理システム、情報処理用プログラムおよび同プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2012213624A (ja) * 2011-03-30 2012-11-08 Denso It Laboratory Inc 身体能力判定装置及び身体能力判定方法
CN111582675A (zh) * 2020-04-22 2020-08-25 北京启安智慧科技有限公司 Natech事件的关键特征分析系统及方法
CN111582675B (zh) * 2020-04-22 2024-02-20 北京启安智慧科技有限公司 Natech事件的关键特征分析系统及方法
WO2022160498A1 (zh) * 2021-01-29 2022-08-04 深圳市科曼医疗设备有限公司 基于血液细胞分析仪的自动检测方法和装置

Also Published As

Publication number Publication date
JP4952162B2 (ja) 2012-06-13
JP2008073077A (ja) 2008-04-03

Similar Documents

Publication Publication Date Title
WO2008035611A1 (fr) Dispositif de traitement de données, procédé de traitement de données et programme de traitement de données
US7892189B2 (en) Movement analysis display apparatus and movement analyzing method
WO2012051300A2 (en) Method and system for detection and rejection of motion/noise artifacts in physiological measurements
JP6843855B2 (ja) 動的生体信号解析において適応的にノイズを定量化するためのシステム及び方法
Drapała et al. Two stage EMG onset detection method
US20110152704A1 (en) Sleep determination device and sleep determination method
JP5544365B2 (ja) マルチパラメーターモニタリングにおける改善又はマルチパラメーターモニタリングに関する改善
EP3603514A1 (en) Emotion estimating apparatus
CN111860102A (zh) 用于在含噪环境中分析系统的状态的设备和方法
CN107320097A (zh) 利用肌电信号边际谱熵提取肌肉疲劳特征的方法和装置
JP5046286B2 (ja) ストレス評価装置、ストレス評価システムおよびストレス評価プログラム
JP2008188092A (ja) データ処理方法,データ処理装置及びデータ処理プログラム
KR101744691B1 (ko) 심탄도 신호를 이용한 심박 검출 방법 및 그 장치
Shi et al. A Novel Method of sEMG Signal Segmentation
WO2004098409A1 (ja) 心拍や呼吸等の生体信号の抽出法及び装置
CN104305958B (zh) 一种极短时自主神经状态的光电容积波多参量分析方法
CN113616194B (zh) 一种监测手部震颤频率和强度的装置及方法
JPWO2019130840A1 (ja) 信号処理装置、解析システム、信号処理方法および信号処理プログラム
Hu et al. Multiscale entropy: recent advances
JP2010131061A (ja) 眠気検出装置
JP4725218B2 (ja) 脳機能測定装置
CN113876316A (zh) 下肢屈伸活动异常检测系统、方法、装置、设备及介质
CN105395191A (zh) 基于年龄段检测的人员生理状态检测装置
JP7001806B2 (ja) 呼吸音解析装置、呼吸音解析方法、コンピュータプログラム及び記録媒体
JP5419157B2 (ja) 脳波測定装置および脳波測定方法

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: 07807239

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: 07807239

Country of ref document: EP

Kind code of ref document: A1