CN104027110B - 心房颤动解析装置、心房颤动解析系统及解析方法 - Google Patents
心房颤动解析装置、心房颤动解析系统及解析方法 Download PDFInfo
- Publication number
- CN104027110B CN104027110B CN201310739868.9A CN201310739868A CN104027110B CN 104027110 B CN104027110 B CN 104027110B CN 201310739868 A CN201310739868 A CN 201310739868A CN 104027110 B CN104027110 B CN 104027110B
- Authority
- CN
- China
- Prior art keywords
- auricular fibrillation
- activity
- user
- pulse wave
- coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
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
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/361—Detecting fibrillation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1118—Determining activity level
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
Abstract
本发明涉及心房颤动解析装置、心房颤动解析系统以及心房颤动解析方法。其中,心房颤动解析装置,其特征在于,具备:取得部,取得表示心电或者脉波的检测结果的时间变化的检测波形信号;RR间隔算出部,根据通过所述取得的检测波形信号的频率解析得到的各帧的频谱,对应于每帧算出相当于平均RR间隔的参数,并算出表示该参数的时间变化的RR波形信号;功率算出部,算出表示所述RR波形信号的频谱中预先确定的频带的功率的时间变化;变动系数算出部,算出所述平均RR间隔的变动系数;解析部,根据所述算出的功率和变动系数的组合,对应于每一个时刻解析有无心房颤动;以及测量部,测量用户的活动量。
Description
技术领域
本发明涉及解析心房颤动的技术。
背景技术
在有关心脏疾病的医疗领域中具有解析心房颤动的技术。在专利文献1中,公开了一种测量由每一拍的心电图得到的RR间隔,根据其标准偏差和度数分布解析心房颤动的技术。在非专利文献1中,记载有RR间隔不规则,如果进行心房颤动的心拍的频率解析则存在1/fβ成分,由于其波动而成为白色杂音状的情况。
在先技术文献
专利文献
专利文献1:日本专利特开2009-89883号公报
非专利文献
非专利文献1:Hayano J,Yamasaki F,Sakata S,Okada A,Mukai S,Fujinami T“Spectral characteristics of ventricular response to atrial fibrillation.”Am.J.Physiol.1997;273:H2811-H2816
在上述专利文献1和非专利文献1中,为了进行准确的心房颤动的解析,需要准确地测量每一拍的RR间隔。在测量该RR间隔时,虽然也可以根据测定心电得到的心电图的波形信号进行测量,但也可以根据测定脉波得到的脉波信号进行测量。
可是,在测定脉波时,由于在测定中多数情况下受验者能够自由移动,因而脉波信号中也容易包含人体运动噪声的影响。在测定心电的情况下,虽然与测定脉波的情况相比具有程度的差异,但是心电图的波形信号中有时也包含人体运动噪声的影响。这样在受到人体运动噪声的影响时,准确地测量每一拍的RR间隔是非常困难的。
因此,在如专利文献1以及非专利文献1中公开的技术那样,将测量每一拍的准确的RR间隔作为前提时,不能够使用含有人体运动噪声的影响的信号进行心房颤动解析。另外,已知心房颤动的病症有时也受到患者的活动量的影响,但活动量与心房颤动的关系也存在个体差异而无法充分掌握。
发明内容
鉴于上述情况,本发明的目的之一在于,即使在脉波信号和心电图的波形信号等可测量RR间隔的信号中含有人体运动噪声的影响,也根据该信号解析心房颤动。另外,本发明的目的之一在于,对这样解析后的心房颤动和患者的活动量进行比较。
为解决上述课题,本发明提供一种心房颤动解析装置,具有:取得部,取得表示心电或者脉波的检测结果的时间变化的检测波形信号;波形存储部,存储所述取得的检测波形信号;RR间隔算出部,根据通过所述取得的检测波形信号的频率解析得到的各单位时间的频谱,对应于每一个该单位时间算出相当于平均RR间隔的参数;功率算出部,算出表示由所述RR间隔算出部算出的平均RR间隔的时间变化的RR波形信号的频谱中预先确定的频带的功率;变动系数算出部,算出所述RR波形信号中的所述平均RR间隔的变动系数;解析部,根据所述功率和所述变动系数的组合,对应于每一个时刻解析有无心房颤动;测量部,测量用户的活动量;以及第一存储部,使所述测量部测量的活动量和测量该活动量的时刻各自相对应地进行存储。
根据该心房颤动解析装置,即使在脉波信号和心电图的波形信号等可测量RR间隔的信号中含有人体运动噪声的影响,也能够根据该信号解析心房颤动,并能够对心房颤动和用户的活动量进行比较。
在优选的方式中,可以具有第二存储部,所述第二存储部对应于所述用户的行动的每一个时刻,存储表示该行动的行动数据。
根据该心房颤动解析装置,能够对心房颤动的发病时刻和用户的行动的时刻进行比较。
另外,在优选的方式中,也可以具有:范围存储部,存储表示所述用户的活动量能够变动的范围的范围信息;以及转换部,根据由所述范围存储部存储的范围信息,将由所述测量部测量的活动量转换为表示该用户的负荷的负荷量。
根据该心房颤动解析装置,当对心房颤动和用户的活动量进行比较时,能够反映用户的个体差异。
另外,在优选的方式中,可以具有第三存储部,所述第三存储部对应于所述用户自觉的所述心房颤动的每一个自觉的时刻,存储表示所述心房颤动的自觉症状数据。
根据该心房颤动解析装置,能够对心房颤动的发病时刻和用户自觉症状的时刻进行比较。
另外,在优选的方式中,可以具有:被所述用户操作的操作部;以及通知部,当所述测量部测量的活动量超过预先确定的阈值时,将该活动量超过该阈值的内容通知给该用户,以便所述用户对所述操作部进行操作,所述第二存储部在所述用户操作所述操作部时存储所述行动数据。
根据该心房颤动解析装置,能够根据用户的活动量,促使该用户进行用于存储行动数据的操作。
另外,在优选的方式中,可以具有评价部,所述评价部在由所述解析部解析为存在的所述心房颤动的发病期间,评价由所述测量部测量的活动量。
根据该心房颤动解析装置,能够进行对应于有无心房颤动的用户的活动量的评价。
另外,本发明提供一种心房颤动解析系统,具有上述心房颤动解析装置;加工装置,取得所述测量部测量的活动量,并根据由所述解析部进行解析的结果对表示该活动量的数据进行加工;以及提示装置,提示由所述加工装置加工的数据。
根据该心房颤动解析系统,即使在脉波信号和心电图的波形信号等可测量平均RR间隔的信号中含有人体运动噪声的影响,也能够根据该信号解析心房颤动,并能够对心房颤动和用户的活动量进行比较。
另外,本发明提供一种心房颤动解析方法,具有:取得表示心电或者脉波的检测结果的时间变化的检测波形信号的步骤;存储所述取得的检测波形信号的步骤;根据通过所述取得的检测波形信号的频率解析得到的各单位时间的频谱,对应于每一个该单位时间算出相当于平均RR间隔的参数的步骤;算出表示所述算出的平均RR间隔的时间变化的RR波形信号的频谱中预先确定的频带的功率的步骤;算出所述RR波形信号中的所述平均RR间隔的变动系数的步骤;根据所述功率和所述变动系数的组合,对应于每一个时刻解析有无心房颤动的步骤;测量用户的活动量的步骤;以及使所述测量的活动量和测量该活动量的时刻各自相对应地进行存储的步骤。
根据该心房颤动解析方法,即使在脉波信号和心电图的波形信号等可测量平均RR间隔的信号中含有人体运动噪声的影响,也能够根据该信号解析心房颤动,并能够对心房颤动和用户的活动量进行比较。
另外,本发明提供一种程序,用于使计算机执行以下步骤:取得表示心电或者脉波的检测结果的时间变化的检测波形信号的步骤;存储所述被取得的检测波形信号的步骤;根据通过所述被取得的检测波形信号的频率解析得到的各单位时间的频谱,对应于每一个该单位时间算出相当于平均RR间隔的参数的步骤;算出表示所述算出的平均RR间隔的时间变化的RR波形信号的频谱中预先确定的频带的功率的步骤;算出所述RR波形信号中的所述平均RR间隔的变动系数的步骤;根据所述功率和所述变动系数的组合,对应于每一个时刻解析有无心房颤动的步骤;测量用户的活动量的步骤;以及使所述测量的活动量和测量该活动量的时刻各自相对应地进行存储的步骤。
根据该程序,即使在脉波信号和心电图的波形信号等可测量平均RR间隔的信号中含有人体运动噪声的影响,也能够根据该信号解析心房颤动,并能够对心房颤动和用户的活动量进行比较。
附图说明
图1是用于说明实施方式中的心房颤动解析系统的结构的图。
图2是说明实施方式中的脉波测定装置的外观的图。
图3是说明实施方式中的脉波测定装置的结构的图。
图4是说明实施方式中的心房颤动解析装置的功能结构的图。
图5是说明进行检测波形信号的频率解析时的帧的图。
图6是说明RR间隔算出部的频率特性的图。
图7是说明进行RR波形信号的频率解析时的帧的图。
图8是表示由本实施方式进行的判定结果的图。
图9是心房颤动解析处理的流程图。
图10是说明采用了心电RR间隔的心房颤动的解析方法的图。
图11是说明采用了平均脉波RR间隔的心房颤动的解析方法的问题点的图。
图12是表示步骤A410中的心房颤动解析处理的细节的图。
图13是例示功率波形信号和变动系数信号的图。
图14是例示移动平均处理后的功率和变动系数的图。
图15是例示功率和变动系数的关系的图。
图16是例示图15的数据的聚类后的结果的图。
图17是表示变形例的心房颤动解析系统的结构的图。
图18是表示在可对比的形态下表现了数据的曲线图的一个示例的图。
符号说明
1,1a…脉波测定装置、10…装置主体、11…CPU、12…RAM、13…ROM、14…操作部、140…操作按钮、141…行动按钮、142…自觉症状数据、15…显示部、16…测量部、17…计时电路、171…振荡电路、18…A/D转换电路、19…放大电路、20…脉波检测部、30…电缆、40…固定垫、50…手腕套、100…心房颤动解析装置、101…通信部、111…噪声降低部、112…RR间隔算出部、113…功率算出部、114…解析部、115…显示控制部、116…变动系数算出部、117…转换部、118…评价部、119…通知部、121…检测波形信号存储区域、122…RR波形信号存储区域、123…功率波形信号存储区域、6…加工装置、7…提示装置、9,9a…心房颤动解析系统、A…患者、B…医生、C…加工中心、1000…手。
具体实施方式
实施方式
概要
现有技术中,根据心电图解析心房颤动时,采用了每一拍的RR间隔(称为“心电RR间隔”)。与此相对,在本实施方式中,根据脉波来解析心房颤动。与心电图不同,采用脉波时有时难以准确地确定每一拍的RR间隔。因此,在本实施方式中,采用表示在某一单位时间(以下,称为帧)内的RR间隔的平均的值(称为“平均脉波RR间隔”)。在此,首先,说明有关采用了心电RR间隔的心房颤动的解析方法,接着,说明采用平均脉波RR间隔时的问题点,最后,说明本实施方式中的解析方法的概要。
图10是说明采用了心电RR间隔的心房颤动的解析方法的图。图10是对表示心电RR间隔的变动的波形信号,将480秒作为1帧,在该帧中,在0.01Hz~0.2Hz的波段进行频率解析,对峰值频率和功率进行对数转换后表示的曲线图。功率是指谱密度。图10(a)表示心房颤动未发病时采用心电RR间隔的情况,图10(b)表示心房颤动发病时采用心电RR间隔的情况。图中的直线表示由绘制的数据得到的线性回归线。如果根据这些曲线图算出线性回归线的斜率β和相关系数γ,则得到如下所示结果。
在图10(a)所示的心房颤动未发病的情况下,γ=-0.72,β=-1.29。另外,在图10(b)所示的心房颤动发病的情况下,γ=-0.07,β=-0.13。由此可知,心房颤动发病的情况下,相关消失而出现白色杂音状,而且斜率β接近于“0”。这样,在采用心电RR间隔时,能够根据峰值频率和功率的绘图中的线性回归线的斜率β和相关系数γ来判定有无心房颤动。
图11是说明采用了平均脉波RR间隔的心房颤动的解析方法的问题点的图。图11是对表示平均脉波RR间隔的变动的波形信号,将480秒作为1帧,在该帧中,在0.01Hz~0.2Hz的波段进行频率解析,对峰值频率和功率进行对数转换后表示的曲线图。图11(a)表示心房颤动未发病的示例,图11(b)表示心房颤动发病的示例。图中的直线表示由绘制的数据得到的线性回归线。
在图11(a)所示的心房颤动未发病的情况下,γ=-0.68,β=-1.40。另外,在图11(b)所示的心房颤动发病的情况下,γ=-0.41,β=-1.02。这样,在采用平均脉波RR间隔时,如图11所示,不存在由有无心房颤动的发病引起的γ、β的显著性差异,通过与采用心电RR间隔的情况相同的方法难以判定有无心房颤动。
在此,再次比较图10(a)和(b)则可知,心房颤动发病的情况下,在高频带一侧功率增大。例如图10中,如果对0.2Hz附近的频带比较功率,则心房颤动未发病时的功率为“1.59”,心房颤动发病时的功率为“4.97”。心房颤动发病时,与未发病时相比,该频带的功率增加至数倍而出现显著性差异。
在采用平均脉波RR间隔时也观测到该功率的增加。在图11中,如果对0.2Hz附近的频带比较功率,则心房颤动未发病时的功率为“0.05”,心房颤动发病时的功率为“0.30”。这样,即使采用平均脉波RR间隔,心房颤动发病时,与未发病时相比,该频带的功率也增加至数倍而出现显著性差异。在本实施方式中,将该功率的增大用作一个指标来判定有无心房颤动。
在本实施方式中,采用变动系数作为有无心房颤动的另一个指标。变动系数是表示相对于平均脉波RR间隔的平均的偏差程度的参数。如果心房颤动发病,则引起RR间隔的不整齐。即,每一拍的时间间隔变得不规则。对于平均脉波RR间隔也是同样的,能够将不整齐状态(相对于平均的偏差)作为心房颤动的指标。在本实施方式中,将功率和变动系数作为指标来判定有无心房颤动。下面,将详细地说明本实施方式的装置结构以及动作。
心房颤动解析系统9的结构
图1是用于说明实施方式中的心房颤动解析系统9的结构的图。患者A是接受心房颤动的诊断的人。医生B是患者A的主治医生,为了诊断心房颤动将脉波测定装置1出借给患者A。患者A是佩戴有出借的脉波测定装置1,利用该脉波测定装置1的用户。加工中心C是接收从患者A佩戴的脉波测定装置1得到的数据的输出,按照规定的程序加工该数据的设备。医生B接受在加工中心C中被加工的数据的提示而进行有关患者A的心房颤动的诊断。
以下说明的实施方式是脉波测定装置1具有加工中心C中的加工功能的情况。
脉波测定装置1的结构
图2是说明实施方式中的脉波测定装置1的外观的图。如图2(a)所示,本发明的实施方式中的脉波测定装置1具有如手表那样被佩戴于作为用户的检测对象的手1000中的手腕部分(腕)上的装置主体10,以及安装于检测部位而检测脉波的脉波检测部20。装置主体10和脉波检测部20由电缆30连接。电缆30将从脉波检测部20输出的脉波信号(以下,称为检测波形信号L)供给装置主体10,另一方面,将来自装置主体10的电能供给脉波检测部20。
在装置主体10上安装有腕套50。装置主体10通过将腕套50缠绕在用户的手腕上而被安装在手腕上。在装置主体10上设有操作部14和显示部15。操作部14具有用户用于向脉波测定装置1输入功能选择的指示等的按钮开关等操作件,通过这些操作件接受用户的操作。操作部14在装置主体10的侧面左右分别各具有两个按钮开关。在这些按钮中,通过腕套50将装置主体10安装于用户的手腕部分时,设于小指侧的按钮开关是行动按钮141,设于拇指侧的按钮开关是自觉症状按钮142。另外,在这些按钮中,设于用户的胳膊肘(未图示)侧的两个按钮是具有其他功能的操作按钮140。操作部14也可以含有设于显示部15上的触摸传感器等。显示部15是液晶显示器、有机EL显示器等显示设备。
如图2(b)所示,在该示例中,安装有脉波检测部20的检测部位是手1000中的从食指的指根到第二指关节之间的一部分,但只要是能够检测脉波的部位则可以是任何部位。脉波检测部20通过由固定带40固定而安装于检测部位。此时,固定带40是覆盖脉波检测部20的状态,并且是以来自固定带40的外侧的光不会到达脉波检测部20的光接收部的方式进行遮光的结构。
脉波检测部20如下所示检测脉波,并输出表示检测结果的检测波形信号L。脉波检测部20具有发光部(例如绿色LED(Light Emitting Diode:发光二极管))和光接收部。脉波检测部20中由发光部照射对应于经由电缆30而从装置主体10供给的电能的光。脉波检测部20通过光接收部接收从发光部发出的光中被毛细血管中的血红蛋白反射的光,将对应于光接收电平的信号作为检测波形信号L后经由电缆30供给装置主体10。
图3是说明实施方式中的脉波测定装置1的结构的图。脉波测定装置1具有装置主体10和脉波检测部20,其中装置主体10具有:CPU(Central Processing Unit:中央处理器)11、RAM(Random Access Memory:随机存储器)12、ROM(Read Only Memory:只读存储器)13、操作部14、显示部15、测量部16、计时电路17、A/D转换电路18、放大电路19。除放大电路19和脉波检测部20以外的各结构都经由总线连接。
CPU11按照存储于ROM13的控制程序,进行各部控制以及数据传送。RAM12临时存储检测波形信号L等生物体信息以及CPU11中的控制程序在执行中产生的各种数据。CPU11通过执行控制程序,实现心房颤动解析功能,使脉波测定装置1作为心房颤动解析装置而起作用。此外,CPU11也可以通过执行控制程序实现心房颤动解析功能以外的各种功能。这些功能例如可以通过用户操作操作部14而实现。
如上所述,操作部14具有用于将用户的指示输入脉波测定装置1的按钮开关等。如果用户操作操作部14,则向CPU11输出表示操作内容的操作信号。操作部14具有的按钮开关中,行动按钮141是用户进行预先决定的行动时按压的操作件。另外,自觉症状按钮142是用户感觉到心房颤动的自觉症状时按压的操作件。
如上所述,显示部15具有液晶显示器、有机EL显示器等显示设备,显示内容由CPU11控制。该显示内容例如为表示时刻显示、各种菜单画面、脉波测定结果、心房颤动的解析结果等的各种图像。
测量部16是安装于装置主体10的壳体内部、检测作为该安装部位(本实施方式中为手腕)的动作的人体运动、并输出对应于检测结果的信号的传感器。该信号由于表示人体运动的程度,因而根据该信号可以求出用户的活动量。即,测量部16相当于测量用户的活动量的测量部。此外,活动量是指表示生物体活动程度的量,例如,随着该生物体的活动而消耗的热量等。
具体而言,测量部16具备静电电容式加速度传感器或者压电电阻式加速度传感器等,分别检测3轴方向的加速度。然后,测量部16根据例如检测出的上述3轴方向的加速度而算出上述安装部位的加速度的绝对值和用户的活动量,并输出至CPU11。
计时电路17具备振荡电路171。振荡电路171将作为控制基础的时钟信号供给CPU11。计时电路17通过CPU11的控制测量时间。
在行动按钮141或者自觉症状按钮142被按压时,CPU11参照由计时电路17测量的时间,将这些定时(timing)分别存储于RAM12。
放大电路19将经由电缆30从脉波检测部20供给的检测波形信号L放大。放大时的增益通过CPU11的控制来设定。
A/D转换电路18将在放大电路19中放大的模拟信号的检测波形信号L转换为数字信号。在该示例中,采样频率为100Hz,是与根据脉波得到的RR间隔相比十分高的频率。另外,在该例中,量子化以10比特进行。此外,关于采样频率、量子化比特,也可以根据需要的精度确定为不同的值。
接着,说明由CPU11实现的心房颤动解析装置的功能结构(心房颤动解析功能以及存储功能)。
功能结构
图4是说明实施方式中的心房颤动解析装置100的功能结构的图。心房颤动解析装置100具有噪声降低部111、RR间隔算出部112、功率算出部113、解析部114、显示控制部115、变动系数算出部116、转换部117、评价部118、通知部119,并通过作为各种数据的存储区域的检测波形信号存储区域121、RR波形信号存储区域122、功率波形信号存储区域123、变动系数波形存储区域124、活动量存储区域125、行动数据存储区域126、自觉症状数据存储区域127、范围存储区域128以及阈值存储区域129的各功能结构来实现。
检测波形信号存储区域121是设于存储有由A/D转换电路18转换为数字信号的检测波形信号L的RAM12上的区域。
噪声降低部111进行从存储于检测波形信号存储区域121的检测波形信号L中降低相当于RR间隔的频带以外的人体运动噪声成分的滤波器处理并输出。作为滤波处理,是由例如高通滤波器、带通滤波器、自适应滤波器等进行的处理。在噪声降低部111中降低人体运动噪声成分后的检测波形信号L,也可以暂且存储在RAM12中。检测波形信号存储区域121以及噪声降低部111作为取得在RR间隔算出部112中用于频率解析的检测波形信号L的取得部而起作用。
此外,在该处理中,虽然人体运动噪声成分被降低,其影响从检测波形信号L降低,但是在作为背景技术示出的技术(专利文献1、非专利文献1)中,越能够进行精密的心房颤动的解析,越不能够测量准确的RR间隔。
RR间隔算出部112对在噪声降低部111中降低人体运动噪声成分后的检测波形信号L,按照每次采样切出帧,通过短时间内的频率解析(STFT(Short-Time Fouriertransform:短时傅里叶变换)解析)算出频谱。然后,RR间隔算出部112根据算出的频谱,按照每帧算出相当于RR间隔的参数,并将表示该参数的时间变化的RR波形信号FRR存储于设于RAM12上的区域RR波形信号存储区域122。此外,RR波形信号FRR是表示该参数随时间变化的数据的集合。
算出的参数在该例中是表示帧内的RR间隔的平均的值(平均脉波RR间隔),例如是成为频谱的最大峰值的频率。因此,RR波形信号FRR表示平均脉波RR间隔的时间变化。通过在RR间隔算出部112中的处理,即使在噪声降低部111中人体运动噪声未被完全除去,也能够大幅度地降低RR波形信号FRR中含有的人体运动噪声的影响。
图5是说明进行检测波形信号L的频率解析时的帧的图。图5所示的波形是检测波形信号L的波形的示例。如图5所示,各帧的时间在该示例中是4秒,按照每一秒采样并进行频率解析。即,各帧逐秒错开设定,与下一帧重叠3秒钟。这样,由于采样定时、帧被设定,因而平均脉波RR间隔是RR间隔在4秒内的平均值,RR波形信号FRR表示每一秒的平均脉波RR间隔的变化。
图6是说明RR间隔算出部112中的频率特性的图。在RR间隔算出部112中,以如上所述设定的帧进行频率解析与移动平均处理中的频率特性重叠是等同的。图6所示的频率特性成为如下所示的频率特性:在相当于帧的时间4秒的0.25Hz及其整数倍的频率处产生波谷,作为连结波峰的顶点那样的整体的倾向,频率越高则电平越低,即具有负的斜率。帧的时间越长则该斜率越陡。反之,帧的时间越短则斜率越接近于“0”,但是导致检测波形信号L中的人体运动噪声成分的残留增多。因此,帧时间最好为1秒以上5秒以下,优选为2秒以上4秒以下。
功率算出部113对存储于RR波形信号存储区域122的RR波形信号FRR进行短时间内的频率解析(STFT解析),根据所得到的频谱算出一部分频带(以下,称为算出频带)的功率(以下,称为波段功率)。功率算出部113将表示算出的波段功率的时间变化的功率波形信号Pa存储于设于RAM12的区域的功率波形信号存储区域123。此外,功率波形信号Pa是表示波段功率随时间变化的数据的集合。
图7是说明进行RR波形信号FRR的频率解析时的帧的图。图7所示的波形是RR波形信号FRR的波形的示例。如图7所示,各帧的时间在该示例中是120秒,每60秒进行采样并进行频率解析。即,各帧各错开60秒设定,与下一个帧重叠60秒。
另外,在功率算出部113中波段功率被算出的上述算出频带预先被确定,在该示例中,是0.25Hz至0.5Hz的波段。这被确定为图6所示的频率特性的2个波谷(0.25Hz、0.5Hz的波谷)之间。这由于波谷部分的功率受到抑制而对于有无心房颤动的判定几乎没有帮助,因此以有助于判定有无心房颤动的部分为中心而确定算出频带。即,如0.3Hz至0.45Hz的波段这样,可以以进一步除去频率特性的波谷的部分,仅包括波峰的部分而将其范围变窄的方式设定算出频带。
在此,在该例中,算出频带的最低频率(下限)和最高频率(上限)根据RR间隔算出部112中的频率特性、即在RR间隔算出部112的频率解析中所使用的帧的时间而决定。另一方面,上下限频率中的任意一方或者双方的频率也可以不必一定根据帧的时间而决定。
如图10、图11所示,算出频带的最低频率最好为功率的变化明确的0.1Hz以上,优选为0.2Hz以上。此时,如上所述,最低频率进一步优选为RR间隔算出部112的频率解析中所使用的帧的时间的倒数以上。
另外,考虑到由奈奎斯特频率所产生的影响,算出频带的最高频率优选为RR间隔算出部112中的频率解析的采样频率的1/2以下。此时,如上所述,最高频率进一步优选为RR间隔算出部112的频率解析中所使用的帧的时间的倒数的两倍以下。
返回图4继续进行说明。变动系数算出部116根据存储于RR波形信号存储区域122的RR波形信号FRR(平均脉波RR间隔),按照下式(1)算出变动系数CVRR。
CVRR=σRR/aveRR…(1)
此外,σRR和aveRR分别表示每帧期间中的平均脉波RR间隔的标准偏差和平均值。即,变动系数CVRR是表示相对于平均的偏差程度的参数。变动系数算出部116将算出的变动系数CVRR存储于变动系数波形存储区域124。由于按照每帧算出变动系数CVRR,因而在变动系数波形存储区域124存储有表示变动系数CVRR的时间变化的信号(称为“变动系数信号Sc”)。此外,变动系数信号Sc是表示变动系数CVRR随时间变化的数据的集合。
解析部114根据存储于功率波形信号存储区域123的功率波形信号Pa以及存储于变动系数波形存储区域124的变动系数信号Sc,按照每一时刻解析有无心房颤动。具体而言,对这些信号判定是否满足特定的判定条件,并输出对应于判定结果的信息。有关特定的判定条件在后面说明。
如果判定为是心房颤动,则解析部114将表示该判定结果的信息输出至显示控制部115。从解析部114输出的信息只要是例如表示是心房颤动的判定的信息等关于有无心房颤动的信息即可。显示控制部115根据从解析部114输出的信息,控制显示部15的显示内容,使表示判定为是心房颤动的图像显示。用户通过观看该显示内容,就能够确认是否是心房颤动。此外,该显示内容既可以是实时地表示心房颤动的判定结果的显示,也可以是表示判定为是心房颤动的期间的显示。
以上是心房颤动解析装置100的功能结构的说明。接下来,使用图9来说明有关心房颤动解析装置100的解析动作(心房颤动解析处理)。
心房颤动解析处理
图9是说明实施方式中的心房颤动解析处理的流程图。首先,如果用户操作操作部14输入开始心房颤动的解析处理的指示,则CPU11开始图9所示的流程。CPU11判定是否用户操作操作部14输入结束解析处理的指示(步骤S110)。在结束解析处理的指示被输入时(步骤S110:“是”),CPU11结束心房颤动的解析处理。
在结束解析处理的指示未被输入时(步骤S110:“否”),CPU11在脉波检测部20中检测脉波并测定检测波形信号L(步骤S120),由噪声降低部111进行人体运动噪声降低处理(步骤S130)。此时,CPU11虽然将检测波形信号L存储于RAM12的检测波形信号存储区域121,但也可以存储已进行人体运动噪声降低处理的检测波形信号L。
CPU11判定已进行人体运动噪声降低处理的波形信号是否已在RAM12中积蓄一帧(步骤S140)。在尚未积蓄一帧时(步骤S140:“否”),CPU11返回步骤S110继续进行处理。反之,在已经积蓄一帧时(步骤S140:“是”),CPU11通过RR间隔算出部112算出平均脉波RR间隔(步骤S210)。
CPU11将通过RR间隔算出部112算出的平均脉波RR间隔存储于RR波形信号存储区域122(步骤S220)。存储于该存储区域的平均脉波RR间隔的时间变化成为RR波形信号FRR。
CPU11判定存储于RR波形信号存储区域122的RR波形信号FRR是否已积蓄一帧(步骤S230)。在尚未积蓄一帧时(步骤S230:“否”),CPU11返回步骤S110继续进行处理。反之,在已经积蓄一帧时(步骤S230:“是”),CPU11通过功率算出部113算出波段功率(步骤S310)。
CPU11将通过功率算出部113算出的波段功率存储于功率波形信号存储区域123(步骤S320)。存储于该存储区域的波段功率的时间变化成为功率波形信号Pa。
CPU11通过变动系数算出部116算出变动系数(步骤S330)。CPU11将通过变动系数算出部116算出的变动系数存储于变动系数波形存储区域124(步骤S340)。
CPU11参照存储的功率波形信号Pa和变动系数信号Sc,通过解析部114判定功率波形信号Pa和变动系数信号Sc是否满足规定的判定条件(步骤S410)。
图12是表示步骤S410中的心房颤动解析处理的细节的图。在步骤S500中,CPU11从RAM12读取功率波形信号Pa和变动系数信号Sc。
图13是例示功率波形信号Pa和变动系数信号Sc的图。图13示出从对某位患者历时24小时测量的脉波信号中得到的功率波形信号Pa[msec2]和变动系数信号Sc[%]。此外,该患者在测定期间中心房颤动发病。
再次参照图12。在步骤S501中,CPU11对功率波形信号Pa和变动系数信号Sc进行移动平均处理。为了将微细的变动(短时间内的变动)平均而分别对功率以及变动系数CVRR进行移动平均处理。在该示例中,使用20点的数据(即,从20分钟的脉波测定得到的数据)进行移动平均处理。
图14是例示移动平均处理后的功率波形信号Pa和变动系数信号Sc的图。通过移动平均处理,微细的变动被平均。以下,移动平均处理后的数据作为表示某一时刻的功率和变动系数的数据来处理。由于每60秒进行测定,因而在24小时的测定中可以得到1440点的数据。
再次参照图12。在步骤S502中,CPU11通过规定的算法(例如,作为聚类的方法而被广为人知的k-means法),将这些数据聚类为两个簇群。
图15是例示功率和变动系数的关系的图。纵轴表示功率[msec1],横轴表示变动系数CVRR[%]。如已经说明的,心房颤动发病时与平时相比,功率以及变动系数都显示出相对高的值。因此,可以认为图15的右上附近的绘图对应于心房颤动发病时测定的数据。在本实施方式中,通过聚类方法将数据划分为两个簇群,根据这两个簇群在变动系数-功率空间中的位置关系,解析有无心房颤动。
图16是例示图15的数据聚类后的结果的图。这样,通过聚类方法(在该例中为k-means法),数据被分割为两个簇群。将这两个簇群分别称为簇群C0和簇群C1。簇群C0是相对高变动系数且高功率的数据群,簇群C1是低变动系数且低功率的数据群。簇群C0的数据用圆圈(○)表示,簇群C1用三角(△)表示。另外,根据k-means法计算各簇群在变动系数-功率空间中的重心的坐标。图16同时图示出簇群C0以及C1的重心位置。
再次参照图12。在步骤S503中,CPU11判定簇群C0的重心位置是否以簇群C1的重心位置为基准而在规定的范围例如±30%的范围内。当判定为簇群C0的重心位置在距离簇群C1的重心位置的±30%的范围以外时(S503:“否”),CPU11将处理转移至步骤S504。当判定为簇群C0的重心位置在距离簇群C1的重心位置的±30%的范围内时(S503:“是”),CPU11将处理转移至步骤S505。
在步骤S504中,CPU11判定为簇群C0是心房颤动发病时的数据。
当簇群C0的重心位置在距离簇群C1的重心位置的±30%的范围内时,判定为数据不能够分离为两个簇群。此时作为可能性,可以考虑在整个测定期间心房颤动一直未发病的情况,以及在整个测定期间心房颤动持续发病的情况。此时,通过功率和变动系数的值解析有无心房颤动。在步骤S505以后,进行用于该解析的处理。
在步骤S505中,CPU11计算簇群C0的重心与簇群C1的重心的平均位置(以下称为“平均重心”)的坐标(以下称为“平均重心坐标”)。平均重心例如是簇群C0的重心与簇群C1的重心的简单平均(即簇群C0的重心与簇群C1的重心的中点)。或者,平均重心也可以是由簇群C0的重心与簇群C1的重心的数据点数得到的加权重心(即,所有测定点的重心)。
在步骤S506中,CPU11判定平均重心坐标是否在规定的范围内(例如,变动系数10.0以上且功率0.5以上)。当判定为平均重心坐标在规定的范围内时(S506:“是”),CPU11判定为在整个测定期间心房颤动一直在发病(步骤S507)。当判定为平均重心坐标不在规定的范围内时(S506:“否”),CPU11判定为在整个测定期间心房颤动一直未发病(步骤S508)。
再次参照图9。当判定为心房颤动一直未发病时(步骤S410:“否”),CPU11返回至步骤S110继续进行处理。反之,当判定为心房颤动发病时(步骤S410:“是”),CPU11通过显示控制部115使是心房颤动的判定结果显示在显示部15上(步骤S420),然后返回步骤S110继续进行处理。
此外,CPU11也可以不管步骤S140的判定如何,都反复执行步骤S110至步骤S140的处理。此时,CPU11也可以每次在步骤S140中为“是”时,与步骤S110至步骤S140的处理并行地执行步骤S210以后的处理。此时,在步骤S230中“否”的情况或者步骤S410中“否”的情况下,可以结束并行执行的步骤S210以后的处理。
以上是心房颤动解析处理的说明。
图8是表示基于本实施方式的解析结果的图。为了便于比较,同时图示出通过利用霍尔特式心电仪进行心电图解析而解析为心房颤动在发病的期间。图中粗线表示的部分是通过本实施方式解析为心房颤动在发病的期间,阴影部分是通过利用霍尔特式电仪进行心电图解析而解析为心房颤动在发病的期间。对于利用霍尔特式心电仪解析为心房颤动在发病的期间,在本实施方式中几乎没有被遗漏,也被解析为心房颤动在发病。在本实施方式中,虽然作为心房颤动而检测出一部分噪声,但在例如比规定的阈值短的期间被解析为心房颤动时,也可以将该解析结果放弃(即,在比阈值长的时间,连续地解析为心房持续颤动的情况下,采用是心房颤动这样的解析结果)。
这样,在本实施方式的脉波测定装置1中,通过测量平均脉波RR间隔取代测量每一拍的脉波RR间隔,能够减少人体运动噪声的影响,并进行心房颤动的解析。
存储处理
接下来,对基于心房颤动解析装置100的各种数据的存储处理进行说明。
如图4所示,测量部16向心房颤动解析装置100输出表示用户的活动量的信息。该输出既可以周期性地进行,也可以根据心房颤动解析装置100的要求进行。另外,该输出也可以在测量部16测量的活动量超过例如预先确定的阈值时等满足规定的条件时进行。由测量部16输出的活动量通过参照计时电路17测量的时间,与测量该活动量的时刻一起被存储于活动量存储区域125。
范围存储区域128存储有表示用户的活动量能够变动的范围的范围信息。用户的活动量由于该用户的年龄、体质和生活习惯等各种因素而存在个体差异。范围存储区域128存储的范围信息是表示佩戴有脉波测定装置1的用户的活动量的上限和下限的信息。活动量的下限是例如用户安静时的活动量,活动量的上限是例如用户通过使身体运动而感到难以承受的负荷时的活动量。
如果用户操作操作部14的操作按钮140等,则CPU11根据该操作而使用户的范围信息存储于范围存储区域128。此外,范围存储区域128可以事先使多个用户的范围信息与各用户的识别信息的每一个建立对应而存储。这种情况下,例如CPU11可以通过操作部14受理用户的操作并进行认证处理,识别该用户,从范围存储区域128提取出对应于该用户的范围信息。
转换部117根据存储于范围存储区域128的范围信息,将存储于活动量存储区域125的活动量转换为表示用户的负荷的负荷量。负荷量既可以例如用以多阶段划分由范围信息表示的活动量的范围的数值表示,也可以用由该活动消耗的消耗热量表示。消耗热量可以根据例如由范围信息表示的活动量的范围和对每个用户设定的基础代谢等的值来推算。
阈值存储区域129存储有与用户的活动量或者负荷量进行比较的阈值。在此,阈值存储区域129存储与用户的负荷量进行比较的阈值。在将转换部117转换而得到的用户的负荷量与存储于阈值存储区域129的该用户的阈值进行比较,用户的负荷量超过阈值时,通知部119将该内容通知给用户,以便用户进行对操作部14的操作。具体而言,通知部119如果判定为用户的负荷量超过上述的阈值,则向显示控制部115输出确定的信号。显示控制部115如果接受到该信号,则使例如催促用户按压行动按钮141的画面显示在显示部15上。
操作部14具有的行动按钮141如果被用户按压,则将分配给行动按钮141的固有信号输出至心房颤动解析装置100。心房颤动解析装置100如果接受到该信号,则将该信号作为“行动数据”存储于行动数据存储区域126。心房颤动解析装置100参照由计时电路17测量的时间,按照每一个行动按钮141被按压的时刻将行动数据存储于行动数据存储区域126。
操作部14具有的行动按钮142如果被用户按压,则将分配给行动按钮142的固有信号输出至心房颤动解析装置100。心房颤动解析装置100如果接受到该信号,则将该信号作为“自觉症状数据”存储于自觉症状数据存储区域127。心房颤动解析装置100参照由计时电路17测量的时间,按照每一个自觉症状按钮142被按压的时刻将自觉症状数据存储于自觉症状数据存储区域127。
评价部118按照每一个通过测量部16测量对应于由转换部117转换的负荷量的活动量的时刻,取得该负荷量。另外,评价部118按照每一个心房颤动产生的时刻或者消失的时刻,取得表示由解析部114解析的有无心房颤动的信息。然后,评价部118在由解析部114解析为有心房颤动(发病)的期间(发病期间)内,评价由测量部16测量的活动量或者由该活动量转换而来的负荷量。该评价例如根据对整个期间的负荷量的平均值和由在解析为有心房颤动的期间内测量的活动量转换而来的负荷量的平均值进行比较等,通过算出并提示被解析为有心房颤动的期间中的负荷量相比于所有的负荷量增加了百分之几等而进行。另外,该评价也可以通过将被解析为有心房颤动的期间中的负荷量的曲线图用与其以外的期间中的负荷量的曲线图的颜色不同的颜色来表示而进行。也就是说,只要通过该评价部118进行的评价,以根据心房有无颤动区别被测量的活动量或者对应于该活动量的负荷量的方式加工信息即可。
此外,评价部118既可以从行动数据存储区域126读取行动数据,也可以从自觉症状数据存储区域127读取自觉症状数据。而且,评价部118可以按照时间序列排列由测量部16测量的活动量或者对应于该活动量的负荷量,以与行动按钮141被按压的时刻或者自觉症状按钮142被按压的时刻等可对比的方式来表现。
这样,本发明的实施方式中的脉波测定装置1由于能够表示心房颤动的发病与活动量的关系,因而例如看到该显示的医生B能够将被借给脉波测定装置1的患者A的症状与活动量建立关系而进行诊断。
变形例
以上,虽然说明了本发明的实施方式,但是本发明如下所示能够以各种各样的方式来实施。
变型例1
在上述实施方式中,检测波形信号L是表示在脉波检测部20中检测脉波的结果的信号,但也可以是作为检测心电的结果而得到的波形信号。即,只要是可以取得相当于RR间隔的参数的波形信号即可。
变形例2
在上述实施方式中,虽然作为心房颤动解析装置100的功能结构而设有噪声降低部111,但也可以不一定设置。在该情况下,RR间隔算出部112可以从检测波形信号存储区域121取得进行频率解析的检测波形信号L。
变形例3
在上述实施方式中,虽然心房颤动解析装置100在脉波测定装置1中被实现,但也可以使其在个人电脑等信息处理装置中被实现。在该情况下,信息处理装置可以从外部装置取得预先测定的检测波形信号L,并将其存储于检测波形信号存储区域121。然后,信息处理装置可以通过心房颤动解析处理来解析该检测波形信号L,进而解析有无心房颤动。
变形例4
在上述实施方式中,装置主体10和脉波检测部20虽然通过电缆30有线连接,但也可以通过无线连接。在该情况下,装置主体10和脉波检测部20可以通过无线通信交换脉波检测部20的控制所需要的控制信号以及在脉波检测部20中生成的检测波形信号L等各种信号。另外,只要在装置主体10和脉波检测部20中分别具有可供电的电池等结构即可。
变形例5
在上述实施方式中,虽然心房颤动的解析结果显示于显示部15而被通知给用户,但也可以通过声音、振动等进行通知。例如,在通过声音通知用户时,只要设置扬声器以及根据来自解析部114的信息控制扬声器的放送内容的声音控制部即可。另外,在通过振动通知用户时,只要设置振动致动器以及根据来自解析部114的信息控制振动致动器的振动内容的振动控制部即可。这样,实施方式中的显示控制部115和显示部15也能够定义为根据心房颤动的解析结果通知用户的通知部。
变形例6
在实施方式中说明的各种参数例如簇群分离的阈值(±30%)、相对于平均重心坐标的规定的范围(变动系数10.0以上且功率0.5以上)、移动平均处理的数据点数(20点)以及帧期间(120秒)是一个示例,这些参数的值并不局限于此。另外,聚类的算法并不局限于k-means法。也可以通过k-means法以外的算法,使数据群分离为两个簇群。另外,有无心房颤动的具体的解析方法不限于图12中说明过的。也可以通过例如将功率和变动系数中的至少一方与阈值比较等在图12中说明的方法以外的方法来解析有无心房颤动。
变形例7
上述实施方式中的控制程序能够在存储于磁记录介质(磁带、磁盘等)、光记录介质(光盘等)、光磁记录介质、半导体存储器等计算机可读取的记录介质的状态下提供。另外,脉波测定装置1也可以经过网络下载各程序。
变形例8
在上述实施方式中,心房颤动解析系统9的加工中心C中的加工功能虽然被包含在脉波测定装置1中,但心房颤动解析系统也可以与脉波测定装置1分离而利用加工中心C。图17是表示变形例中的心房颤动解析系统9a的结构的图。心房颤动解析系统9a具有脉波测定装置1a、加工装置6以及提示装置7。加工装置6以及提示装置7都由个人电脑等构成。
医院D拥有脉波测定装置1a,在例如10天等确定的期间内将该脉波测定装置1a出借给患者A。被借给脉波测定装置1a的患者A如图2所示将该脉波测定装置1a佩戴于自身的手腕而进行脉波的测定以及心房颤动的解析。在经过了上述的期间之后,脉波测定装置1a存储的各种数据向加工中心C所有的加工装置6发送数据。数据的发送可以通过专用线路和通过公共线路进行,并且也可以通过有线和无线进行。在该情况下,脉波测定装置1a具备图3虚线所示的通信部101,经由该通信部101进行数据的发送。通信部101是在脉波测定装置1a和加工装置6之间进行信息通信的接口,例如为各种调制解调器或者无线通信线路,或者串行接口等。另外,加工中心C的职员可以将由患者A邮寄的脉波测定装置1a连接于加工装置6而使上述数据读入至该加工装置6。
加工装置6具有与上述实施方式中的脉波测定装置1的评价部118同样的功能,加工从脉波测定装置1a取得的数据。此外,在该情况下,脉波测定装置1a也可以不具有脉波测定装置1的评价部118。脉波测定装置1a将由转换部117转换的负荷量的时间序列数据和由解析部114解析的表示有无心房颤动的信息的时间序列数据向加工装置6发送。加工装置6评价接收到的负荷量中的、在由解析部114解析为有心房颤动的期间中由测量部16测量的活动量转换后的负荷量。然后,加工装置6根据负荷量的评价结果,以根据有无心房颤动区别每一时刻的负荷量(或者对应于该负荷量的活动量)的方式加工数据。
由加工装置6加工的数据被发送至提示装置7。提示装置7具有液晶等显示画面,例如通过将对应于加工装置6发送的数据的图像显示在显示画面上,从而向医生B提示该被加工的数据。
这样,由于心房颤动解析系统9a能够提示心房颤动的发病与活动量的关系,因而利用该系统的医生B能够使患者A的心房颤动的发病与活动量建立关系而进行诊断。具体而言,医生B通过利用心房颤动解析系统9a,能够诊断患者A的症状是在夜间等活动量低时心房颤容易发病的类型,还是在进行剧烈的运动等而活动量高时心房颤动容易发病的类型。
此外,加工装置6从脉波测定装置1a既可以取得行动数据,也可以取得自觉症状数据。而且,加工装置6可以按照时间序列排列由测量部16测量的活动量或者对应于该活动量的负荷量,进而以与行动按钮141被按压的时刻或者自觉症状按钮142被按压的时刻等可对比的方式来表现。图18是表示以可对比的方式表现这些数据的曲线图的一个示例的图。该图所示的曲线图通过横轴表示经过时间,心房颤动的解析结果、活动量、行动数据以及自觉症状数据按照时间序列以可对比的方式排列。被提示如此加工的数据的医生B确认例如活动量的上升期间与被解析为心房颤动发病的期间一致的程度,进而诊断患者A的症状的类型。另外,医生B通过比较行动按钮141被按压的时刻或者自觉症状按钮142被按压的时刻和被解析为心房颤动发病的时刻,能够诊断患者A在进行哪一种行动时应该注意发病,或者诊断患者A对症状的自觉的灵敏度。
变形例9
在上述实施方式或者变形例中,脉波测定装置1的评价部118或者加工装置6虽然取得由转换部117转换的负荷量的数据并评价了该负荷量,但也可以取得转换前的活动量的数据而评价活动量。在该情况下,脉波测定装置1可以不具有转换部17。
另外,在上述实施方式或者变形例中,脉波测定装置1的操作部14虽然具有在被操作时向心房颤动解析系统100输出用于使其存储“行动数据”的信号的行动按钮141,但也可以不具有行动按钮141。在该情况下,用户可以操作设于外部的输入装置向脉波测定装置1进行行动数据的输入,脉波测定装置1可以通过上述通信部101从用户操作的该输入装置接受行动数据的输入。该通信部101作为从输入装置接受行动数据的输入的接受部而起作用。此外,即使是该情况,通知部119也比较转换部117转换而得到的用户的负荷量与被存储于阈值存储区域129的该用户的阈值,在用户的负荷量超过阈值时,将该内容通知用户,以便用户对输入装置进行操作。
Claims (8)
1.一种心房颤动解析装置,其特征在于,
具有:
取得部,取得表示心电或者脉波的检测结果的时间变化的检测波形信号;
波形存储部,存储所述取得的检测波形信号;
RR间隔算出部,根据通过所述取得的检测波形信号的频率解析得到的各单位时间的频谱,对应于每一个所述单位时间算出相当于平均RR间隔的参数;
功率算出部,算出表示由所述RR间隔算出部算出的平均RR间隔的时间变化的RR波形信号的频谱中预先确定的频带的功率;
变动系数算出部,算出所述RR波形信号中的所述平均RR间隔的变动系数;
解析部,根据所述功率和所述变动系数的组合,对应于每一个时刻判定有无心房颤动;
测量部,测量用户的活动量;
第一存储部,使所述测量部测量的活动量和测量所述活动量的时刻各自相对应地进行存储;以及
评价部,在由所述解析部解析为存在的所述心房颤动的发病期间,评价由所述测量部测量的活动量。
2.根据权利要求1所述的心房颤动解析装置,其特征在于,
具有第二存储部,所述第二存储部对应于所述用户的行动的每一个时刻,存储表示所述行动的行动数据。
3.根据权利要求1或2所述的心房颤动解析装置,其特征在于,
具有:
范围存储部,存储表示所述用户的活动量能够变动的范围的范围信息;以及
转换部,根据由所述范围存储部存储的范围信息,将由所述测量部测量的活动量转换为表示该用户的负荷的负荷量。
4.根据权利要求1或2所述的心房颤动解析装置,其特征在于,
具有第三存储部,所述第三存储部对应于所述用户自觉所述心房颤动的每一个时刻,存储表示所述用户自觉的所述心房颤动的自觉症状数据。
5.根据权利要求2所述的心房颤动解析装置,其特征在于,
具有:
操作部,被所述用户操作;以及
通知部,当由所述测量部测量的活动量超过预先确定的阈值时,将所述活动量超过所述阈值的内容通知所述用户,以便所述用户对所述操作部进行操作,
所述第二存储部在所述用户操作了所述操作部时存储所述行动数据。
6.根据权利要求2所述的心房颤动解析装置,其特征在于,
具有:
接受部,接受来自输入装置的所述行动数据的输入;以及
通知部,当由所述测量部测量的活动量超过预先确定的阈值时,将所述活动量超过所述阈值的内容通知所述用户,以便所述用户从所述输入装置进行所述行动数据的输入,
所述第二存储部在所述接受部接受了所述行动数据的输入时存储所述行动数据。
7.一种心房颤动解析系统,其特征在于,
具有:
权利要求1至6中任一项所述的心房颤动解析装置;
加工装置,取得所述测量部测量的活动量,并根据由所述解析部进行解析的结果加工表示所述活动量的数据;以及
提示装置,提示由所述加工装置加工后的数据。
8.一种心房颤动解析方法,其特征在于,
具有:
取得表示心电或者脉波的检测结果的时间变化的检测波形信号的步骤;
存储所述取得的检测波形信号的步骤;
根据通过所述取得的检测波形信号的频率解析得到的各单位时间的频谱,对应于每一个所述单位时间算出相当于平均RR间隔的参数的步骤;
算出表示所述算出的平均RR间隔的时间变化的RR波形信号的频谱中预先确定的频带的功率的步骤;
算出所述RR波形信号中的所述平均RR间隔的变动系数的步骤;
根据所述功率和所述变动系数的组合,对应于每一个时刻解析有无心房颤动的步骤;
测量用户的活动量的步骤;
使所述被测量的活动量和测量所述活动量的时刻各自相对应地进行存储的步骤;以及
在由对应于每一个时刻解析有无心房颤动的步骤解析为存在的所述心房颤动的发病期间,对由测量用户的活动量的步骤所测量的活动量进行评价的步骤。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013-046648 | 2013-03-08 | ||
JP2013046648A JP2014171660A (ja) | 2013-03-08 | 2013-03-08 | 心房細動解析装置、心房細動解析システム、心房細動解析方法およびプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104027110A CN104027110A (zh) | 2014-09-10 |
CN104027110B true CN104027110B (zh) | 2018-03-30 |
Family
ID=51458273
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310739868.9A Active CN104027110B (zh) | 2013-03-08 | 2013-12-26 | 心房颤动解析装置、心房颤动解析系统及解析方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9504400B2 (zh) |
JP (1) | JP2014171660A (zh) |
CN (1) | CN104027110B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014171589A (ja) * | 2013-03-07 | 2014-09-22 | Seiko Epson Corp | 心房細動解析装置およびプログラム |
CN106037702A (zh) * | 2015-04-17 | 2016-10-26 | 精工爱普生株式会社 | 生物体信息处理系统、生物体信息处理装置及终端装置 |
JP6668606B2 (ja) * | 2015-04-17 | 2020-03-18 | セイコーエプソン株式会社 | 生体情報処理システム、生体情報処理装置及び解析結果情報の生成方法 |
US20160302736A1 (en) * | 2015-04-17 | 2016-10-20 | Seiko Epson Corporation | Biological information processing system, biological information processing device, and method for generating analysis result information |
US9910964B2 (en) * | 2015-06-25 | 2018-03-06 | Analytics For Life | Methods and systems using mathematical analysis and machine learning to diagnose disease |
JP6836264B2 (ja) * | 2016-12-09 | 2021-02-24 | デルタ工業株式会社 | 生体状態推定装置、生体状態推定方法、コンピュータプログラム及び記録媒体 |
JP7159163B2 (ja) * | 2016-12-20 | 2022-10-24 | コーニンクレッカ フィリップス エヌ ヴェ | 肺及び全身の血行動態の促進 |
CN108245148B (zh) * | 2016-12-28 | 2022-08-23 | 三星电子株式会社 | 自适应生物信号特征组合设备及方法 |
CN106901724A (zh) * | 2017-04-26 | 2017-06-30 | 河南省中医院(河南中医药大学第二附属医院) | 一种心电图测试装置及测试方法 |
JP7327816B2 (ja) * | 2018-07-10 | 2023-08-16 | 国立大学法人 香川大学 | 脈波信号の解析装置、脈波信号の解析方法およびコンピュータプログラム |
JP6964174B1 (ja) * | 2020-11-30 | 2021-11-10 | 真旭 徳山 | 情報処理装置、情報処理方法、及びプログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101534783A (zh) * | 2006-11-14 | 2009-09-16 | 皇家飞利浦电子股份有限公司 | 对运动的敏感性降低的cpr指导设备 |
JP3159276U (ja) * | 2010-02-24 | 2010-05-13 | 株式会社メディリンク | 生体情報記録紙 |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01221140A (ja) * | 1988-02-27 | 1989-09-04 | Suzuken:Kk | 長時間心電図解析および表示方法並びに装置 |
US7738936B1 (en) * | 1999-11-10 | 2010-06-15 | Pacesetter, Inc. | Methods and systems for reducing data acquisition, power and/or processing for pulse oximetry applications |
JP3846844B2 (ja) | 2000-03-14 | 2006-11-15 | 株式会社東芝 | 身体装着型生活支援装置 |
US7123953B2 (en) | 2001-12-26 | 2006-10-17 | Mediwave Star Technology Inc. | Method and system for evaluating arrhythmia risk with QT-RR interval data sets |
JP4023429B2 (ja) * | 2003-10-07 | 2007-12-19 | 株式会社デンソー | 携帯型生体情報モニタ装置 |
JP2006061416A (ja) * | 2004-08-26 | 2006-03-09 | Matsushita Electric Works Ltd | 生体電位検知システム |
KR100601981B1 (ko) * | 2005-01-14 | 2006-07-18 | 삼성전자주식회사 | 활동패턴 감시 방법 및 장치 |
US20060224072A1 (en) * | 2005-03-31 | 2006-10-05 | Cardiovu, Inc. | Disposable extended wear heart monitor patch |
JP4752673B2 (ja) | 2005-10-06 | 2011-08-17 | コニカミノルタセンシング株式会社 | 脈波データ解析方法、システム、プログラム |
JP4862730B2 (ja) | 2007-04-05 | 2012-01-25 | コニカミノルタセンシング株式会社 | 脈波データ解析方法、システム、プログラム |
US20090076559A1 (en) * | 2007-09-14 | 2009-03-19 | Corventis, Inc. | Adherent Device for Cardiac Rhythm Management |
JP2009089883A (ja) | 2007-10-09 | 2009-04-30 | Konica Minolta Sensing Inc | 心房細動検出装置、システムおよび方法 |
US9320448B2 (en) * | 2008-04-18 | 2016-04-26 | Pacesetter, Inc. | Systems and methods for improved atrial fibrillation (AF) monitoring |
JP5413887B2 (ja) * | 2009-03-31 | 2014-02-12 | フクダ電子株式会社 | 心電図解析装置および心電図解析プログラム |
WO2010132311A1 (en) * | 2009-05-13 | 2010-11-18 | Monitoring Information Technologies, Inc. | Systems and methods for heart and activity monitoring |
JP5576234B2 (ja) | 2010-10-14 | 2014-08-20 | 株式会社デンソー | 生体状態監視装置 |
US20130060154A1 (en) * | 2011-09-07 | 2013-03-07 | Seiko Epson Corporation | Atrial fibrillation decision apparatus, and method and program for deciding presence of atrial fibrillation |
JP5760876B2 (ja) | 2011-09-07 | 2015-08-12 | セイコーエプソン株式会社 | 心房細動判定装置、心房細動判定方法およびプログラム |
JP6060563B2 (ja) * | 2012-08-24 | 2017-01-18 | セイコーエプソン株式会社 | 心房細動判定装置、心房細動判定方法およびプログラム |
-
2013
- 2013-03-08 JP JP2013046648A patent/JP2014171660A/ja not_active Withdrawn
- 2013-12-26 CN CN201310739868.9A patent/CN104027110B/zh active Active
-
2014
- 2014-03-04 US US14/197,050 patent/US9504400B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101534783A (zh) * | 2006-11-14 | 2009-09-16 | 皇家飞利浦电子股份有限公司 | 对运动的敏感性降低的cpr指导设备 |
JP3159276U (ja) * | 2010-02-24 | 2010-05-13 | 株式会社メディリンク | 生体情報記録紙 |
Also Published As
Publication number | Publication date |
---|---|
US20140257123A1 (en) | 2014-09-11 |
JP2014171660A (ja) | 2014-09-22 |
CN104027110A (zh) | 2014-09-10 |
US9504400B2 (en) | 2016-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104027110B (zh) | 心房颤动解析装置、心房颤动解析系统及解析方法 | |
Schneider | Evaluations of stressful transactions: what's in an appraisal? | |
CN104602604B (zh) | 生物体信息处理系统、可穿戴装置、服务器系统及生物体信息处理系统的控制方法 | |
US20120143022A1 (en) | Physiological signal detection system capable of showing emotions, device and emotional display method | |
US9504401B2 (en) | Atrial fibrillation analyzer and program | |
US20160242672A1 (en) | Vital signal measuring apparatus and method for estimating contact condition | |
US20090149778A1 (en) | Depression detection system | |
CN106419841A (zh) | 睡眠评估方法、装置及系统 | |
KR20170087855A (ko) | 적어도 부분적으로 펄스 파형에 기반하는 자동화 진단 | |
US10524676B2 (en) | Apparatus and method for determining a health parameter of a subject | |
CN105939658A (zh) | 用于传感器的最佳定位的方法、系统和装置 | |
CN105982642A (zh) | 一种基于体震信号的睡眠检测方法及检测系统 | |
JP2009148372A (ja) | ストレス判定システムおよびストレス改善システム | |
KR20170045204A (ko) | 우울 상태 판정 방법, 및 우울 상태 판정 장치 | |
CN102462494A (zh) | 一种新型智能心电测试保健器 | |
JP2019504747A (ja) | 体状態の分類 | |
CN111588384A (zh) | 获得血糖检测结果的方法、装置及设备 | |
Mahmud et al. | SensoRing: An integrated wearable system for continuous measurement of physiological biomarkers | |
JP6060563B2 (ja) | 心房細動判定装置、心房細動判定方法およびプログラム | |
CN201894645U (zh) | 一种新型智能心电测试保健器 | |
JP4882052B2 (ja) | 自己組織化マップを用いた脈波診断システム並びに自己組織化マップの生成プログラム及び生成方法 | |
Lamarque et al. | A new concept of virtual patient for real-time ECG analyzers | |
Androutsou et al. | A smart computer mouse with biometric sensors for unobtrusive office work-related stress monitoring | |
US20210236039A1 (en) | Mobile ecg game controller apparatus | |
CN107714039A (zh) | 一种基于电子秤检测人体动脉血管硬化的方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |