WO2022149215A1 - 信号解析装置、信号解析方法及びプログラム - Google Patents

信号解析装置、信号解析方法及びプログラム Download PDF

Info

Publication number
WO2022149215A1
WO2022149215A1 PCT/JP2021/000209 JP2021000209W WO2022149215A1 WO 2022149215 A1 WO2022149215 A1 WO 2022149215A1 JP 2021000209 W JP2021000209 W JP 2021000209W WO 2022149215 A1 WO2022149215 A1 WO 2022149215A1
Authority
WO
WIPO (PCT)
Prior art keywords
electrocardiogram
information
cumulative distribution
signal analysis
function
Prior art date
Application number
PCT/JP2021/000209
Other languages
English (en)
French (fr)
Inventor
信吾 塚田
Original Assignee
日本電信電話株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 日本電信電話株式会社 filed Critical 日本電信電話株式会社
Priority to JP2022573841A priority Critical patent/JPWO2022149215A1/ja
Priority to US18/267,688 priority patent/US20240050018A1/en
Priority to PCT/JP2021/000209 priority patent/WO2022149215A1/ja
Priority to EP21917443.0A priority patent/EP4275606A1/en
Priority to PCT/JP2021/033138 priority patent/WO2022149307A1/ja
Priority to US18/039,553 priority patent/US20240016433A1/en
Priority to EP21917602.1A priority patent/EP4238498A1/en
Priority to JP2022573946A priority patent/JP7457267B2/ja
Priority to EP21917603.9A priority patent/EP4241691A1/en
Priority to US18/265,663 priority patent/US20240057924A1/en
Priority to JP2022573945A priority patent/JP7457266B2/ja
Priority to PCT/JP2021/044316 priority patent/WO2022149381A1/ja
Priority to PCT/JP2021/044327 priority patent/WO2022149382A1/ja
Publication of WO2022149215A1 publication Critical patent/WO2022149215A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/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
    • 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/35—Detecting specific parameters of the electrocardiograph cycle by template matching
    • 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 pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • 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 pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Definitions

  • the present invention relates to a signal analysis device, a signal analysis method and a program.
  • the electrocardiogram is useful information for grasping the state of the heart. For example, if the electrocardiogram is used, it can be determined whether or not the subject is in a state in which there is a high possibility of heart failure (Non-Patent Document 1).
  • the information obtained from the ECG waveform may not always be sufficient to grasp the state of the heart.
  • the way of developing a disease related to the heart may be different.
  • the information obtained from the waveform of the electrocardiogram may not always be sufficient for grasping the state of the heart, it may be difficult to obtain other information depending on the disease.
  • the onset can be suppressed by observing the state of the heart from the waveform of an electrocardiogram in daily life.
  • Such a situation is not limited to the case of grasping the state of the heart based on the waveform of the electrocardiogram.
  • Such a situation is also common when grasping the state of the heart based only on the time-series biological information regarding the heartbeat.
  • the time-series biological information regarding the heartbeat shows, for example, a waveform showing a change in the electrocardiographic potential, a waveform showing a change in the image of the heart, a waveform showing a change in the pressure of the heart, and a change in blood flow. It is a waveform or a waveform showing a change in heartbeat.
  • the waveform of the electrocardiogram is also an example of time-series biological information regarding the heartbeat.
  • the image of the heart is, for example, an MRI (Magnetic Resonance Image) image of the heart or a CT (Computed Tomography) image of the heart.
  • One aspect of the present invention shows a biological information acquisition unit that acquires time-series biological information regarding the beat of the heart to be analyzed, and a candidate for a channel activity timing distribution that is a distribution of the timing of ion channel activity in the heart. It is a signal analysis device including an analysis unit that acquires information indicating the state of activity of the ion channel based on the biological information and the distribution candidate information, using the information as distribution candidate information.
  • the flowchart which shows an example of the flow of the process executed by the signal analysis apparatus 1 in embodiment.
  • FIG. 3 is a third diagram showing an example of an analysis result of the signal analysis device 1 in the embodiment.
  • FIG. 4 is a fourth diagram showing an example of the analysis result of the signal analysis device 1 in the embodiment.
  • FIG. 1 is a diagram showing an example of a hardware configuration of the signal analysis device 1 of the embodiment.
  • the signal analysis device 1 will be described by taking at least a case where analysis is performed based on the waveform of an electrocardiogram as an example.
  • the signal analysis device 1 is not limited to the electrocardiogram, and can perform the same analysis based on time-series biological information regarding the heartbeat.
  • the time-series biological information related to the heartbeat is, for example, a waveform showing a change in cardiac potential, a waveform showing a change in a heart image, a waveform showing a change in heart pressure, or a waveform showing a change in blood flow.
  • the signal analysis device 1 may use a waveform indicating the cardiac cycle instead of the waveform of the electrocardiogram.
  • the signal analyzer 1 may use a waveform showing a change in the image of the heart instead of the waveform of the electrocardiogram.
  • the signal analyzer 1 may use a waveform showing a change in heart pressure instead of the waveform of the electrocardiogram.
  • the signal analyzer 1 may use a waveform showing a change in blood flow rate instead of the waveform of the electrocardiogram.
  • the signal analysis device 1 may use a waveform showing a change in heart sounds instead of the waveform of the electrocardiogram.
  • the electrocardiogram is also an example of time-series biological information regarding the heartbeat.
  • the time-series biological information regarding the heartbeat may be the time-series biological information regarding the periodic heartbeat.
  • the image of the heart is, for example, an MRI (Magnetic Resonance Image) image of the heart or a CT (Computed Tomography) image of the heart.
  • the signal analysis device 1 acquires the waveform of the electrocardiogram of the heart to be analyzed (hereinafter referred to as "target heart”).
  • the signal analyzer 1 is based on the acquired ECG waveform, and is referred to as information showing the state of ion channel activity between the myocardium of the outer layer and the myocardium of the inner layer of the target heart (hereinafter referred to as "ion channel activity information”. ).
  • an electromotive force dipole model (Reference 1) that explains the relationship between myocardial movement and electrocardiogram is known.
  • the cardiac electromotive force dipole model the myocardium is modeled with two layers, an outer layer (hereinafter referred to as “outer myocardial layer”) and an inner layer (hereinafter referred to as “inner myocardial layer”).
  • the outer and inner myocardial layers are modeled as different sources of electromotive force.
  • the synthetic wave of the epicardial side myocardial action potential and the endocardial side myocardial action potential substantially coincides with the time change of the body surface potential observed on the body surface.
  • the graph showing the time change of the electric potential on the body surface is the waveform of the electrocardiogram.
  • the epicardial myocardial action potential is the result of directly measuring the change in electromotive force caused by the pulsation of the outer myocardium by inserting a catheter.
  • the endocardial myocardial action potential is the result of directly measuring the change in electromotive force caused by the pulsation of the myocardial inner layer by inserting a catheter.
  • the outer layer of the myocardium in the electromotive force dipole model is a collection of cells. Therefore, the timing of beats of cells in the outer layer of myocardium in one beat of the outer layer of myocardium is not necessarily the same in all cells, and there may be a distribution in the timing of beats.
  • This also applies to the inner layer of the myocardium. That is, therefore, the timing of the beats of the cells in the inner layer of the myocardium in one beat of the inner layer of the myocardium is not necessarily the same in all the cells, and there may be a distribution in the timing of the beats.
  • the possibility that such a distribution exists in the timing of cell beats is not assumed in the electromotive force dipole model.
  • multiple types of ion channels act at the timing corresponding to the beat timing of each cell for each cell.
  • the timing of the activity of the ion channel is the timing according to the type of the ion channel and is not necessarily the same.
  • sodium channels and potassium channels have different timings of activity. Therefore, if there is a distribution in the timing of cell beats, the timing of ion channel activity in one myocardial beat also has a distribution according to the distribution of cell beat timings.
  • the timing of ion channel activity is the timing according to the type of ion channel, if there is a distribution in the timing of cell beats, the channel activity timing distribution will be different for ion channels that are active at different timings. Appears at the timing.
  • the channel activity timing distribution is the distribution of ion channel activity timing.
  • the channel activity timing distribution is, for example, a Gaussian distribution.
  • the channel activity timing distribution may satisfy the condition that the cumulative distribution function of the probability density function representing the channel activity timing distribution is a function representing the sigmoid curve.
  • the function representing the sigmoid curve may be, for example, a sigmoid function, a cumulative normal distribution function, a Gompertz function, or a Gudermannian function.
  • the channel activity timing distribution may be, for example, a logistic distribution.
  • the signal analysis device 1 acquires ion channel activity information based on the acquired waveform of the electrocardiogram of the target heart using information indicating candidates for the channel activity timing distribution (hereinafter referred to as “distribution candidate information”).
  • the distribution candidate information is represented by, for example, a cumulative distribution function (hereinafter referred to as “channel cumulative distribution function”) in which a function representing a candidate for a channel activity timing distribution is used as a probability density function.
  • the function representing the candidate of the channel activity timing distribution is a function in which the value representing the shape of the channel activity timing distribution such as the variance and the average value of the channel activity timing distribution is expressed by a parameter.
  • the parameter representing the shape of the channel activity timing distribution is referred to as a shape parameter. Therefore, the cumulative distribution function that expresses the distribution candidate information is a functional having one or a plurality of shape parameters.
  • ion channel activity information acquisition process the process of acquiring ion channel activity information based on the acquired ECG waveform of the target heart using the distribution candidate information is referred to as ion channel activity information acquisition process.
  • the acquired ion channel activity information is an example of the analysis result of the waveform of the electrocardiogram of the target heart by the signal analysis device 1.
  • the signal analysis device 1 will be described by taking the case where the distribution candidate information is the channel cumulative distribution function as an example.
  • the ion channel activity information acquisition process is performed by fitting the waveform of the electrocardiogram of the target heart using one or a plurality of channel cumulative distribution functions, and the shape parameter possessed by the channel cumulative distribution function. It is a process to determine the value of.
  • the channel cumulative distribution function is represented by, for example, the following equation (1) or equation (2).
  • Equation (1) is a normal distribution (Gaussian distribution).
  • Equation (2) is the cumulative distribution function of Equation (1).
  • "erf” is an error function (sigmoid function).
  • the mean is expressed as ⁇ a and ⁇ b , respectively, and the variance is expressed as ⁇ 2 a and ⁇ 2 b , respectively.
  • the difference is expressed by the following equation (3).
  • the first function described in the image G1 of FIG. 3 described later is represented by the following equation (4).
  • the average and variance of the cumulative distribution function of the first function and the second function can be obtained.
  • the process of approximating the third function to the R wave or T wave of the electrocardiographic potential is a process of minimizing the difference between the electrocardiographic potential and the third function by, for example, the least squares method.
  • the domain of the channel cumulative distribution function used for fitting does not necessarily have to be the same as the domain of the function representing the waveform of the electrocardiogram to be fitted. Also, the domains of the channel cumulative distribution function used for fitting do not necessarily have to be the same. Not necessarily the same means that at least one may be different or all may be the same.
  • the fitting in the ion channel activity information acquisition process is a process of optimizing the difference between the channel cumulative distribution function and the ECG waveform within the definition range of the channel cumulative distribution function, specifically using the shape parameter as a variable. ..
  • the channel cumulative distribution function is the result of the integral of the probability density function and has a one-to-one correspondence with the probability density function.
  • the probability density function is a function representing the channel activity timing distribution. Therefore, when the value of the shape parameter is fixed, the channel cumulative distribution function is a function representing the channel activity timing distribution. Therefore, when the value of the shape parameter is not fixed, the channel cumulative distribution function is a function representing candidates for the channel activity timing distribution as many as the number of combinations of the value of the shape parameter.
  • fitting is a process of optimizing using the value of a parameter representing the shape of a function as a variable. Therefore, executing fitting using the channel cumulative distribution function for fitting is a process of determining the channel activity timing distribution that minimizes the difference from the fitting target from the channel activity timing distribution candidates represented by the channel cumulative distribution function. be. Therefore, performing a fitting using the channel cumulative distribution function for the fitting is equivalent to acquiring information indicating how the activity of the ion channel is related to the fitting target. Therefore, when the fitting target is an electrocardiogram waveform, the fitting acquires information indicating how the activity of the ion channel is related to the electrocardiogram waveform.
  • the relationship between the activity of ion channels and the waveform of the electrocardiogram is shown, for example, by the value of the shape parameter of the fitting result. Therefore, the value of the shape parameter of the fitting result is an example of ion channel activity information.
  • each of the plurality of channel cumulative distribution functions used for fitting in such cases represents different types of channel activity timing distributions.
  • information showing how the activities of a plurality of types of ion channels are related to the waveform of the electrocardiogram is acquired.
  • the difference in the type of channel activity timing distribution means that one or both of the type of ion channel and the difference in the operation of opening and closing the ion channel are different.
  • the domain of the channel cumulative distribution function used for fitting does not necessarily have to be the same as the domain of the function representing the waveform of the electrocardiogram to be fitted. Further, when a plurality of channel cumulative distribution functions are used for fitting, each domain of each channel cumulative distribution function does not necessarily have to be the same. Not necessarily the same means that at least one may be different or all may be the same.
  • FIG. 2 is a diagram showing an example of the result of fitting the waveform of the electrocardiogram of the target heart in the embodiment by the four channel cumulative distribution functions.
  • the horizontal axis of FIG. 2 represents time, and the vertical axis represents potential. Both the horizontal axis and the vertical axis are arbitrary units (arbitrary units).
  • the four channel cumulative distribution functions in FIG. 2 are a first cumulative distribution function, a second cumulative distribution function, a third cumulative distribution function, and a fourth cumulative distribution function. At least one domain of the first cumulative distribution function, the second cumulative distribution function, the third cumulative distribution function, and the fourth cumulative distribution function in FIG. 2 is different. Specifically, the domain of the first cumulative distribution function and the domain of the second cumulative distribution function are the same, and are from time 0 to time T1. The domain of the third cumulative distribution function is from time T3 to time T5. The domain of the fourth cumulative distribution function is from time T2 to time T4. The domain of the waveform of the electrocardiogram to be fitted includes at least the period from time 0 to time T5.
  • the first cumulative distribution function is a candidate for the channel activity timing distribution of the sodium channel existing in the inner layer of the myocardium and represents a candidate for the channel activity timing distribution at the timing of the start of myocardial excitement.
  • the second cumulative distribution function is a candidate for the channel activity timing distribution of the potassium channel existing in the outer layer of the myocardium and represents a candidate for the channel activity timing distribution at the timing of the start of myocardial excitement.
  • the third cumulative distribution function is a candidate for the channel activity timing distribution of the sodium channel existing in the inner layer of the myocardium, and represents a candidate for the channel activity timing distribution at the timing of the start of suppression of myocardial excitement.
  • the fourth cumulative distribution function is a candidate for the channel activity timing distribution of potassium channels existing in the outer layer of the myocardium, and represents a candidate for the channel activity timing distribution at the start timing of suppression of myocardial excitement.
  • the "first fitting result” shows the result of fitting the ECG waveform by the first cumulative distribution function.
  • the fitting to the waveform of the electrocardiogram by the first cumulative distribution function is referred to as the first fitting.
  • the “second fitting result” shows the result of fitting the ECG waveform by the second cumulative distribution function.
  • the fitting to the waveform of the electrocardiogram by the second cumulative distribution function is referred to as the second fitting.
  • the “third fitting result” shows the result of fitting the ECG waveform by the third cumulative distribution function.
  • the fitting to the waveform of the electrocardiogram by the third cumulative distribution function is referred to as the third fitting.
  • the "fourth fitting result” shows the result of fitting the ECG waveform by the fourth cumulative distribution function.
  • the fitting to the waveform of the electrocardiogram by the fourth cumulative distribution function is referred to as the fourth fitting.
  • the “body surface potential” represents the waveform of the electrocardiogram to be fitted.
  • the first cumulative distribution function, the second cumulative distribution function, the third cumulative distribution function, and the fourth cumulative distribution function are all cumulative Gaussian distributions.
  • the "first fitting result” and the “second fitting result” are the results of fitting to the R wave.
  • the “third fitting result” and the “fourth fitting result” are the results of fitting to the T wave.
  • the period from time T1 to time T2 does not belong to the domain of any channel cumulative distribution function from the first cumulative distribution function to the fourth cumulative distribution function. Therefore, among the waveforms of the electrocardiogram shown in FIG. 2, the waveform in the period from time T1 to time T2 is not fitted.
  • the period during which fitting is not performed is represented by a constant function.
  • the period during which fitting is not performed is not limited to a constant function, and may be expressed by any predetermined function such as a linear function.
  • the result obtained by subtracting the "second fitting result” from the “first fitting result” shows substantially the same waveform as the waveform during the rising period of the R wave.
  • the result obtained by subtracting the "fourth fitting result” from the “third fitting result” shows substantially the same waveform as the waveform during the rising period of the T wave.
  • the fitting result for the R wave in FIG. 2 is under the condition that the result of subtracting the "second fitting result” from the “first fitting result” shows substantially the same waveform as the waveform during the rising period of the R wave. It may be the result of fitting.
  • the fitting result for the T wave in FIG. 2 is under the condition that the result of subtracting the "fourth fitting result” from the “third fitting result” shows substantially the same waveform as the waveform during the rising period of the T wave. It may be the result of fitting.
  • the difference between the functions may be fitted to the ECG waveform with respect to the overlapping domains. That is, for two channel cumulative distribution functions having an overlap of domains, fitting may be performed to minimize the difference between the functions and the waveform of the electrocardiogram.
  • FIG. 3 is an explanatory diagram illustrating that the difference between the two channel cumulative distribution functions in the embodiment can be fitted to a waveform that is substantially the same as the waveform at the rising edge of the R wave.
  • FIG. 3 shows four images, an image G1, an image G2, an image G3, and an image G4.
  • Each image of the image G1 to the image G4 shows a graph in which the horizontal axis is the time and the vertical axis is the potential.
  • the units of the horizontal axis and the vertical axis of each of the images G1 to G4 in FIG. 3 are arbitrary units.
  • the "first function” in FIG. 3 is an example of a cumulative distribution function.
  • the “second function” in FIG. 3 is an example of a cumulative distribution function and is a cumulative distribution function different from the “second function”.
  • the “third function” in FIG. 3 represents a function obtained by subtracting the "second function" from the "first function”.
  • the “third function” in FIG. 3 indicates that the shape is substantially the same as the shape of the R wave.
  • the "fourth function” in FIG. 3 is an example of a cumulative distribution function.
  • the “fifth function” in FIG. 3 is an example of a cumulative distribution function and is a cumulative distribution function different from the “fourth function”.
  • the “sixth function” in FIG. 3 represents a function obtained by subtracting the "fifth function” from the "fourth function”.
  • the “sixth function” in FIG. 3 shows that the shape is substantially the same as the shape of the R wave.
  • the "seventh function” in FIG. 3 is an example of a cumulative distribution function.
  • the "eighth function” in FIG. 3 is an example of a cumulative distribution function and is a cumulative distribution function different from the "seventh function”.
  • the "9th function” in FIG. 3 represents a function obtained by subtracting the "7th function” from the "8th function”.
  • the “9th function” in FIG. 3 indicates that the shape is substantially the same as the shape of the R wave.
  • the "10th function” in FIG. 3 is an example of a cumulative distribution function.
  • the “11th function” in FIG. 3 is an example of a cumulative distribution function and is a cumulative distribution function different from the “10th function”.
  • the "12th function” in FIG. 3 represents a function obtained by subtracting the "11th function” from the "10th function”.
  • the "12th function” in FIG. 3 shows that the shape is substantially the same as the shape of the R wave.
  • the function of the difference between the two cumulative distribution functions can express a shape that is substantially the same as the shape of the R wave.
  • This also applies to the shape of the T wave. That is, the function of the difference between the two cumulative distribution functions can express substantially the same shape as the shape of the T wave.
  • the function of the difference between the two cumulative distribution functions can express different waves of the width in the vertical axis direction and the width in the horizontal axis direction, such as "third function", "sixth function", and "ninth function". be.
  • the function of the difference between the two cumulative distribution functions can express a negative wave like the "12th function".
  • the signal analysis device 1 fits the waveform of the electrocardiogram of the target heart by the channel cumulative distribution function. Then, the signal analysis device 1 acquires ion channel activity information based on the fitting result.
  • the signal analysis device 1 includes a control unit 11 including a processor 91 such as a CPU connected by a bus and a memory 92, and executes a program.
  • the signal analysis device 1 functions as a device including a control unit 11, an input unit 12, a communication unit 13, a storage unit 14, and an output unit 15 by executing a program.
  • the processor 91 reads out the program stored in the storage unit 14, and stores the read program in the memory 92.
  • the signal analysis device 1 functions as a device including a control unit 11, an input unit 12, a communication unit 13, a storage unit 14, and an output unit 15.
  • the control unit 11 controls the operation of various functional units included in the signal analysis device 1.
  • the control unit 11 executes, for example, an ion channel activity information acquisition process.
  • the control unit 11 controls, for example, the operation of the output unit 15, and causes the output unit 15 to output the acquisition result of the ion channel activity information acquisition process.
  • the control unit 11 records, for example, various information generated by executing the ion channel activity information acquisition process in the storage unit 14.
  • the input unit 12 includes an input device such as a mouse, a keyboard, and a touch panel.
  • the input unit 12 may be configured as an interface for connecting these input devices to the signal analysis device 1.
  • the input unit 12 receives input of various information to the signal analysis device 1.
  • Information indicating, for example, the number of channel cumulative distribution functions used for fitting (hereinafter referred to as “number designation information”) is input to the input unit 12.
  • number designation information information indicating the number of channel cumulative distribution functions used for fitting
  • distribution shape designation information information indicating the shape of the distribution represented by each channel cumulative distribution function
  • domain designation information information indicating a domain for each channel cumulative distribution function used for fitting is input to the input unit 12.
  • the storage unit 14 may have stored a part of the number designation information, the distribution shape designation information, or the domain designation information in advance. In such a case, it is not necessary to input from the input unit 12 about the information stored in the storage unit 14 and a part of the number designation information, the distribution shape designation information, or the definition area designation information.
  • the signal analysis device 1 will be described by taking as an example a case where the storage unit 14 has stored the number designation information, the distribution shape designation information, or the domain designation information in advance.
  • the communication unit 13 includes a communication interface for connecting the signal analysis device 1 to an external device.
  • the communication unit 13 communicates with an external device via wired or wireless.
  • the external device is, for example, a device from which the waveform of the electrocardiogram of the target heart is transmitted.
  • the device from which the ECG waveform of the target heart is transmitted is, for example, an electrocardiogram measuring device.
  • the communication unit 13 acquires an electrocardiogram waveform from the electrocardiogram measuring device by communication.
  • the waveform of the electrocardiogram may be input to the input unit 12.
  • the storage unit 14 is configured by using a non-temporary computer-readable storage medium device such as a magnetic hard disk device or a semiconductor storage device.
  • the storage unit 14 stores various information about the signal analysis device 1.
  • the storage unit 14 stores information input via, for example, the input unit 12 or the communication unit 13.
  • the storage unit 14 stores an electrocardiogram input via, for example, the input unit 12 or the communication unit 13.
  • the storage unit 14 stores various information generated by, for example, executing an ion channel activity information acquisition process.
  • the output unit 15 outputs various information.
  • the output unit 15 includes, for example, a display device such as a CRT (Cathode Ray Tube) display, a liquid crystal display, and an organic EL (Electro-Luminence) display.
  • the output unit 15 may be configured as an interface for connecting these display devices to the signal analysis device 1.
  • the output unit 15 outputs, for example, the information input to the input unit 12.
  • the output unit 15 may display, for example, an electrocardiogram input to the input unit 12 or the communication unit 13.
  • the output unit 15 may display, for example, the execution result of the ion channel activity information acquisition process.
  • FIG. 4 is a diagram showing an example of the functional configuration of the control unit 11 in the embodiment.
  • the control unit 11 includes an electrocardiogram acquisition unit 110, a fitting information acquisition unit 120, an analysis unit 130, and a recording unit 140.
  • the electrocardiogram acquisition unit 110 acquires the waveform of the electrocardiogram of the target heart input to the input unit 12 or the communication unit 13.
  • the fitting information acquisition unit 120 acquires number designation information, distribution shape designation information, and domain designation information. When the number designation information, the distribution shape designation information and the definition area designation information are stored in the storage unit 14, the fitting information acquisition unit 120 reads the number designation information, the distribution shape designation information and the definition area designation information from the storage unit 14. ..
  • the analysis unit 130 includes a fitting unit 131 and an ion channel activity information acquisition unit 132.
  • the fitting unit 131 performs fitting on the waveform of the electrocardiogram acquired by the electrocardiogram acquisition unit 110 by using the channel cumulative distribution function indicated by the number designation information, the distribution shape designation information, and the domain designation information.
  • the ion channel activity information acquisition unit 132 acquires ion channel activity information based on the fitting result by the fitting unit 131.
  • the ion channel activity information acquisition unit 132 acquires, for example, the value of the shape parameter for each channel cumulative distribution function.
  • the value of the shape parameter acquired by the ion channel activity information acquisition unit 132 is an example of the ion channel activity information.
  • the analysis unit 130 acquires ion channel activity information based on the waveform of the electrocardiogram of the target heart and the distribution candidate information.
  • the storage unit 14 records various information generated by the processing executed by the control unit 11 in the storage unit 14.
  • FIG. 5 is a flowchart showing an example of the flow of processing executed by the signal analysis device 1 in the embodiment.
  • the electrocardiogram acquisition unit 110 acquires an electrocardiogram of the target heart via the input unit 12 or the communication unit 13 (step S101).
  • the fitting information acquisition unit 120 acquires the number designation information, the distribution shape designation information, and the domain designation information (step S102).
  • the fitting unit 131 performs fitting on the waveform of the electrocardiogram acquired in step S101 by using the channel cumulative distribution function indicated by the number designation information, the distribution shape designation information, and the domain designation information (step S103).
  • the ion channel activity information acquisition unit 132 acquires ion channel activity information based on the fitting result (step S104).
  • the acquired ion channel activity information is output to the output unit 15 (step S105).
  • step S105 a graph of the fitting result of each channel cumulative distribution function may be displayed. Further, the process of step S102 may be executed before the execution of the process of step S103, or may be executed before the execution of step S101.
  • the processing of step S103 and step S104 is an example of the processing executed by the analysis unit 130.
  • FIG. 6 is a first diagram showing an example of the analysis result of the signal analysis device 1 in the embodiment. More specifically, FIG. 6 is an example of the analysis result by the signal analysis device 1 for the waveform of the electrocardiogram of the target heart whose operation is normal.
  • the horizontal axis of FIG. 6 represents time, and the vertical axis represents potential.
  • the unit on the vertical axis is an arbitrary unit.
  • FIG. 6 shows the results of performing the first fitting and the second fitting on the R waveform of the depolarization of the electrocardiogram of the target heart whose movement is normal. Further, FIG. 6 shows an example of the results of performing the fourth fitting and the third fitting on the T waveform of the repolarized phase of the electrocardiogram of the target heart whose operation is normal.
  • FIG. 6 shows that in normal cardiac depolarization, the inner myocardium begins to activate ion channels faster than the outer myocardium and progresses rapidly.
  • FIG. 6 shows that the outer myocardium begins its activity slightly later than the timing at which the activity of the ion channels in the inner myocardium begins.
  • FIG. 6 shows that the difference between the timing at which the activity of the ion channel in the inner layer of the myocardium starts and the timing at which the activity of the outer layer of the myocardium starts shows a positive and sharp R waveform.
  • the mean value and variance of the inner layer side cumulative distribution function which is the first fitting result and the outer layer side cumulative distribution function which is the second fitting result represent the timing of the activity and the progress of the activity, respectively.
  • FIG. 6 shows that in the repolarized phase of the normal heart, the outer myocardium begins to inactivate the activity of ion channels earlier than the inner myocardium.
  • FIG. 6 shows that the inactivation of the myocardial inner layer lags behind it, and both proceed slowly. The difference is positive and shows a gradual T waveform.
  • the mean value and variance of the inner layer side cumulative distribution function, which is the fourth fitting result, and the outer layer side cumulative distribution function, which is the third fitting result, represent the timing of inactivation and the progress of inactivation, respectively.
  • the averaging and dispersion of the inner layer side cumulative distribution function and the outer layer side cumulative distribution function obtained by the first to fourth fittings are the collective activation of the depolarized ion channels and the ions of the repolarized phase. Represents the timing and progress of collective inactivation of channels.
  • FIG. 6 shows that ion channel activity information can be acquired only from the electrocardiogram without inserting a catheter by the signal analyzer 1.
  • FIG. 7 is a second diagram showing an example of the analysis result of the signal analysis device 1 in the embodiment. More specifically, FIG. 7 is an example of the analysis result by the signal analysis device 1 for the waveform of the electrocardiogram of the target heart whose operation is normal.
  • FIG. 7 shows three results of graph G5, graph G6 and result G7.
  • the “inner layer side cumulative distribution function” shows the fitting result of a function representing a candidate for a collective channel activity timing distribution of ion channels existing in the inner layer of the myocardium.
  • the “outer layer side cumulative distribution function” shows the fitting result of a function representing a candidate for a collective channel activity timing distribution of ion channels existing in the outer layer of the myocardium.
  • the “body surface potential” is a function representing the time change of the body surface potential, and is a waveform of an electrocardiogram. The horizontal axis of FIG. 7 represents time, and the vertical axis represents potential.
  • Both the horizontal axis and the vertical axis are arbitrary units (arbitrary units).
  • the length of time represented by the interval of one scale on each horizontal axis of FIGS. 7 to 9 is the same. Further, in each of the vertical axes of FIGS. 7 to 9, 1 represents the maximum value of the cumulative Gaussian distribution.
  • Graph G5 represents all of the ECG waveforms generated by one beat.
  • Graph G6 is a part of graph G5 and shows an enlarged view of a T wave region.
  • the region of the T wave is the region shown as region A1 in FIG.
  • Result G7 shows a statistic of two Gaussian distributions.
  • Each value of the result G7 represents a statistic of two Gaussian distributions.
  • the two Gaussian distribution statistics are, specifically, the mean and variance of the Gaussian distribution of the inner layer side cumulative distribution function and the outer layer side cumulative distribution function.
  • FIG. 8 is a third diagram showing an example of the analysis result of the signal analysis device 1 in the embodiment. More specifically, FIG. 8 is an example of the analysis result by the signal analysis device 1 for the waveform of the electrocardiogram of the target heart whose movement is T extension type 3.
  • FIG. 8 shows three results of graph G8, graph G9 and result G10.
  • the “inner layer side cumulative distribution function” shows the fitting result of a function representing a candidate for a collective channel activity timing distribution of channels existing in the inner layer of the myocardium.
  • the “outer layer side cumulative distribution function” shows the fitting result of a function representing a candidate for a collective channel activity timing distribution of channels existing in the outer layer of the myocardium.
  • the “body surface potential” is a function representing the time change of the body surface potential, and is a waveform of an electrocardiogram.
  • the horizontal axis of FIG. 8 represents time, and the vertical axis represents potential. Both the horizontal axis and the vertical axis are arbitrary units (arbitrary units).
  • Graph G8 represents all of the ECG waveforms generated by one beat.
  • Graph G9 is a part of Graph G8 and shows an enlarged view of a region of a T wave.
  • the region of the T wave is the region shown as region A2 in FIG.
  • Result G10 shows a statistic of two Gaussian distributions.
  • Each value of the result G10 represents a statistic of two Gaussian distributions, that is, the mean and variance of the Gaussian distribution of the inner layer side cumulative distribution function and the outer layer side cumulative distribution function.
  • FIG. 8 shows that ion channel activity information can be acquired only from the electrocardiogram without inserting a catheter by the signal analyzer 1.
  • FIG. 9 is a fourth diagram showing an example of the analysis result of the signal analysis device 1 in the embodiment. More specifically, FIG. 9 is an example of the analysis result by the signal analysis device 1 for the waveform of the electrocardiogram of the target heart whose movement is QT extension type 1.
  • FIG. 9 shows three results of graph G11, graph G12 and result G13.
  • the “inner layer side cumulative distribution function” shows the fitting result of a function representing a candidate for a collective channel activity timing distribution of channels existing in the inner layer of the myocardium.
  • the “outer layer side cumulative distribution function” shows the fitting result of a function representing a candidate for a collective channel activity timing distribution of channels existing in the outer layer of the myocardium.
  • the “body surface potential” is a function representing the time change of the body surface potential, and is a waveform of an electrocardiogram.
  • the horizontal axis of FIG. 9 represents time, and the vertical axis represents potential. Both the horizontal axis and the vertical axis are arbitrary units (arbitrary units).
  • Graph G11 represents all of the ECG waveforms generated by one beat.
  • Graph G12 is a part of graph G11 and shows an enlarged view of a T wave region.
  • the region of the T wave is the region shown as region A3 in FIG.
  • Result G13 shows a statistic of two Gaussian distributions.
  • Each value of the result G13 represents a statistic of two Gaussian distributions.
  • the two Gaussian distribution statistics are, specifically, the mean and variance of the Gaussian distribution of the inner layer side cumulative distribution function and the outer layer side cumulative distribution function.
  • FIG. 9 shows that the ion channel activity information can be acquired only from the electrocardiogram without inserting the catheter by the signal analyzer 1. It is shown that the signal analysis device 1 acquires ion channel activity information only from the electrocardiogram.
  • FIG. 9 is also an example of the estimation result by the signal analysis device 1 of the channel current characteristic related to sudden death.
  • FIGS. 10 to 12 it is shown that ion channel activity information is acquired only from the electrocardiogram by the signal analyzer 1 even for the electrocardiogram of ventricular extrasystole.
  • the horizontal axis of FIGS. 10 to 12 represents time (seconds), and the vertical axis represents potential (mV).
  • FIG. 10 is a first explanatory diagram of an example in which an electrocardiogram of ventricular extrasystole is analyzed by the signal analysis device 1 of the embodiment.
  • FIG. 11 is a second explanatory view of an example in which an electrocardiogram of ventricular extrasystole is analyzed by the signal analysis device 1 of the embodiment.
  • FIG. 12 is a third explanatory diagram of an example in which an electrocardiogram of ventricular premature ventricular contraction is analyzed by the signal analysis device 1 of the embodiment.
  • FIG. 10 shows the electrocardiographic potential of the body surface. That is, FIG. 10 shows a normal single heartbeat and two consecutive ventricular extrasystoles recorded on the electrocardiogram. More specifically, FIG. 11 shows the inner layer side cumulative distribution function of the depolarization of the ventricular extrasystole analyzed by the signal analyzer 1, the outer layer side cumulative distribution function, and the inner layer side cumulative distribution function of the repolarized phase. Is shown. More specifically, FIG. 12 shows an example of the waveform of the premature ventricular contraction of the actually measured electrocardiogram.
  • FIG. 11 shows that the inner layer cumulative distribution function of depolarization precedes the outer layer side cumulative distribution function, and the variance of both is larger than that of a normal heartbeat and the spread of excitement is gradual. This analysis result is consistent with the characteristics of the wide-tailed R-wave waveform.
  • FIG. 11 shows that in the repolarized phase, the inner layer side cumulative distribution function starts inactivation earlier than the outer layer side cumulative distribution function, which is shown as the average value of the cumulative distribution function. There is. These are consistent with the characteristics of the large negative T wave in the repolarization phase, and as shown in FIG. 12, the waveform obtained by the difference in the cumulative distribution function is abbreviated as the waveform of the premature ventricular contraction of the actually measured electrocardiogram. Match.
  • FIGS. 10 to 12 show that the signal analyzer 1 corresponds to an example in which a rogue wave or a negative potential is generated due to a change conduction of myocardial excitement, early repolarization, or a delay in repolarization.
  • the average ⁇ of the inner layer side cumulative distribution function of the depolarized phase is 0.46
  • the variance ⁇ 2 is 1.43.
  • the average ⁇ of the outer layer side cumulative distribution function of depolarization is -1
  • the variance ⁇ 2 is 0.99.
  • the average ⁇ of the inner layer side cumulative distribution function of the repolarized phase is 0.46
  • the variance ⁇ 2 is 1.43.
  • the average ⁇ of the outer layer side cumulative distribution function of the repolarized phase is -1, and the variance ⁇ 2 is 0.99.
  • FIGS. 13 to 15 it is shown that ion channel activity information is acquired only from the electrocardiogram by the signal analyzer 1 for the electrocardiogram of the target heart in the depolarization stage of Brugada syndrome type 1.
  • the vertical axis of FIGS. 13 to 15 represents the potential in millivolts.
  • FIG. 13 is a first explanatory diagram in which the electrocardiogram of the target heart in the depolarization period of Brugada syndrome type 1 is analyzed by the signal analysis device 1 in the embodiment.
  • FIG. 14 is a second explanatory diagram in which the electrocardiogram of the target heart in the depolarization phase of Brugada syndrome type 1 is analyzed by the signal analysis device 1 in the embodiment.
  • FIG. 15 is a third explanatory diagram in which the electrocardiogram of the target heart in the depolarization period of Brugada syndrome type 1 is analyzed by the signal analysis device 1 in the embodiment.
  • FIG. 13 shows an electrocardiogram of the second lead of the chest of Brugada syndrome.
  • the inner frame W1 shows the depolarization phase
  • the inner frame W2 shows the repolarization phase.
  • FIGS. 14 and 15 That is, also in FIGS. 14 and 15, the inner frame W1 represents the depolarized phase, and the inner frame W2 represents the repolarized phase.
  • FIG. 14 shows the depolarization and the inner layer side cumulative distribution function and the outer layer side cumulative distribution function in the repolarized phase analyzed by the signal analyzer 1.
  • the repolarized phase of the inner layer side cumulative distribution function starts following the depolarized phase, which is characteristic of early repolarization.
  • a difference is observed between the depolarized phase and the repolarized phase in the potential amplitude of the outer layer side cumulative distribution function.
  • FIG. 14 shows the gap and anisotropy between the depolarized phase and the repolarized phase of the outer layer side cumulative distribution function.
  • the signal analyzer 1 determines the early repolarization, which is a characteristic of the waveform shown by the electrocardiogram of the target heart of Brugada syndrome, and the anisotropy between depolarization and repolarization, with the inner layer side cumulative distribution function and the outer layer side cumulative distribution function. It can be expressed by the average, dispersion, and ratio of the depolarized phase and the repolarized phase of the distribution function.
  • FIG. 15 is a comparison between the analysis result and the measured value. More specifically, FIG. 15 shows the difference between the inner layer side cumulative distribution function and the outer layer side cumulative distribution function in the depolarization of FIG. 14 and the repolarized phase. FIG. 15 also shows the measured values of the electrocardiogram. Except for the tail, the analysis results and the measured values are almost the same. The tail means a potential at a late time.
  • the average ⁇ of the inner layer side cumulative distribution function of the depolarized phase is 15, and the variance ⁇ 2 is 0.15. Further, in FIGS. 13 to 15, the average ⁇ of the outer layer side cumulative distribution function of depolarization is 14, and the variance ⁇ 2 is 0.25. Further, in FIGS. 13 to 15, the ratio of the inner / outer layers of depolarization is 0.45. Further, in FIGS. 13 to 15, the average ⁇ of the inner layer side cumulative distribution function of the repolarized phase is 25, and the variance ⁇ 2 is 0.25. Further, in FIGS. 13 to 15, the average ⁇ of the outer layer side cumulative distribution function of the repolarized phase is 20, and the variance ⁇ 2 is 0.5. Further, in FIGS. 13 to 15, the ratio of the inner and outer layers of depolarization is 0.75. Further, in FIGS. 13 to 15, the polarization / repolarization ratio is 1.4.
  • FIGS. 16 to 54 show the results of analysis by the signal analysis device 1 using the publicly available electrocardiogram data library https://physionet.org/about/database/. 16 to 54 show the inner layer side cumulative distribution function, the outer layer side cumulative distribution function, and the fitting result obtained as a result of the execution of the analysis of the electrocardiographic potential by the signal analysis device 1.
  • 16 to 54 are diagrams showing an example of analyzing an electrocardiogram by the signal analysis device 1 in the embodiment, respectively.
  • the determined points shown in each figure represent the Q point, the R point, the S point, the T start point, and the T end point of the electrocardiographic potential, respectively, in order from the left side of the figure.
  • the determined point was determined by the inflection detection and peak detection algorithms.
  • 16 to 54 show the inner layer side cumulative distribution function and the outer layer side cumulative distribution in each section obtained for the depolarized phase (QRS wave) section and the repolarized phase (T wave) section. Indicates a function.
  • FIGS. 16 to 54 show that the signal analyzer 1 can show substantially the same shape for various QRS waves and T waves by adjusting the average and variance of the inner layer side cumulative distribution function and the outer layer side cumulative distribution function. Show that.
  • the figures below FIGS. 16 to 54 show the original waveform of the electrocardiogram and the fitting result.
  • the results of FIGS. 16 to 54 are the results of sampling at 300 Hz. Therefore, the horizontal axis of each of FIGS. 16 to 54 represents 0 seconds at the origin and 3.33 milliseconds at the value 1.
  • the signal analysis device 1 configured in this way fits the waveform of the electrocardiogram of the target heart using one or a plurality of channel cumulative distribution functions, and obtains ion channel activity information based on the shape of each channel cumulative distribution function of the fitting result. get.
  • Each channel cumulative distribution function used for fitting is information indicating candidates for channel activity timing distribution.
  • the channel activity timing distribution represents the distribution of ion channel activity timings, and the different channel cumulative distribution functions represent different types of channel activity timing distributions. Therefore, each shape of each channel cumulative distribution function of the fitting result indicates information of the channel activity timing distribution represented by each channel cumulative distribution function. Such information cannot be obtained by conventional electrocardiogram waveform analysis. Therefore, the signal analysis device 1 can increase the information obtained from the waveform of the electrocardiogram.
  • the number of channel cumulative distribution functions in which the domains overlap does not necessarily have to be two, and may be three or more.
  • fitting may be performed to minimize the difference between the linear sum of the cumulative distribution functions of a plurality of channels having the overlap of the domain and satisfying a predetermined condition regarding the weight and the waveform of the electrocardiogram. ..
  • the predetermined condition regarding the weight may be any condition as long as the condition that at least one has a different sign from the other weight is included.
  • the electrocardiogram is preferably an induced electrocardiogram close to the electromotive force vector.
  • the lead electromotive force close to the electromotive force vector for example, an electrocardiogram capable of acquiring three-dimensional information on the electrocardiogram such as lead II, lead V4, and lead V5 is preferable. Since the amount of information increases as the number of channels of the electrocardiogram increases, it is desirable that the number of channels of the electrocardiogram increases.
  • the signal analysis device 1 may be mounted by using a plurality of information processing devices connected so as to be communicable via a network.
  • each functional unit included in the signal analysis device 1 may be distributed and mounted in a plurality of information processing devices.
  • the electrocardiogram acquisition unit 110 is an example of the biological information acquisition unit.
  • each function of the signal analysis device 1 may be realized by ASIC (Application Specific Integrated Circuit), PLD (Programmable Logic Device), FPGA (Field Programgable Gate Array), or the like. ..
  • the program may be recorded on a computer-readable recording medium.
  • the computer-readable recording medium is, for example, a flexible disk, a magneto-optical disk, a portable medium such as a ROM or a CD-ROM, or a storage device such as a hard disk built in a computer system.
  • the program may be transmitted over a telecommunication line.
  • Electrocardiogram acquisition unit 120 ... Fitting information acquisition unit, 130 ... Analysis unit, 131 ... Fitting section, 132 ... Ion channel activity information acquisition section, 140 ... Recording section, 91 ... Processor, 92 ... Memory

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Physiology (AREA)
  • Signal Processing (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本発明の一態様は、解析対象の心臓の拍動に関する時系列の生体情報を取得する生体情報取得部と、前記心臓におけるイオンチャネルの活動のタイミングの分布であるチャネル活動タイミング分布の候補を示す情報を分布候補情報として、前記生体情報と前記分布候補情報とに基づき、前記イオンチャネルの活動の様子を示す情報を取得する解析部と、を備える信号解析装置である。

Description

信号解析装置、信号解析方法及びプログラム
　本発明は、信号解析装置、信号解析方法及びプログラムに関する。
　心電図は心臓の状態を把握するために有用な情報であり、例えば心電図を用いれば対象に心不全が起きる可能性が高い状態にあるか否かを判定することができる（非特許文献１）。
田中博、「心電図逆問題における方法論」、医用電子と生体工学、１９８５年、２３巻、３号、ｐ．１４７－１５８ 清水渉、「Ｔ波の成因におけるＭ細胞の役割」、ＪＰＮ．Ｊ．ＥＬＥＣＴＲＯＣＡＲＤＩＯＬＯＧＹ、Ｖｏｌ．２１、Ｎｏ．２、２００１年、ｐ.１０１～１０８ 田中義文、野土信司、山崎康夫、仲宗根ありさ、「体表心電図を心筋活動電位波形より構築するプログラム」、麻酔・集中治療とテクノロジー２０１１（１）：p.９１－９５、２０１１
　しかしながら、心電図の波形から得られる情報は、心臓の状態の把握に必ずしも充分では無い場合がある。例えば、心電図の波形が似たような波形であっても心臓に関わる病の発症の仕方が異なる場合がある。このように、心電図の波形から得られる情報は心臓の状態の把握に必ずしも充分では無い場合があるものの、病によっては他の情報を取得するのが難しい場合がある。例えば心不全であれば、その抑制には日常生活の中で心電図の波形により心臓の状態を観測することで発症を抑制することができる。発症を抑制する精度をより一層上げるには血液を採取する等、他の技術を用いて他の情報を取得することが考えられるが、日常生活でそれを行うことは現実的ではない。そのため、病によっては実質的に心電図の波形だけに基づいて心臓の状態を把握しなければならない場合がある。
　またこのような事情は、心電図の波形に基づいて心臓の状態を把握する場合に限らない。このような事情は、心臓の拍動に関する時系列の生体情報だけに基づいて心臓の状態を把握する場合にも共通する。なお、心臓の拍動に関する時系列の生体情報は、例えば心電位の変化を示す波形や、心臓の画像の変化を示す波形や、心臓の圧の変化を示す波形や、血流量の変化を示す波形や、心音の変化を示す波形である。なお心電図の波形も心臓の拍動に関する時系列の生体情報の一例である。なお心臓の画像は、例えば心臓のＭＲＩ（Magnetic Resonance Image)画像や心臓のＣＴ（Computed Tomography）画像である。
　上記事情に鑑み、本発明は、心臓の拍動に関する時系列の生体情報から得られる情報を増やす技術を提供することを目的としている。
　本発明の一態様は、解析対象の心臓の拍動に関する時系列の生体情報を取得する生体情報取得部と、前記心臓におけるイオンチャネルの活動のタイミングの分布であるチャネル活動タイミング分布の候補を示す情報を分布候補情報として、前記生体情報と前記分布候補情報とに基づき、前記イオンチャネルの活動の様子を示す情報を取得する解析部と、を備える信号解析装置である。
　本発明により、心臓の拍動に関する時系列の生体情報から得られる情報を増やす技術を提供することが可能となる。
実施形態の信号解析装置１のハードウェア構成の一例を示す図。 実施形態における対象心臓の心電図の波形を４つのチャネル累積分布関数によってフィッティングした結果の一例を示す図。 実施形態における２つのチャネル累積分布関数の差分がＲ波の立ち上がりの波形に略同一の波形にフィッティング可能であることを説明する説明図。 実施形態における制御部１１の機能構成の一例を示す図。 実施形態における信号解析装置１が実行する処理の流れの一例を示すフローチャート。 実施形態における信号解析装置１の解析結果の一例を示す第１の図。 実施形態における信号解析装置１の解析結果の一例を示す第２の図。 実施形態における信号解析装置１の解析結果の一例を示す第３の図。 実施形態における信号解析装置１の解析結果の一例を示す第４の図。 実施形態の信号解析装置１によって心室性期外収縮の心電図を解析した一例の第１の説明図。 実施形態の信号解析装置１によって心室性期外収縮の心電図を解析した一例の第２の説明図。 実施形態の信号解析装置１によって心室性期外収縮の心電図を解析した一例の第３の説明図。 実施形態における信号解析装置１によってブルガダ症候群１型の脱分極期における対象心臓の心電図を解析した第１の説明図。 実施形態における信号解析装置１によってブルガダ症候群１型の脱分極期における対象心臓の心電図を解析した第２の説明図。 実施形態における信号解析装置１によってブルガダ症候群１型の脱分極期における対象心臓の心電図を解析した第３の説明図。 実施形態における信号解析装置１によって心電図を解析した第１の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第４の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第５の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第６の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第７の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第８の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第９の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１０の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１１の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１２の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１３の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１４の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１５の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１６の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１７の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１８の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第１９の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２０の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２１の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２２の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２３の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２４の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２５の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２６の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２７の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２８の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第２９の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３０の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３１の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３２の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３３の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３４の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３５の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３６の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３７の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３８の例を示す図。 実施形態における信号解析装置１によって心電図を解析した第３９の例を示す図。
（第１実施形態）
　図１は、実施形態の信号解析装置１のハードウェア構成の一例を示す図である。以下、説明の簡単のため、少なくとも心電図の波形に基づき解析を行う場合を例に信号解析装置１を説明する。しかしながら、信号解析装置１は心電図に限らず心臓の拍動に関する時系列の生体情報に基づけば同様の解析が可能である。なお心臓の拍動に関する時系列の生体情報は、例えば心電位の変化を示す波形や、心臓の画像の変化を示す波形や、心臓の圧の変化を示す波形や、血流量の変化を示す波形や、心音の変化を示す波形である。そのため、信号解析装置１は心電図の波形に代えて心周期を示す波形を用いてもよい。信号解析装置１は心電図の波形に代えて心臓の画像の変化を示す波形を用いてもよい。
　信号解析装置１は心電図の波形に代えて心臓の圧の変化を示す波形を用いてもよい。信号解析装置１は心電図の波形に代えて血流量の変化を示す波形を用いてもよい。信号解析装置１は心電図の波形に代えて心音の変化を示す波形を用いてもよい。なお、心電図も心臓の拍動に関する時系列の生体情報の一例である。なお、心臓の拍動に関する時系列の生体情報は、心臓の周期的な拍動に関する時系列の生体情報であってもよい。なお心臓の画像は、例えば心臓のＭＲＩ（Magnetic Resonance Image)画像や心臓のＣＴ（Computed Tomography）画像である。
　信号解析装置１は、解析対象の心臓（以下「対象心臓」という。）の心電図の波形を取得する。信号解析装置１は取得した心電図の波形に基づき、対象心臓の心筋の外側の層の心筋と内側の層の心筋とのイオンチャネルの活動の様子を示す情報（以下「イオンチャネル活動情報」という。）を取得する。
　ここで、イオンチャネルと心電図の波形との関係について説明する。医学の分野では心起電力双極子モデル（参考文献１）という心筋の動きと心電図と関係を説明するモデルが知られている。心起電力双極子モデルによれば、心筋は外側の層（以下「心筋外層」という。）と内側の層（以下「心筋内層」という。）の２層でモデル化される。
　参考文献１：田中義文「成り立ちから理解する心電図波形　心筋の活動電位を読み解く」　学研メディカル秀潤社（２０１２）
　心起電力双極子モデルでは、心筋外層と心筋内層とはそれぞれ異なる起電力の発生源としてモデル化される。心起電力双極子モデルによれば、心外膜側心筋活動電位と心内膜側心筋活動電位との合成波が体表面で観測される体表面の電位の時間変化に略一致する。体表面の電位の時間変化を表すグラフが心電図の波形である。心外膜側心筋活動電位は、カテーテルの挿入を行うことで心筋外層の拍動によって生じた起電力の変化を直接測定した結果である。心内膜側心筋活動電位は、カテーテルの挿入を行うことで心筋内層の拍動によって生じた起電力の変化を直接測定した結果である。ここまでが心起電力双極子モデルの概略の説明である。
　ところで、心起電力双極子モデルにおける心筋外層は細胞の集合である。そのため心筋外層の１回の拍動における心筋外層の細胞の拍動のタイミングは必ずしも全ての細胞で同一ではなく、拍動のタイミングには分布が存在する可能性がある。このことは心筋内層についても同様である。すなわち、そのため心筋内層の１回の拍動における心筋内層の細胞の拍動のタイミングは必ずしも全ての細胞で同一ではなく、拍動のタイミングには分布が存在する可能性がある。しかしながらこのような、細胞の拍動のタイミングに分布が存在する可能性は心起電力双極子モデルでは想定されていない。
　また、１回の心筋の拍動では細胞ごとに各細胞の拍動のタイミングに応じたタイミングで複数種類のイオンチャネルが活動を行う。イオンチャネルの活動のタイミングは、イオンチャネルの種類の応じたタイミングであり必ずしも同一ではない。例えばナトリウムチャネルとカリウムチャネルとは活動のタイミングが異なる。そのため、細胞の拍動のタイミングに分布が存在すれば、１回の心筋の拍動におけるイオンチャネルの活動のタイミングにも細胞の拍動のタイミングの分布に応じた分布が存在する。
　そして、イオンチャネルの活動のタイミングはイオンチャネルの種類の応じたタイミングであるため、細胞の拍動のタイミングに分布が存在すれば、異なるタイミングで活動するイオンチャネルについては、チャネル活動タイミング分布が異なるタイミングで現れる。チャネル活動タイミング分布は、イオンチャネルの活動のタイミングの分布である。
　チャネル活動タイミング分布は、例えばガウス分布である。チャネル活動タイミング分布は、チャネル活動タイミング分布を表す確率密度関数の累積分布関数がシグモイド曲線を表す関数である、という条件を満たすものであってもよい。シグモイド曲線を表す関数は、例えばシグモイド関数であってもよいし、累積正規分布関数であってもよいし、ゴンペルツ関数であってもよいし、グーデルマン関数であってもよい。チャネル活動タイミング分布は、例えばロジスティック分布であってもよい。
　信号解析装置１は、チャネル活動タイミング分布の候補を示す情報（以下「分布候補情報」という。）を用いて、取得した対象心臓の心電図の波形に基づきイオンチャネル活動情報を取得する。分布候補情報は、例えばチャネル活動タイミング分布の候補を表す関数を確率密度関数とする累積分布関数（以下「チャネル累積分布関数」という。）で表現される。チャネル活動タイミング分布の候補を表す関数は、具体的にはチャネル活動タイミング分布の分散や平均値等のチャネル活動タイミング分布の形状を表す値がパラメータで表現された関数である。以下、チャネル活動タイミング分布の形状を表すパラメータを形状パラメータという。そのため、分布候補情報を表現する累積分布関数は１又は複数の形状パラメータを有する汎関数である。
　以下、分布候補情報を用いて、取得した対象心臓の心電図の波形に基づきイオンチャネル活動情報を取得する処理を、イオンチャネル活動情報取得処理という。取得されたイオンチャネル活動情報が、信号解析装置１による対象心臓の心電図の波形の解析結果の一例である。
　以下、説明の簡単のため分布候補情報がチャネル累積分布関数である場合を例に、信号解析装置１を説明する。このような場合イオンチャネル活動情報取得処理は具体的には、１又は複数のチャネル累積分布関数を用いて対象心臓の心電図の波形に対してフィッティングを行うことで、チャネル累積分布関数が有する形状パラメータの値を決定する処理である。
　チャネル累積分布関数は、例えば以下の式（１）又は式（２）で表される。
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000002
　式（１）は正規分布（ガウス分布）である。式（２）は式（１）の累積分布関数である。式（２）において"erf"は誤差関数（シグモイド関数）である。
　少なくとも平均、分散のいずれか一方が異なる２つの累積分布関数について、平均をそれぞれμａ、μｂ、と表現し、分散をそれぞれσ２ ａ、σ２ ｂと表現すると、２つの累積分布関数の差は以下の式（３）で表される。
Figure JPOXMLDOC01-appb-M000003
　例えば後述の図３の画像Ｇ１に記載の第１関数は以下の式（４）で表される。
Figure JPOXMLDOC01-appb-M000004
　例えば後述の図３の第２関数は以下の式（５）で表される。
Figure JPOXMLDOC01-appb-M000005
　例えば後述の図３の第３関数は以下の式（６）で表される。
Figure JPOXMLDOC01-appb-M000006
　この第３関数を心電位のＲ波やＴ波に近似させることで、第１関数と第２関数との累積分布関数の平均及び分散が得られる。第３関数を心電位のＲ波やＴ波に近似させる処理は、例えば最小二乗法により心電位と第３関数との差を最小化する処理である。
　フィッティングに用いられるチャネル累積分布関数の定義域は必ずしも、フィッティング対象である心電図の波形を表す関数の定義域と同一でなくてもよい。また、フィッティングに用いられるチャネル累積分布関数の各定義域は必ずしも、同一で無くてもよい。必ずしも同一で無い、とは少なくとも１つが異なってもよいし、全てが同一であってもよい、ことを意味する。
　イオンチャネル活動情報取得処理におけるフィッティングは、具体的には形状パラメータを変数として、チャネル累積分布関数の定義域内においてチャネル累積分布関数と心電図の波形との違いを最小にする最適化を行う処理である。
（チャネル累積分布関数を用いたフィッティングの結果の意味について）
　ここで、チャネル累積分布関数を用いたフィッティング結果の意味について説明する。形状パラメータの値が固定されている場合、チャネル累積分布関数は確率密度関数の積分の結果であり確率密度関数に１対１に対応する。そして、形状パラメータの値が固定されている場合、確率密度関数はチャネル活動タイミング分布を表す関数である。したがって、形状パラメータの値が固定されている場合、チャネル累積分布関数はチャネル活動タイミング分布を表す関数である。そしてそのため、形状パラメータの値が固定されていない場合、チャネル累積分布関数は、形状パラメータの値の組合せの数だけのチャネル活動タイミング分布の候補を表す関数である。
　ところで一般に、フィッティングとは関数の形状を表すパラメータの値を変数として最適化を行う処理である。したがって、フィッティングにチャネル累積分布関数を用いてフィッティングを実行することは、チャネル累積分布関数が表すチャネル活動タイミング分布の候補から、フィッティング対象との違いを最小にするチャネル活動タイミング分布を決定する処理である。したがって、フィッティングにチャネル累積分布関数を用いてフィッティングを実行することは、イオンチャネルの活動がどのようにフィッティング対象と関わるかを示す情報を取得することに相当する。このため、フィッティング対象が心電図の波形である場合には、フィッティングにより、イオンチャネルの活動が心電図の波形とどのように関わるかを示す情報が取得される。
　イオンチャネルの活動と心電図の波形との関わりは、例えばフィッティング結果の形状パラメータの値によって示される。そのため、フィッティング結果の形状パラメータの値は、イオンチャネル活動情報の一例である。
　ここまでフィッティングにおいてチャネル累積分布関数が１つ用いられる場合を説明したが、フィッティングには、複数のチャネル累積分布関数が用いられてもよい。このような場合においてフィッティングに用いられる複数のチャネル累積分布関数それぞれは、互いに異なる種類のチャネル活動タイミング分布を表す。このような場合、複数種類のイオンチャネルの活動について心電図の波形とどのように関わるかを示す情報が取得される。なお、チャネル活動タイミング分布の種類が異なるとは、イオンチャネルの種類と、イオンチャネルの開閉の動作の違いとのいずれか一方又は両方が異なることを意味する。
　フィッティングに用いられるチャネル累積分布関数の定義域は必ずしも、フィッティング対象である心電図の波形を表す関数の定義域と同一でなくてもよい。また、複数のチャネル累積分布関数がフィッティングに用いられる場合、チャネル累積分布関数それぞれの各定義域は必ずしも、同一で無くてもよい。必ずしも同一で無い、とは少なくとも１つが異なってもよいし、全てが同一であってもよい、ことを意味する。
　図２は、実施形態における対象心臓の心電図の波形を４つのチャネル累積分布関数によってフィッティングした結果の一例を示す図である。図２の横軸は時刻を表し、縦軸は電位を表す。横軸と縦軸との単位はいずれも任意単位（arbitrary unit）である。
　図２における４つのチャネル累積分布関数は、具体的には、第１累積分布関数、第２累積分布関数、第３累積分布関数及び第４累積分布関数である。図２における第１累積分布関数、第２累積分布関数、第３累積分布関数及び第４累積分布関数の定義域は少なくとも１つが異なる。具体的には、第１累積分布関数の定義域と第２累積分布関数の定義域とは同一であり、時刻０から時刻Ｔ１である。第３累積分布関数の定義域は、時刻Ｔ３から時刻Ｔ５である。第４累積分布関数の定義域は時刻Ｔ２から時刻Ｔ４である。なお、フィッティング対象である心電図の波形の定義域は、少なくとも時刻０から時刻Ｔ５の期間を含む。
　第１累積分布関数は、心筋内層に存在するナトリウムチャネルのチャネル活動タイミング分布の候補であって心筋の興奮の開始のタイミングにおけるチャネル活動タイミング分布の候補を表す。
　第２累積分布関数は、心筋外層に存在するカリウムチャネルのチャネル活動タイミング分布の候補であって心筋の興奮の開始のタイミングのチャネル活動タイミング分布の候補を表す。
　第３累積分布関数は、心筋内層に存在するナトリウムチャネルのチャネル活動タイミング分布の候補であって心筋の興奮の抑制の開始のタイミングのチャネル活動タイミング分布の候補を表す。
　第４累積分布関数は、心筋外層に存在するカリウムチャネルのチャネル活動タイミング分布の候補であって心筋の興奮の抑制の開始のタイミングのチャネル活動タイミング分布の候補を表す。
　図２において“第１フィッティング結果”は、第１累積分布関数による心電図の波形へのフィッティングの結果を示す。以下、第１累積分布関数による心電図の波形へのフィッティングを、第１フィッティングという。図２において“第２フィッティング結果”は、第２累積分布関数による心電図の波形へのフィッティングの結果を示す。以下、第２累積分布関数による心電図の波形へのフィッティングを、第２フィッティングという。図２において“第３フィッティング結果”は、第３累積分布関数による心電図の波形へのフィッティングの結果を示す。以下、第３累積分布関数による心電図の波形へのフィッティングを、第３フィッティングという。
　図２において“第４フィッティング結果”は、第４累積分布関数による心電図の波形へのフィッティングの結果を示す。以下、第４累積分布関数による心電図の波形へのフィッティングを、第４フィッティングという。図２において“体表面の電位”は、フィッティング対象の心電図の波形を表す。図２において第１累積分布関数、第２累積分布関数、第３累積分布関数及び第４累積分布関数はいずれも、具体的には累積ガウス分布である。
　“第１フィッティング結果”と“第２フィッティング結果”とは、Ｒ波に対するフィッティングの結果である。“第３フィッティング結果”と“第４フィッティング結果”とは、Ｔ波に対するフィッティングの結果である。
　なお、時刻Ｔ１から時刻Ｔ２の期間は第１累積分布関数から第４累積分布関数までのいずれのチャネル累積分布関数の定義域にも属さない。そのため、図２に示す心電図の波形のうち時刻Ｔ１から時刻Ｔ２の期間の波形についてはフィッティングが行われない。図２においてフィッティングが行われない期間については、定数関数で表現されている。フィッティングが行われない期間については、定数関数に限らず一次関数等の予め定められた関数であればどのような関数で表現されてもよい。
　図２に示すフィッティングの結果のうち、“第１フィッティング結果”から“第２フィッティング結果”を引き算した結果はＲ波の立ち上がりの期間の波形に略同一の波形を示す。図２に示すフィッティングの結果のうち、“第３フィッティング結果”から“第４フィッティング結果”を引き算した結果はＴ波の立ち上がりの期間の波形に略同一の波形を示す。
　なお、図２のＲ波に対するフィッティング結果は、“第１フィッティング結果”から“第２フィッティング結果”を引き算した結果はＲ波の立ち上がりの期間の波形に略同一の波形を示すという条件、の下にフィッティングが行われた結果であってもよい。なお、図２のＴ波に対するフィッティング結果は、“第３フィッティング結果”から“第４フィッティング結果”を引き算した結果はＴ波の立ち上がりの期間の波形に略同一の波形を示すという条件、の下にフィッティングが行われた結果であってもよい。
　このように定義域の重なりを有する２つのチャネル累積分布関数については、重なる定義域に関して、関数間の差が心電図の波形にフィッティングされてもよい。すなわち、定義域の重なりを有する２つのチャネル累積分布関数については関数間の差と心電図の波形との違いを最小にするフィッティングが行われてもよい。
　２つのチャネル累積分布関数の差がＲ波の立ち上がりの波形に略同一であることについて、図３を用いて説明する。
　図３は、実施形態における２つのチャネル累積分布関数の差分がＲ波の立ち上がりの波形に略同一の波形にフィッティング可能であることを説明する説明図である。図３は、画像Ｇ１、画像Ｇ２、画像Ｇ３及び画像Ｇ４の４つの画像を示す。画像Ｇ１～画像Ｇ４の各画像は、横軸を時刻、縦軸を電位とするグラフを示す。図３の画像Ｇ１～画像Ｇ４の各画像の横軸と縦軸との単位はいずれも任意単位（arbitrary unit）である。
　図３の“第１関数”は、累積分布関数の一例である。図３の“第２関数”は、累積分布関数の一例であり“第２関数”とは異なる累積分布関数である。図３の“第３関数”は、“第１関数”から“第２関数”を引き算した関数を表す。図３の“第３関数”は、Ｒ波の形状に略同一の形状であることを示す。
　図３の“第４関数”は、累積分布関数の一例である。図３の“第５関数”は、累積分布関数の一例であり“第４関数”とは異なる累積分布関数である。図３の“第６関数”は、“第４関数”から“第５関数”を引き算した関数を表す。図３の“第６関数”は、Ｒ波の形状に略同一の形状であることを示す。
　図３の“第７関数”は、累積分布関数の一例である。図３の“第８関数”は、累積分布関数の一例であり“第７関数”とは異なる累積分布関数である。図３の“第９関数”は、“第８関数”から“第７関数”を引き算した関数を表す。図３の“第９関数”は、Ｒ波の形状に略同一の形状であることを示す。
　図３の“第１０関数”は、累積分布関数の一例である。図３の“第１１関数”は、累積分布関数の一例であり“第１０関数”とは異なる累積分布関数である。図３の“第１２関数”は、“第１０関数”から“第１１関数”を引き算した関数を表す。図３の“第１２関数”は、Ｒ波の形状に略同一の形状であることを示す。
　このように、２つの累積分布関数の差の関数は、Ｒ波の形状に略同一の形状を表現可能である。このことは、Ｔ波の形状についても同様である。すなわち、２つの累積分布関数の差の関数は、Ｔ波の形状に略同一の形状を表現可能である。２つの累積分布関数の差の関数は、“第３関数”、“第６関数”、“第９関数”のように縦軸方向の幅と横軸方向の幅との異なる波を表現可能である。また２つの累積分布関数の差の関数は、“第１２関数”のように陰性の波を表現可能である。
　このように、信号解析装置１は、チャネル累積分布関数によって対象心臓の心電図の波形をフィッティングする。そして信号解析装置１は、フィッティング結果に基づき、イオンチャネル活動情報を取得する。
　図１の説明に戻る。信号解析装置１は、バスで接続されたＣＰＵ等のプロセッサ９１とメモリ９２とを備える制御部１１を備え、プログラムを実行する。信号解析装置１は、プログラムの実行によって制御部１１、入力部１２、通信部１３、記憶部１４及び出力部１５を備える装置として機能する。
　より具体的には、プロセッサ９１が記憶部１４に記憶されているプログラムを読み出し、読み出したプログラムをメモリ９２に記憶させる。プロセッサ９１が、メモリ９２に記憶させたプログラムを実行することによって、信号解析装置１は、制御部１１、入力部１２、通信部１３、記憶部１４及び出力部１５を備える装置として機能する。
　制御部１１は、信号解析装置１が備える各種機能部の動作を制御する。制御部１１は、例えばイオンチャネル活動情報取得処理を実行する。制御部１１は、例えば出力部１５の動作を制御し、出力部１５にイオンチャネル活動情報取得処理の取得結果を出力させる。制御部１１は、例えばイオンチャネル活動情報取得処理の実行により生じた各種情報を記憶部１４に記録する。
　入力部１２は、マウスやキーボード、タッチパネル等の入力装置を含んで構成される。入力部１２は、これらの入力装置を信号解析装置１に接続するインタフェースとして構成されてもよい。入力部１２は、信号解析装置１に対する各種情報の入力を受け付ける。
　入力部１２には、例えばフィッティングに用いるチャネル累積分布関数の数を示す情報（以下「数指定情報」という。）が入力される。入力部１２には、例えばフィッティングに用いる各チャネル累積分布関数について、各チャネル累積分布関数が表す分布の形状を示す情報（以下「分布形状指定情報」という。）が入力される。入力部１２には、例えばフィッティングに用いる各チャネル累積分布関数について定義域を示す情報（以下「定義域指定情報」という。）が入力される。
　なお、数指定情報、分布形状指定情報又は定義域指定情報の一部は予め記憶部１４が記憶済みであってもよい。このような場合、記憶部１４が記憶済みの情報であって数指定情報、分布形状指定情報又は定義域指定情報の一部の情報については、入力部１２から入力される必要は無い。以下、説明の簡単のため記憶部１４が数指定情報、分布形状指定情報又は定義域指定情報を予め記憶済みの場合を例に信号解析装置１を説明する。
　通信部１３は、信号解析装置１を外部装置に接続するための通信インタフェースを含んで構成される。通信部１３は、有線又は無線を介して外部装置と通信する。外部装置は、例えば対象心臓の心電図の波形の送信元の装置である。対象心臓の心電図の波形の送信元の装置は、例えば心電図測定装置である。通信部１３は、例えば外部装置が心電図測定装置である場合、通信により心電図測定装置から心電図の波形を取得する。なお、心電図の波形は入力部１２に入力されてもよい。
　記憶部１４は、磁気ハードディスク装置や半導体記憶装置などの非一時的コンピュータ読み出し可能な記憶媒体装置を用いて構成される。記憶部１４は信号解析装置１に関する各種情報を記憶する。記憶部１４は、例えば入力部１２又は通信部１３を介して入力された情報を記憶する。記憶部１４は、例えば入力部１２又は通信部１３を介して入力された心電図を記憶する。記憶部１４は、例えばイオンチャネル活動情報取得処理の実行により生じた各種情報を記憶する。
　出力部１５は、各種情報を出力する。出力部１５は、例えばＣＲＴ（Ｃａｔｈｏｄｅ　Ｒａｙ　Ｔｕｂｅ）ディスプレイや液晶ディスプレイ、有機ＥＬ（Ｅｌｅｃｔｒｏ－Ｌｕｍｉｎｅｓｃｅｎｃｅ）ディスプレイ等の表示装置を含んで構成される。出力部１５は、これらの表示装置を信号解析装置１に接続するインタフェースとして構成されてもよい。出力部１５は、例えば入力部１２に入力された情報を出力する。出力部１５は、例えば入力部１２又は通信部１３に入力された心電図を表示してもよい。出力部１５は、例えばイオンチャネル活動情報取得処理の実行結果を表示してもよい。
　図４は、実施形態における制御部１１の機能構成の一例を示す図である。制御部１１は、心電図取得部１１０、フィッティング情報取得部１２０、解析部１３０及び記録部１４０を備える。
　心電図取得部１１０は、入力部１２又は通信部１３に入力された対象心臓の心電図の波形を取得する。フィッティング情報取得部１２０は、数指定情報、分布形状指定情報及び定義域指定情報を取得する。数指定情報、分布形状指定情報及び定義域指定情報が記憶部１４に記憶されている場合、フィッティング情報取得部１２０は、数指定情報、分布形状指定情報及び定義域指定情報を記憶部１４から読み出す。
　解析部１３０は、フィッティング部１３１及びイオンチャネル活動情報取得部１３２を備える。
　フィッティング部１３１は、数指定情報、分布形状指定情報及び定義域指定情報によって示されるチャネル累積分布関数を用いて、心電図取得部１１０が取得した心電図の波形に対してフィッティングを行う。
　イオンチャネル活動情報取得部１３２は、フィッティング部１３１によるフィッティングの結果に基づきイオンチャネル活動情報を取得する。イオンチャネル活動情報取得部１３２は、例えば各チャネル累積分布関数について形状パラメータの値を取得する。イオンチャネル活動情報取得部１３２が取得した形状パラメータの値がイオンチャネル活動情報の一例である。
　このように解析部１３０は、対象心臓の心電図の波形と分布候補情報とに基づき、イオンチャネル活動情報を取得する。
　記憶部１４は、制御部１１が実行した処理によって生じた各種情報を記憶部１４に記録する。
　図５は、実施形態における信号解析装置１が実行する処理の流れの一例を示すフローチャートである。心電図取得部１１０が、入力部１２又は通信部１３を介して対象心臓の心電図を取得する（ステップＳ１０１）。次にフィッティング情報取得部１２０が、数指定情報、分布形状指定情報及び定義域指定情報を取得する（ステップＳ１０２）。次にフィッティング部１３１が数指定情報、分布形状指定情報及び定義域指定情報によって示されるチャネル累積分布関数を用いて、ステップＳ１０１で取得された心電図の波形に対してフィッティングを行う（ステップＳ１０３）。次に、イオンチャネル活動情報取得部１３２が、フィッティング結果に基づき、イオンチャネル活動情報を取得する（ステップＳ１０４）。取得されたイオンチャネル活動情報が出力部１５に出力される（ステップＳ１０５）。
　ステップＳ１０５では、各チャネル累積分布関数のフィッティング結果のグラフを表示してもよい。また、ステップＳ１０２の処理は、ステップＳ１０３の処理の実行前に実行されればよく、ステップＳ１０１の実行前であってもよい。ステップＳ１０３及びステップＳ１０４の処理が解析部１３０の実行する処理の一例である。
　図６は、実施形態における信号解析装置１の解析結果の一例を示す第１の図である。より具体的には、図６は、動作が正常である対象心臓の心電図の波形に対する信号解析装置１による解析結果の一例である。図６の横軸は時刻を表し、縦軸は電位を表す。縦軸の単位は、任意単位（arbitrary unit）である。
　図６は、動作が正常である対象心臓の心電図の脱分極のＲ波形に対して第１フィッティング及び第２フィッティングを行った結果を示す。また図６は、動作が正常である対象心臓の心電図の再分極相のＴ波形に対して第４フィッティング及び第３フィッティングを行った結果の一例を示す。
　図６は、正常な心臓の脱分極においては、心筋内層は心筋外層よりも早くイオンチャネルの活動が開始し急速に進むことを示す。図６は、心筋内層のイオンチャネルの活動が開始するタイミングよりも僅かに遅れて心筋外層が活動を開始する、ことを示す。図６は、心筋内層のイオンチャネルの活動が開始するタイミングと心筋外層が活動を開始するタイミングとの差分は陽性で鋭いＲ波形を示す、ことを示す。第１フィッティング結果である内層側累積分布関数と第２フィッティング結果である外層側累積分布関数の平均値と分散とは、それぞれ活動のタイミングと活動の進み具合を表す。
　図６は、正常な心臓の再分極相においては、心筋外層は心筋内層よりも早くイオンチャネルの活動の不活性化が開始することを示す。図６は、心筋内層の不活性化はそれに遅れ、両者は緩やかに進む、ことを示す。その差分は陽性で緩やかなＴ波形を示す。第４フィッティング結果である内層側累積分布関数と第３フィッティング結果である外層側累積分布関数の平均値と分散は、それぞれ不活性化のタイミングと不活性化の進み具合を表す。以上のように、第１から４のフィッティングにより得られた内層側累積分布関数と外層側累積分布関数の平均と分散は、脱分極のイオンチャネルの集合的な活性化と、再分極相のイオンチャネルの集合的な不活性化のタイミングと進み具合を表す。
　図６の内層側累積分布関数及び外層側累積分布関数の形状は、動作が正常である対象心臓の心筋に対して、カテーテルの挿入を行うことで直接測定した心筋の起電力の計測結果に略一致する。そのため、図６は、信号解析装置１によりカテーテルを挿入することなく、心電図だけからイオンチャネル活動情報が取得されることを示す。
　図７は、実施形態における信号解析装置１の解析結果の一例を示す第２の図である。より具体的には、図７は、動作が正常である対象心臓の心電図の波形に対する信号解析装置１による解析結果の一例である。
　図７は、グラフＧ５、グラフＧ６及び結果Ｇ７の３つの結果を示す。図７において、“内層側累積分布関数”は心筋内層に存在するイオンチャネルの集合的なチャネル活動タイミング分布の候補を表す関数のフィッティング結果を示す。図７において、“外層側累積分布関数”は心筋外層に存在するイオンチャネルの集合的なチャネル活動タイミング分布の候補を表す関数のフィッティング結果を示す。図７において、“体表面の電位”は、体表面の電位の時間変化を表す関数であって、心電図の波形である。図７の横軸は時刻を表し、縦軸は電位を表す。横軸と縦軸との単位はいずれも任意単位（arbitrary unit）である。なお、図７～図９の各横軸の１目盛の間隔が表す時間の長さは同じである。また、図７～図９の縦軸はいずれも、１が累積ガウス分布の最大値を表す。
　グラフＧ５は、１回の拍動で生じる心電図の波形の全てを表す。グラフＧ６は、グラフＧ５の一部であってＴ波の領域の拡大図を示す。Ｔ波の領域は図７において領域Ａ１と示されている領域である。結果Ｇ７は、２つのガウス分布の統計量を示す。結果Ｇ７の各値は２つのガウス分布の統計量、を表す。２つのガウス分布の統計量とは、具体的には、内層側累積分布関数と外層側累積分布関数のガウス分布の平均値及び分散である。
　図８は、実施形態における信号解析装置１の解析結果の一例を示す第３の図である。より具体的には、図８は、動作がＴ延長３型である対象心臓の心電図の波形に対する信号解析装置１による解析結果の一例である。
　図８は、グラフＧ８、グラフＧ９及び結果Ｇ１０の３つの結果を示す。図８において、“内層側累積分布関数”は心筋内層に存在するチャネルの集合的なチャネル活動タイミング分布の候補を表す関数のフィッティング結果を示す。図８において、“外層側累積分布関数”は心筋外層に存在するチャネルの集合的なチャネル活動タイミング分布の候補を表す関数のフィッティング結果を示す。図８において、“体表面の電位”は、体表面の電位の時間変化を表す関数であって、心電図の波形である。図８の横軸は時刻を表し、縦軸は電位を表す。横軸と縦軸との単位はいずれも任意単位（arbitrary unit）である。
　グラフＧ８は、１回の拍動で生じる心電図の波形の全てを表す。グラフＧ９は、グラフＧ８の一部であってＴ波の領域の拡大図を示す。Ｔ波の領域は図８において領域Ａ２と示されている領域である。結果Ｇ１０は、２つのガウス分布の統計量を示す。結果Ｇ１０の各値は２つのガウス分布の統計量すなわち内層側累積分布関数と外層側累積分布関数のガウス分布の平均値及び分散を表す。
　図８の内層側累積分布関数及び外層側累積分布関数の形状は、カテーテルの挿入を行うことで動作がＴ延長３型である対象心臓の心筋外層の拍動によって生じた起電力の変化を直接測定した結果に略一致する。そのため、図８は、信号解析装置１によりカテーテルを挿入することなく、心電図だけからイオンチャネル活動情報が取得されることを示す。
　図９は、実施形態における信号解析装置１の解析結果の一例を示す第４の図である。より具体的には、図９は、動作がＱＴ延長１型である対象心臓の心電図の波形に対する信号解析装置１による解析結果の一例である。
　図９は、グラフＧ１１、グラフＧ１２及び結果Ｇ１３の３つの結果を示す。図９において、“内層側累積分布関数”は心筋内層に存在するチャネルの集合的なチャネル活動タイミング分布の候補を表す関数のフィッティング結果を示す。図９において、“外層側累積分布関数”は心筋外層に存在するチャネルの集合的なチャネル活動タイミング分布の候補を表す関数のフィッティング結果を示す。図９において、“体表面の電位”は、体表面の電位の時間変化を表す関数であって、心電図の波形である。図９の横軸は時刻を表し、縦軸は電位を表す。横軸と縦軸との単位はいずれも任意単位（arbitrary unit）である。
　グラフＧ１１は、１回の拍動で生じる心電図の波形の全てを表す。グラフＧ１２は、グラフＧ１１の一部であってＴ波の領域の拡大図を示す。Ｔ波の領域は図８において領域Ａ３と示されている領域である。結果Ｇ１３は、２つのガウス分布の統計量を示す。結果Ｇ１３の各値は２つのガウス分布の統計量を表す。２つのガウス分布の統計量とは、具体的には、内層側累積分布関数と外層側累積分布関数のガウス分布の平均値及び分散である。
　図９の内層側累積分布関数及び外層側累積分布関数の形状は、カテーテルの挿入を行うことで動作がＱＴ延長３型である対象心臓の心筋外層の拍動によって生じた起電力の変化を直接測定した結果に略一致する。そのため、図９は、信号解析装置１によりカテーテルを挿入することなく、心電図だけからイオンチャネル活動情報が取得されることを示す。信号解析装置１によって、心電図だけからイオンチャネル活動情報が取得されることを示す。
　なお、図９は、突然死と関連するチャネル電流特性の信号解析装置１による推定結果の一例でもある。
　図１０～図１２を用いて、心室性期外収縮の心電図に対しても信号解析装置１によって、心電図だけからイオンチャネル活動情報が取得されることを示す。図１０～図１２の横軸は時刻（秒）を表し、縦軸は電位（ｍＶ）を表す。
　図１０は、実施形態の信号解析装置１によって心室性期外収縮の心電図を解析した一例の第１の説明図である。図１１は、実施形態の信号解析装置１によって心室性期外収縮の心電図を解析した一例の第２の説明図である。図１２は、実施形態の信号解析装置１によって心室性期外収縮の心電図を解析した一例の第３の説明図である。
　より具体的には図１０は、体表面の心電位を示す。すなわち図１０は、心電図に記録された正常な一回の心拍と２連発の心室性期外収縮を示す。より具体的には図１１は、信号解析装置１によって解析した心室性期外収縮の脱分極の内層側累積分布関数と、外層側累積分布関数と、再分極相の内層側累積分布関数と、を示す。より具体的には図１２は、実際に計測された心電図の心室性期外収縮の波形の一例を示す。
　図１１は、脱分極の内層側累積分布関数は外層側累積分布関数に先行し、両者の分散は正常な心拍よりも大きく興奮の広がりが緩やかであることを示す。この分析結果は裾の広いＲ波の波形の特徴と一致する。
　図１１は、再分極相においては内層側累積分布関数が外層側累積分布関数よりも早く不活性化が開始していることを示しており、累積分布関数の平均の値としてそれが示されている。これらは再分極相の大きな陰性のＴ波の特徴と一致し、図１２に示すように、累積分布関数の差分により得た波形は実際に計測された心電図の心室性期外収縮の波形と略一致する。
　なお図１０～図１２の結果は、心筋の興奮の変更伝導、早期再分極もしくは再分極の遅延により、巨大波や陰性電位の生じる例に対して、信号解析装置１は対応することを示している。なお、図１０～図１２において、脱分極相の内層側累積分布関数の平均μは０．４６であり、分散σ２は１．４３である。また、図１０～図１２において、脱分極の外層側累積分布関数の平均μは－１であり、分散σ２は０．９９である。また、図１０～図１２において、再分極相の内層側累積分布関数の平均μは０．４６であり、分散σ２は１．４３である。また、図１０～図１２において、再分極相の外層側累積分布関数の平均μは－１であり、分散σ２は０．９９である。分散１＝σ２は、０．０８７秒間に相当する。
　図１３～図１５を用いて、ブルガダ症候群１型の脱分極期の対象心臓の心電図に対しても信号解析装置１によって、心電図だけからイオンチャネル活動情報が取得されることを示す。図１３～図１５の縦軸はミリボルト単位の電位を表す。
　図１３は、実施形態における信号解析装置１によってブルガダ症候群１型の脱分極期における対象心臓の心電図を解析した第１の説明図である。図１４は、実施形態における信号解析装置１によってブルガダ症候群１型の脱分極期における対象心臓の心電図を解析した第２の説明図である。図１５は、実施形態における信号解析装置１によってブルガダ症候群１型の脱分極期における対象心臓の心電図を解析した第３の説明図である。
　より具体的には図１３は、ブルガダ症候群の胸部第２誘導の心電図を示す。図１３において、内枠Ｗ１は脱分極相、内枠Ｗ２が再分極相の区分を示す。このことは図１４及び図１５についても同様である。すなわち、図１４及び図１５においても内枠Ｗ１は脱分極相を表し、内枠Ｗ２が再分極相を表す。
　より具体的には、図１４は信号解析装置１によって分析した、脱分極と、再分極相における、内層側累積分布関数と外層側累積分布関数とを示す。図１４に示す例においては、内層側累積分布関数の再分極相は脱分極相に続いて開始しており、早期再分極の特徴を示している。一方、図１４の例において外層側累積分布関数の電位振幅については脱分極相と再分極相とに差が認められる。また図１４は、外層側累積分布関数の脱分極相と再分極相とのギャップと異方性とを示している。このように、信号解析装置１は、ブルガダ症候群の対象心臓の心電図が示す波形の特徴である早期再分極や、脱分極と再分極との異方性を、内層側累積分布関数と外層側累積分布関数の脱分極相と再分極相の平均、分散、比率によって表すことができる。
　図１５は解析の結果と実測値との比較である。より具体的には図１５は、図１４の脱分極と、再分極相とにおける、内層側累積分布関数と外層側累積分布関数との差を示す。また図１５は心電図の実測値も示す。最後尾を除き、解析結果と実測値はほぼ一致している。最後尾とは、遅い時間の電位を意味する。
　なお、図１３～図１５において、脱分極相の内層側累積分布関数の平均μは１５であり、分散σ２は０．１５である。また、図１３～図１５において、脱分極の外層側累積分布関数の平均μは１４であり、分散σ２は０．２５である。また、図１３～図１５において、脱分極の内外層比は０．４５である。また、図１３～図１５において、再分極相の内層側累積分布関数の平均μは２５であり、分散σ２は０．２５である。また、図１３～図１５において、再分極相の外層側累積分布関数の平均μは２０であり、分散σ２は０．５である。また、図１３～図１５において、脱分極の内外層比は０．７５である。また、図１３～図１５において、分極再分極比は１．４である。
　図１６から図５４は、公開されている心電図データのライブラリhttps://physionet.org/about/database/を用いて、信号解析装置１による分析を行った結果を示す。図１６から図５４は、心電位に対する、信号解析装置１による解析の実行の結果得られた、内層側累積分布関数と外層側累積分布関数と、フィッティングの結果を示す。
　図１６から図５４はそれぞれ、実施形態における信号解析装置１によって心電図を解析した一例を示す図である。各図に記載の被決定点は図の左から順にそれぞれ、心電位のＱ点、Ｒ点、Ｓ点、Ｔ開始点、Ｔ終了点を表す。被決定点は、変曲点検知及びピーク検出アルゴリズムによって決定された。図１６～図５４の各図は、脱分極相（ＱＲＳ波）の区間と再分極相（Ｔ波）の区間とに対して得られた、各区間における内層側累積分布関数と外層側累積分布関数とを示す。図１６～図５４は、様々なＱＲＳ波とＴ波に対して、信号解析装置１は、内層側累積分布関数と外層側累積分布関数の平均と分散の調節により、略同一の形状を示し得ることを示す。図１６～図５４の下の図は、心電図の原波形とフィッティングの結果を表す。なお、図１６～図５４の各結果は３００Ｈｚのサンプリングの結果である。そのため、図１６～図５４の各図の横軸は原点が０秒を表し、値１が３．３３ミリ秒を表す。
　このように構成された信号解析装置１は、１又は複数のチャネル累積分布関数を用いて対象心臓の心電図の波形をフィッティングし、フィッティング結果の各チャネル累積分布関数の形状に基づきイオンチャネル活動情報を取得する。フィッティングに用いられる各チャネル累積分布関数は、チャネル活動タイミング分布の候補を示す情報である。チャネル活動タイミング分布はイオンチャネルの活動のタイミングの分布を表すものであり、互いに異なるチャネル累積分布関数はそれぞれ異なる種類のチャネル活動タイミング分布を表す。したがって、フィッティング結果の各チャネル累積分布関数の形状それぞれは、各チャネル累積分布関数が表すチャネル活動タイミング分布の情報を示すものである。このような情報は、従来の心電図の波形の解析によっては得られない。そのため、信号解析装置１は、心電図の波形から得られる情報を増やすことができる。
　（変形例）
　なお、定義域が重なるチャネル累積分布関数は必ずしも２つである必要は無く、３以上でもよい。このような場合、定義域の重なりを有する複数のチャネル累積分布関数の線形和であって重みに関する所定の条件を満たす線形和と心電図の波形との違いを最小にするフィッティングが行われてもよい。重みに関する所定の条件は、少なくとも１つが他の重みと異なる符号を有するという条件を少なくとも含めばどのような条件であってもよい。
　なお、心電図は、心起電力ベクトルに近い誘導の心電図が好ましい。心起電力ベクトルに近い誘導の心電図は、例えばＩＩ誘導や、Ｖ４誘導や、Ｖ５誘導などの心電位の立体的な情報を取得可能な心電図が好ましい。心電図が多チャンネルであるほど情報量が増えるので、心電図はチャンネル数が多いほど望ましい。
　なお、信号解析装置１は、ネットワークを介して通信可能に接続された複数台の情報処理装置を用いて実装されてもよい。この場合、信号解析装置１が備える各機能部は、複数の情報処理装置に分散して実装されてもよい。
　なお、心電図取得部１１０は生体情報取得部の一例である。
　なお、信号解析装置１の各機能の全て又は一部は、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）やＰＬＤ（Ｐｒｏｇｒａｍｍａｂｌｅ　Ｌｏｇｉｃ　Ｄｅｖｉｃｅ）やＦＰＧＡ（Ｆｉｅｌｄ　Ｐｒｏｇｒａｍｍａｂｌｅ　Ｇａｔｅ　Ａｒｒａｙ）等のハードウェアを用いて実現されてもよい。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。プログラムは、電気通信回線を介して送信されてもよい。
　以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
　１…信号解析装置、　１１…制御部、　１２…入力部、　１３…通信部、　１４…記憶部、　１５…出力部、　１１０…心電図取得部、　１２０…フィッティング情報取得部、　１３０…解析部、　１３１…フィッティング部、　１３２…イオンチャネル活動情報取得部、　１４０…記録部、　９１…プロセッサ、　９２…メモリ

Claims (6)

  1. 　解析対象の心臓の拍動に関する時系列の生体情報を取得する生体情報取得部と、
    　前記心臓におけるイオンチャネルの活動のタイミングの分布であるチャネル活動タイミング分布の候補を示す情報を分布候補情報として、前記生体情報と前記分布候補情報とに基づき、前記イオンチャネルの活動の様子を示す情報を取得する解析部と、
    　を備える信号解析装置。
  2. 　前記分布候補情報は、チャネル活動タイミング分布の候補を表す関数を確率密度関数とする累積分布関数によって表現される、
    　請求項１に記載の信号解析装置。
  3. 　前記解析部は、前記生体情報に対して１又は複数の前記累積分布関数をフィッティングさせた結果に基づき前記イオンチャネルの活動の様子を示す情報を取得する、
    　請求項２に記載の信号解析装置。
  4. 　前記フィッティングは、定義域の重なりを有する複数の前記累積分布関数の線形和であって少なくとも１つが他の重みと異なる符号を有するという条件を少なくとも含む重みに関する所定の条件を満たす線形和と前記生体情報との違いを最小にするフィッティングである、
    　請求項３に記載の信号解析装置。
  5. 　解析対象の心臓の拍動に関する時系列の生体情報を取得する生体情報取得ステップと、
    　前記心臓におけるイオンチャネルの活動のタイミングの分布であるチャネル活動タイミング分布の候補を示す情報を分布候補情報として、前記生体情報と前記分布候補情報とに基づき、前記イオンチャネルの活動の様子を示す情報を取得する解析ステップと、
    　を有する信号解析方法。
  6. 　請求項１から４のいずれか一項に記載の信号解析装置としてコンピュータを機能させるためのプログラム。
PCT/JP2021/000209 2021-01-06 2021-01-06 信号解析装置、信号解析方法及びプログラム WO2022149215A1 (ja)

Priority Applications (13)

Application Number Priority Date Filing Date Title
JP2022573841A JPWO2022149215A1 (ja) 2021-01-06 2021-01-06
US18/267,688 US20240050018A1 (en) 2021-01-06 2021-01-06 Signal analyzing apparatus, signal analyzing method and program
PCT/JP2021/000209 WO2022149215A1 (ja) 2021-01-06 2021-01-06 信号解析装置、信号解析方法及びプログラム
EP21917443.0A EP4275606A1 (en) 2021-01-06 2021-01-06 Signal analysis device, signal analysis method, and program
PCT/JP2021/033138 WO2022149307A1 (ja) 2021-01-06 2021-09-09 信号解析装置、信号解析方法及びプログラム
US18/039,553 US20240016433A1 (en) 2021-01-06 2021-12-02 Signal analyzing apparatus, signal analyzing method and program
EP21917602.1A EP4238498A1 (en) 2021-01-06 2021-12-02 Signal analysis device, signal analysis method, and program
JP2022573946A JP7457267B2 (ja) 2021-01-06 2021-12-02 信号解析装置、信号解析方法及びプログラム
EP21917603.9A EP4241691A1 (en) 2021-01-06 2021-12-02 Signal analysis device, signal analysis method, and program
US18/265,663 US20240057924A1 (en) 2021-01-06 2021-12-02 Signal analyzing apparatus, signal analyzing method and program
JP2022573945A JP7457266B2 (ja) 2021-01-06 2021-12-02 信号解析装置、信号解析方法及びプログラム
PCT/JP2021/044316 WO2022149381A1 (ja) 2021-01-06 2021-12-02 信号解析装置、信号解析方法及びプログラム
PCT/JP2021/044327 WO2022149382A1 (ja) 2021-01-06 2021-12-02 信号解析装置、信号解析方法及びプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/000209 WO2022149215A1 (ja) 2021-01-06 2021-01-06 信号解析装置、信号解析方法及びプログラム

Publications (1)

Publication Number Publication Date
WO2022149215A1 true WO2022149215A1 (ja) 2022-07-14

Family

ID=82357851

Family Applications (2)

Application Number Title Priority Date Filing Date
PCT/JP2021/000209 WO2022149215A1 (ja) 2021-01-06 2021-01-06 信号解析装置、信号解析方法及びプログラム
PCT/JP2021/033138 WO2022149307A1 (ja) 2021-01-06 2021-09-09 信号解析装置、信号解析方法及びプログラム

Family Applications After (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/033138 WO2022149307A1 (ja) 2021-01-06 2021-09-09 信号解析装置、信号解析方法及びプログラム

Country Status (4)

Country Link
US (1) US20240050018A1 (ja)
EP (1) EP4275606A1 (ja)
JP (1) JPWO2022149215A1 (ja)
WO (2) WO2022149215A1 (ja)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018504199A (ja) * 2015-01-09 2018-02-15 延世大学校 産学協力団Ｉｎｄｕｓｔｒｙ−Ａｃａｄｅｍｉｃ Ｃｏｏｐｅｒａｔｉｏｎ Ｆｏｕｎｄａｔｉｏｎ，Ｙｏｎｓｅｉ Ｕｎｉｖｅｒｓｉｔｙ 不整脈治療剤の効果評価システムおよび方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3057500B1 (en) * 2013-10-14 2020-12-23 Kaoskey Pty Limited Decomposition of non-stationary signals into functional components

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018504199A (ja) * 2015-01-09 2018-02-15 延世大学校 産学協力団Ｉｎｄｕｓｔｒｙ−Ａｃａｄｅｍｉｃ Ｃｏｏｐｅｒａｔｉｏｎ Ｆｏｕｎｄａｔｉｏｎ，Ｙｏｎｓｅｉ Ｕｎｉｖｅｒｓｉｔｙ 不整脈治療剤の効果評価システムおよび方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
HIROSHI TANAKA: "Methodology in Electrocardiographic Inverse Problems", MEDICAL ELECTRONICS AND BIOTECHNOLOGY, vol. 23, no. 3, 1985, pages 147 - 158
TANAKA HIROSHI: "On the Inverse Solution of Electrocardiology", IYO DENSHI TO SEITAI KOGAKU. JAPANESE JOURNAL OF MEDICAL ELECTRONICS AND BIOLOGICAL ENGINEERING, vol. 23, no. 3, 30 June 1985 (1985-06-30), JP , pages 147 - 158, XP009538878, ISSN: 0021-3292, DOI: 10.11239/jsmbe1963.23.147 *
TANAKA, YOSHIFUMI: "The ECG Recognition and Treatment of Antiarrhythmic Drugs", THE JOURNAL OF JAPAN SOCIETY FOR CLINICAL ANESTHESIA, vol. 30, no. 4, 28 October 2010 (2010-10-28), JP , pages 506 - 516, XP009538877, ISSN: 0285-4945, DOI: 10.2199/jjsca.30.506 *
WATARU SHIMIZU: "The role of M cells in the origin of T waves", JPN. J. ELECTROCARDIOLOGY, vol. 21, no. 2, 2001, pages 101 - 108
YOSHIFUMI TANAKASHINJI NODOYASUO YAMAZAKIARISA NAKASONE: "Program for constructing body surface electrocardiogram from cardiac action potential waveform", ANAESTHESIA/INTENSIVE CARE AND TECHNOLOGY, vol. 2011, no. 1, 2011, pages 91 - 95

Also Published As

Publication number Publication date
WO2022149307A1 (ja) 2022-07-14
EP4275606A1 (en) 2023-11-15
US20240050018A1 (en) 2024-02-15
JPWO2022149215A1 (ja) 2022-07-14

Similar Documents

Publication Publication Date Title
EP2663228B1 (en) Electroanatomical mapping
JP4759115B2 (ja) 心電図信号の交代を数量化するシステムおよび方法
JP4741168B2 (ja) 異なるセンサからの測定の組み合わせ
US8233972B2 (en) System for cardiac arrhythmia detection and characterization
US11147495B2 (en) Noninvasive method and system for estimating mammalian cardiac chamber size and mechanical function
US8923958B2 (en) System and method for evaluating an electrophysiological signal
KR101910982B1 (ko) 개인화된 생체 신호 패턴을 이용한 생체 신호의 동잡음 제거 방법 및 장치
JP2007500550A (ja) 心臓再同期療法の最適化方法
CN112543617A (zh) 用于量化和去除生物物理信号中的异步噪声的方法和系统
Sassi et al. An estimate of the dispersion of repolarization times based on a biophysical model of the ECG
US20120330557A1 (en) System for Cardiac Condition Analysis Based on Cardiac Operation Patterns
Baumert et al. Quantitative-electrogram-based methods for guiding catheter ablation in atrial fibrillation
Tueche et al. Embedded algorithm for QRS detection based on signal shape
US9549681B2 (en) Matrix-based patient signal analysis
Kotas et al. Averaging of nonlinearly aligned signal cycles for noise suppression
CN111528833B (zh) 一种心电信号的快速识别与处理方法及系统
WO2022149215A1 (ja) 信号解析装置、信号解析方法及びプログラム
WO2022149382A1 (ja) 信号解析装置、信号解析方法及びプログラム
KR101498581B1 (ko) 비침습적 심방신호 추정 시스템 및 방법
WO2023248308A1 (ja) 学習装置、情報提供装置、学習方法、情報提供方法及びプログラム
de Meireles ECG denoising based on adaptive signal processing technique
WO2023243090A1 (ja) 信号合成装置、信号合成方法及びプログラム
Mayapur Detection and classification of heart defects
KR20200055582A (ko) 심전도 신호 진단 파라미터 추출 디바이스 및 추출 방법
CN116615144A (zh) 信号解析装置、信号解析方法和程序

Legal Events

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

Ref document number: 21917443

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2022573841

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 18267688

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2021917443

Country of ref document: EP

Effective date: 20230807