CN103027669B - 用于判断脉象浮沉度的方法和设备 - Google Patents

用于判断脉象浮沉度的方法和设备 Download PDF

Info

Publication number
CN103027669B
CN103027669B CN201110317010.4A CN201110317010A CN103027669B CN 103027669 B CN103027669 B CN 103027669B CN 201110317010 A CN201110317010 A CN 201110317010A CN 103027669 B CN103027669 B CN 103027669B
Authority
CN
China
Prior art keywords
pulse
pressure
amplitude
pulse signal
opt
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
Application number
CN201110317010.4A
Other languages
English (en)
Other versions
CN103027669A (zh
Inventor
米凯云
周敏
王颖
郭松
谭思黎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to CN201110317010.4A priority Critical patent/CN103027669B/zh
Publication of CN103027669A publication Critical patent/CN103027669A/zh
Application granted granted Critical
Publication of CN103027669B publication Critical patent/CN103027669B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明提供了一种用于判断脉象浮沉度的方法,该方法包括以下步骤:在取脉位置上施加落于压力区间中的多个取脉压力,所述压力区间为[0,Pmax];在施加每个所述取脉压力时采集预定时间长度的反应所述脉象的脉搏信号,所述预定时间长度被设置为使得所采集到的脉搏信号至少包括一个脉搏周期;确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt;基于Popt在所述压力区间中的位置来判断所述脉象浮沉度;其中,Pmax为通过施压使所采集的脉搏信号的幅度大体上为零的取脉压力,并且所述多个取脉压力之间的压力差被设置为使得所述幅度最大的脉搏周期能够被采集。另外,本发明还提供了相应的用于判断脉象浮沉度的设备。

Description

用于判断脉象浮沉度的方法和设备
技术领域
本发明涉及用于判断脉象浮沉度的方法和设备。
背景技术
脉诊是根据“脉象”观察、判断病症情况的一种诊断方法,它是中医的基本技术,也是中医临床不可缺少的诊察步骤和内容。脉象是手指感觉脉搏跳动的形象,其随心脏搏动而产生,与心气的盛衰、脉道的通利和气血的盈虚直接有关。脉诊之所以重要,就是由于脉象能传递机体各部分的生理病理信息,是窥视体内功能变化的窗口,可为诊断病证提供重要依据。
脉象的种类很多,文献中常以位、数、形、势四个方面加以分析归纳。位是值脉动部位的浅深;数主要指脉动的频率和节律;形和势是指脉的波动形态和趋势状态。如在《费兆馥中医诊断学讲稿》中所描述的那样,在目前的中医理论中,可将构成各种脉象的主要因素大致归纳为脉象的部位、至数、长度、宽度、力度、流利度、紧张度、均匀度八个方面。一般而言,脉象与人体内外环境的关系非常密切,不但受年龄、性别、形体、生活起居和精神情志的影响,而且随着机体适应内外环境的自身调节,还可以出现各种生理性变异。另外,精神情志也可引起脉象的明显变化。
就脉象的部位而言,在中医脉诊理论中通常用脉象的浮沉度来表示脉动显现部位的浅深,即脉位表浅为浮脉、脉位深沉为沉脉、不浮不沉,脉位居中为平脉。在传统地由医生用手切脉进行诊断时,医生可能根据经验以轻、中、重三种不同的力度来取脉,根据脉动应指力量相应的变化来判断脉象的浮沉度。很明显,这是一种依赖于经验的主观判断方法。对于现有技术中的一些脉诊仪,往往会在其中设定若干个固定的压力范围,根据所取得的脉最明显的位置落在哪个区间来判断该脉象是浮脉、沉脉还是平脉,例如沉脉区间为200~250克。这样的确认方法以绝对压力值为主要判断方法,而被测者的个体差异(例如高矮胖瘦、年龄、性别等)可能导致这种绝对值区分法不足以准确识别脉象的浮沉度,造成很大的诊断误差。
因此,对能够精确判断脉象浮沉度的方法和设备有很大的需求。
发明内容
本发明的目标是提供用于判断脉象浮沉度的方法和设备以得到准确符合中医理论中的浮沉脉象定义的结果,从而为中医脉诊提供可靠的依据。
为了实现上述目标,本发明提供了一种用于判断脉象浮沉度的方法,所述方法包括以下步骤:在取脉位置上施加落于压力区间中的多个取脉压力,所述压力区间为[0,Pmax];在施加每个所述取脉压力时采集预定时间长度的反应所述脉象的脉搏信号,所述预定时间长度被设置为使得所采集到的脉搏信号至少包括一个脉搏周期;确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt;基于Popt在所述压力区间中的位置来判断所述脉象浮沉度;其中,Pmax为通过施压使所采集的脉搏信号的幅度大体上为零的取脉压力,并且所述多个取脉压力之间的压力差被设置为使得所述幅度最大的脉搏周期能够被采集。
根据本发明的一些实施例,所述基于Popt在所述压力区间中的位置来判断所述脉象浮沉度可以包括如果Popt<1/3Pmax,则所述脉象为浮脉,如果Popt>2/3Pmax,则所述脉象为沉脉。
根据本发明的一些实施例,所述多个取脉压力之间的压力差可以是相等的。
根据本发明的一些实施例,所述压力差可以被设置为Pmax的20~30等分。
根据本发明的一些实施例,可以通过控制用于施加取脉压力的取脉装置每次从取脉位置表面下降的距离来控制所述多个取脉压力之间的压力差。
根据本发明的一些实施例,所述取脉装置可以每次下降相等的距离。在一些实施例中,所述相等的距离可以是0.1毫米。
根据本发明的另外一些实施例,所述取脉装置可以第一次下降第一距离并且以后每次下降第二距离,其中所述第二距离小于所述第一距离。在一些实施例中,所述第二距离为所述第一距离的十分之一。
根据本发明的一些实施例,所述预定时间长度可以为5~10秒。
根据本发明的一些实施例,所述确定取脉压力Popt的步骤可以包括对所采集的脉搏信号进行平滑滤波以及去除噪声处理。
根据本发明的一些实施例,所述平滑滤波可以采用一维均值滤波方法或者高斯滤波方法。
根据本发明的一些实施例,所述一维均值滤波方法或者高斯滤波方法的滤波窗口的长度可被设置为所述脉搏信号的采样频率的0.03~0.08倍。
根据本发明的一些实施例,所述确定取脉压力Popt的步骤可以包括对所采集的脉搏信号进行去除基线漂移影响的处理。
根据本发明的一些实施例,所述去除基线漂移影响的处理可以包括确定所述脉搏信号中每个脉搏周期的起始点。
根据本发明的一些实施例,所述去除基线漂移影响的处理可以包括对所述起始点进行插值运算以得到所述基线的幅度值。
根据本发明的一些实施例,所述去除基线漂移影响的处理可以包括从所述脉搏信号中对应地减去所得到的所述基线的幅度值。
根据本发明的一些实施例,所述插值运算可以是三次样条插值运算或者分段三次Hermite插值运算。
根据本发明的一些实施例,所述确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt可以包括确定所述脉搏信号中每个脉搏周期的主峰点,并且将最大主峰点所在的脉搏周期确定为所述幅度最大的脉搏周期。
根据本发明的一些实施例,确定所述主峰点可以包括以下步骤:确定遍历窗口;从所采集的脉搏信号的起始点开始依次以每个点为所述遍历窗口的起点截取一段所述遍历窗口大小的脉搏信号;确定所截取的每段脉搏信号中幅度最大的点,并且记录该点的位置和幅度;比较相邻两个最大幅度值点的位置,如果所述两个位置之差大于所述遍历窗口的大小,则将对应的两个最大幅度值点均确定为主峰点;如果所述两个位置之差小于所述遍历窗口的大小,则去除其中幅度较小的点并且将其中幅度较大的点确定为主峰点。
根据本发明的一些实施例,所述遍历窗口的长度可以被设置为所述脉搏信号的采样频率的0.6~0.8倍。
根据本发明的一些实施例,所述确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt可以包括以下步骤:在每次采集预定时间长度的脉搏信号之后,计算所采集的该段脉搏信号的平均幅度值;并且将最大平均幅度值所在的脉搏周期确定为所述幅度最大的脉搏周期。
根据本发明的一些实施例,所述基于Popt在所述压力区间中的位置来判断所述脉象浮沉度可以包括基于所述主峰点的幅度值或者所述平均幅度值随所述取脉压力变化的曲线中上升段的角度来判断所述脉象浮沉度。
根据本发明的一些实施例,所述取脉位置可以包括寸、关、尺三个部位。
另外,本发明还提供了一种用于判断脉象浮沉度的设备,所述设备包括以下部件:用于在取脉位置上施加落于压力区间中的多个取脉压力的装置,所述压力区间为[0,Pmax];用于在施加每个所述取脉压力时采集预定时间长度的反应所述脉象的脉搏信号的装置,所述预定时间长度被设置为使得所采集到的脉搏信号至少包括一个脉搏周期;用于确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt的装置;以及用于基于Popt在所述压力区间中的位置来判断所述脉象浮沉度的装置;其中,Pmax为通过施压使所采集的脉搏信号的幅度为零的取脉压力,并且所述多个取脉压力之间的压力差被设置为使得所述幅度最大的脉搏周期能够被采集。
本发明所提供的用于判断脉象浮沉度的方法和设备充分考虑了被测者的个体化差异,通过客观有效的测量和分析过程实现了对中医理论中的脉象浮沉度的精确判断。下面将结合具体的实施例来详细描述本发明的方法和设备。
附图说明
本发明的前述和其他目标、特征和优点根据下面对本发明的实施例的更具体的说明将是显而易见的,这些实施例在附图中被示意。
图1示出了根据本发明的第一实施例的判断脉象浮沉度的方法的流程图;
图2示意性地示出了在图1所示的方法中所采集到的脉搏信号随时间变化的曲线;
图3示意性地示出了在图1所示的方法中所施加的取脉压力随时间变化的曲线;
图4示意性地示出了对图2所示的脉搏信号进行滤波以及去除基线漂移影响处理之后得到的结果;
图5示意性地示出了在图1所示的方法中所采集到的脉搏信号中每个脉搏周期的主峰点随取脉压力变化的曲线;
图6示出了根据本发明的第二实施例的确定脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt的方法的流程图。
具体实施方式
以下结合附图和具体实施方式进一步详细说明本发明。需要说明的是,附图中的各结构只是示意性说明,用以使本领域普通技术人员最佳地理解本发明的原理,其不一定按比例绘制。
图1示出了根据本发明的第一实施例的判断脉象浮沉度的方法的流程图。与常规的中医脉诊过程相似,为了获得与患者的脉搏跳动状态有关的信息,必须对患者身体的适当部位,即在取脉位置上施加取脉压力P。在本发明的实现过程中,使用恰当的取脉装置来代替医生实施取脉过程,该取脉装置可以是任何能够在取脉位置上施加一定压力范围内的多个取脉压力并且感测相应的人体脉搏波动信号的装置。由此,可以精确地控制每次所施加的取脉压力的大小并且记录反馈的脉搏信号。就取脉位置而言,在以下的说明中将以中医理论中所定义的寸口脉为例来描述本发明的方法的应用,然而应理解的是本发明能够适用于中医理论中所涉及的所有诊脉部位。所述寸口脉就是指桡骨茎突内侧的一段桡动脉,其分为寸、关、尺三个部分,可以将本发明的方法分别用于这三个部位,以了解各部脉象的浮沉特征。
如图1的流程图所示,在步骤S101中从0开始逐步提高所施加的取脉压力,以获得多个不同取脉压力下的脉搏信号响应。ΔP表示每次所施加的取脉压力之间的压力差,在实践中,这个压力差可以每次都是相同的值,也可以是不同的值。然而,根据本发明所提供的方法,不管每次如何对压力差ΔP取值,所有施加的取脉压力P都将落在压力范围[0,Pmax]内。所述压力范围的下限表示尚未在所述取脉位置上施加任何压力时的状态,此时的脉搏信号输出相应地也必然为0。根据常识可知,随着在寸口脉位置上所施加的取脉压力逐渐变大,由于桡动脉的血流逐渐被抑制而输出的脉搏信号将逐渐变弱,直到重新变为0,Pmax即表示使得所感测的脉搏信号大体上刚好重新回到0位所需要施加的压力值。由此可见,该压力范围[0,Pmax]可以间接地表示脉搏信号的完整范围。另外,根据本发明所提供的方法,不管压力差ΔP如何取值,都应保证脉搏信号响应中幅度最大的脉搏周期能够被采集。
可以通过多种方法来确定步骤S101中的ΔP。在一些实施例中,可以采用等压力差加压的方式,并且可以根据Pmax来确定ΔP。具体而言,可以将取脉装置放置在取脉位置上从刚接触皮肤开始逐步加压直到感测的脉搏信号输出为0,记录在这点上所施加的压力作为Pmax。接着,将0~Pmax划分为N个等分,即ΔP=Pmax/N。N的范围可以在20~30之间,这可以例如根据患者的年龄、性别、身体状况等条件来设置。应注意的是,并不一定需要在最开始的时候专门测定Pmax,而是可以在随后逐步加压的过程中动态地判断是否达到该Pmax
在另外一些实施例中,可以通过控制取脉装置每次从取脉位置表面下降的距离来控制ΔP以产生多个不同的取脉压力。举例来说,可以使得取脉装置从皮肤表面开始每次均匀下降0.1毫米,并且记录相应的压力值。由于人体皮肤的压力并非线性的,因此采用等距离下降的方法可能有时不一定能够保证足够的压力分辨率以获得幅度最大的脉搏周期,或者可能由于选取的距离间隔过小而造成过量的处理负担。在一些实施例中,可以使取脉装置在第一次加压时下降较大的第一距离并且以后每次下降较小的第二距离,该第二距离例如可以是第一距离的十分之一,从而既可以保证足够的压力分辨率,又可以避免不必要的数据冗余。
在步骤S102中,保持施加步骤S101中的取脉压力P,同时以一定的采样频率连续采集预定时间长度的脉搏信号波形。在实践中,采集时间的设置需要满足的条件是保证所获得的脉搏信号波形中至少包括一个完整的脉搏周期。优选地,进行采集的时间长度范围可以在5秒~10秒之间,采样频率可以例如为100Hz。在一些实施例中,可以在每次采集完成之后直接对所采集到的信号进行实时处理,这将在下文中结合图6进行说明。在图1所示的方法中,每次采集完成之后将对所获得的脉搏信号数据进行保存。
在步骤S103中,判断当前所采集的脉搏信号波形是否大体上为0。如上文所述,在施压的情况下脉搏信号呈现为0是由于桡动脉的血流在大的压力下被抑制而使输出的脉搏信号变弱。因此,在步骤S103中进行判断的目的在于监控是否已经达到压力范围[0,Pmax]的上限,也即判断是否已经采集到可供进行有效分析的完整范围的脉搏信号。如果否,则重复上述步骤S101及步骤S102。如果是,则将进入步骤S104以开始对所采集到的脉搏信号数据做进一步处理。
图2示意性地示出了步骤S101至步骤S103的循环结束之后所采集到的脉搏信号随时间变化的曲线。该曲线的横坐标为示出为采样点数,其可以被看作是采样的时间长度,而纵坐标表示所感测到的脉搏信号的幅度。从图2中可以明显地看到随着取脉压力的不断升高,脉搏信号从弱到强,进而又由强转弱的趋势。图中所示的每个脉冲尖峰可以被认为大致表示一个脉搏周期。
图3示意性地示出了在图1所示的方法中所施加的取脉压力随时间变化的曲线。从图3所示的情况中,取脉压力基本上是呈线性增长的,也就是说在步骤S101中每次提高的ΔP为同一固定的值。
继续回到图1所示的流程图,当已经在步骤S103中确定采集完成之后,该过程前进至步骤S104。在步骤S104中,首先对脉搏信号进行平滑滤波并且去除噪声,这可以使用本领域中已知的各种技术来实现。举例来说,可以采用一维均值滤波或者高斯滤波。滤波器窗口的大小可以根据脉搏周期来进行选择,而脉搏信号中所呈现的脉搏周期又取决于所使用的采样频率的大小,因而滤波器窗口大小的范围可以例如是采样频率的0.03~0.08倍之间,这样既能去除噪声,同时又能保证脉搏信号的波形不失真。
接着,在步骤S105中,对经滤波和去噪的脉搏信号进行去除基线漂移影响的处理。在目前的脉搏信号检测过程中,通常存在由于呼吸运动和身体移位所引起的脉搏信号的基线漂移,这种基线漂移现象的存在可能会严重影响脉搏信号的可靠性。在本发明的一些实施例中,首先确定脉搏信号中每个脉搏周期的起始点,进而根据这些起始点的数值进行插值运算以拟合出基线。所述插值运算可以例如采用三次样条插值或者分段三次Hermite插值,但应理解的是本发明可以采用任何能够恰当地拟合出基线的插值算法。最后,从脉搏信号在各个采样点上的采样值中减去对应点上的基线值,从而得到去除基线漂移影响的脉搏信号波形,保证了脉搏信号波形的一致性。图4示意性地示出了对图2所示的脉搏信号进行滤波以及去除基线漂移影响处理之后得到的结果。
在步骤S106中,根据经滤波去噪和去基线处理的脉搏信号来确定每个脉搏周期中的主峰点的位置和幅度。在本发明的一些实施例中,利用窗口遍历的方法来确定所有有效的主峰点。具体而言,从脉搏信号的起始点开始每次截取遍历窗口长度的一段脉搏信号,依次以每个采样点作为该遍历窗口的起点。接着,找到每段脉搏信号中幅度最大的采样点,记录该点的位置和幅度。依次比较相邻两个最大幅度值点的位置,如果这两个位置之差大于所述遍历窗口的长度,则将这两个最大幅度值点都确定为主峰点。相反地,如果这两个位置之差小于所述遍历窗口的长度,则去除其中幅度较小的点,而将幅度较大的点确定为主峰点。优选地,该遍历窗口的长度范围可以在采样频率的0.6~0.8倍之间,这样既能保证检测到每个脉搏周期的主峰点,又能保证在每个脉搏周期中只检测到一个主峰点。例如,在采样频率为100Hz的情况下,该遍历窗口的长度可以大致为采样频率的0.6倍,即大约60个采样点。
在步骤S107中,确定幅度值最大的主峰点所对应的取脉压力Popt。事实上,在这个步骤中可以同时确定每个主峰点分别所对应的取脉压力并且画出每个脉搏周期的主峰点随取脉压力变化的曲线,如图5所示。在图5中,横坐标表示取脉压力的大小,而纵坐标表示主峰点的幅度值,也就是每个脉搏周期中最大的脉搏波动幅度。图中用红色圆点标注的点为各个主峰点中幅度值最大的主峰点,其所对应的取脉压力即为Popt
在步骤S108中,根据Popt在压力范围[0,Pmax]中的相对位置来判断脉象浮沉度。如在上文中所描述的那样,压力范围[0,Pmax]可以表示脉搏信号的完整范围,也就是能够在特定的取脉位置上感测到人体脉搏波动的整个范围。相应地,Popt代表了可以感觉最强的脉搏波动的取脉压力,因而Popt在压力范围[0,Pmax]中所处的相对位置代表了最强的脉搏波动在整个脉搏波动范围中的相对位置。由此可见,本发明的方法所提供的判断标准完全符合中医理论中对脉象浮沉度的定义。以图5为例,该Popt恰好处在整个压力范围大致居中的位置上,因此可以认为所感测到的脉象为平脉。优选地,可以认为当Popt<1/3Pmax时,所感测到的脉象为浮脉,而Popt>2/3Pmax时所感测到的脉象为沉脉,而当Popt落在这两个界限之间时则可认为是平脉。
在本发明的另外一些实施例中,还可以基于所述主峰点的幅度值随所述取脉压力变化的曲线中上升段的角度来判断所述脉象浮沉度。参照图5,所述上升段的角度例如可以用最大主峰点和坐标原点的连线与横轴的夹角来表示。对于浮脉的情况,该角度一般偏大,例如大于45度,而对于沉脉的情况,该角度一般偏小,例如小于45度。
从上述分析中可以看出,压力范围[0,Pmax]对于每个被测对象都可能是不同的,本发明所提供的判断脉象浮沉度的方法充分考虑了个体差异。同样地,利用取脉装置施加多个不同压力并且通过分析比较得到与最大脉搏波动所对应的取脉压力也是针对不同被测对象的客观精确的过程。由此,通过使用本发明所提供的方法能够可靠高效地得到准确符合中医理论中的脉象浮沉度定义的结果,大大有利于中医的发展和推广。
图6示出了根据本发明的第二实施例的确定脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt的方法的流程图。在相对于图1所描述的实施例中,取脉压力Popt及其相对位置的确定是在所有取脉压力下的数据采集都完成之后进行的。在实践中,还可以通过简化的方法实时地确定取脉压力Popt及其相对位置,如图6的流程图所示。
在图6的步骤S201中,在取脉位置上施加初始取脉压力P。与图1所描述的方法相同,该初始取脉压力P以及随后逐渐提高的各个取脉压力也都落在压力范围[0,Pmax]中。
在步骤S202中,保持施加步骤S201中的取脉压力P,同时以一定的采样频率连续采集预定时间长度的脉搏信号波形,这与图1中的步骤S102相似,
不同于在图1所示的方法中将每次采集之后的原始脉搏信号数据全部进行保存,在步骤S203中对所采集的该段脉搏信号中每个采样点的幅度值取平均,得到平均幅度值Hmax并且仅对该值进行保存。
接着,在步骤S204中将初始取脉压力P的压力值提高ΔP,这与图1中的步骤S101相似,其中初始采样压力P以及ΔP的确定均与根据图1所描述的方法相同。
在步骤S205中,在提高的取脉压力下采集大约10秒的脉搏信号波形,并且在步骤S206中同样对该次采集的脉搏信号中所有采样点的幅度值取平均,得到平均幅度值Hnext并且对该值进行保存。
在步骤S207中,将Hnext与所保存的Hmax进行比较。如果Hnext>Hmax,则将Hnext作为Hmax重新保存,如步骤S208所示。如果Hnext<=Hmax,则保持Hmax不变,如步骤S209所示。在此之后重复步骤S204至步骤S209,直到所施加的取脉压力使得所采集的脉搏信号输出大体上为0。在整个逐步加压采集数据的过程结束之后,可以容易地将最后所保存的Hmax所对应的取脉压力确定为图1所示的方法中所定义的Popt,即可以感觉最强的脉搏波动的取脉压力。
同样地,可以画出在整个过程中所获得的每个平均幅度值随取脉压力变化的曲线,这与图5所示的曲线相似。由此,可以采用与根据图1所描述的相同的方法进行后续的脉象浮沉度判断。
图6示意了一种实时操作的简化方法来确定脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt,使用这种方法也同样能够得到准确符合中医理论中的脉象浮沉度定义的结果
应当说明的是,以上具体实施方式仅用以说明本发明的技术方案而非对其进行限制。尽管参照上述具体实施方式对本发明进行了详细的说明,本领域的普通技术人员应当理解,依然可以对本发明的具体实施方式进行修改或对部分技术特征进行等同替换而不脱离本发明的实质,其均涵盖在本发明请求保护的范围中。

Claims (23)

1.一种用于判断脉象浮沉度的方法,其特征在于,所述方法包括以下步骤:
在取脉位置上施加落于压力区间中的多个取脉压力,所述压力区间为[0,Pmax];
在施加每个所述取脉压力时采集预定时间长度的反应所述脉象的脉搏信号,所述预定时间长度被设置为使得所采集到的脉搏信号至少包括一个脉搏周期;
确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt
基于Popt在所述压力区间中的位置来判断所述脉象浮沉度;
其中,Pmax为通过施压使所采集的脉搏信号的幅度大体上为零的取脉压力,并且所述多个取脉压力之间的压力差被设置为使得所述幅度最大的脉搏周期能够被采集,
其中所述确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt包括确定所述脉搏信号中每个脉搏周期的主峰点,并且将最大主峰点所在的脉搏周期确定为所述幅度最大的脉搏周期,
其中确定所述主峰点包括以下步骤:
确定遍历窗口;
从所采集的脉搏信号的起始点开始依次以每个点为所述遍历窗口的起点截取一段所述遍历窗口长度的脉搏信号;
确定所截取的每段脉搏信号中幅度最大的点,并且记录该点的位置和幅度;
比较相邻两个最大幅度值点的位置,如果所述两个位置之差大于所述遍历窗口的长度,则将对应的两个最大幅度值点均确定为主峰点;如果所述两个位置之差小于所述遍历窗口的长度,则去除其中幅度较小的点并且将其中幅度较大的点确定为主峰点。
2.如权利要求1所述的方法,其特征在于,其中所述基于Popt在所述压力区间中的位置来判断所述脉象浮沉度包括如果Popt<1/3Pmax,则所述脉象为浮脉,如果Popt>2/3Pmax,则所述脉象为沉脉。
3.如权利要求1所述的方法,其特征在于,其中所述多个取脉压力之间的压力差相等。
4.如权利要求3所述的方法,其特征在于,其中所述压力差被设置为Pmax的20~30等分。
5.如权利要求1所述的方法,其特征在于,其中通过控制用于施加取脉压力的取脉装置每次从取脉位置表面下降的距离来控制所述多个取脉压力之间的压力差。
6.如权利要求5所述的方法,其特征在于,其中所述取脉装置每次下降相等的距离。
7.如权利要求6所述的方法,其特征在于,其中所述相等的距离为0.1毫米。
8.如权利要求5所述的方法,其特征在于,其中所述取脉装置第一次下降第一距离并且以后每次下降第二距离,所述第二距离小于所述第一距离。
9.如权利要求8所述的方法,其特征在于,其中所述第二距离为所述第一距离的十分之一。
10.如权利要求1所述的方法,其特征在于,其中所述预定时间长度为5~10秒。
11.如权利要求1所述的方法,其特征在于,其中所述确定取脉压力Popt的步骤包括对所采集的脉搏信号进行平滑滤波以及去除噪声处理。
12.如权利要求11所述的方法,其特征在于,其中所述平滑滤波采用一维均值滤波方法或者高斯滤波方法。
13.如权利要求12所述的方法,其特征在于,其中所述一维均值滤波方法或者高斯滤波方法的滤波窗口的长度被设置为所述脉搏信号的采样频率的0.03~0.08倍。
14.如权利要求1所述的方法,其特征在于,其中所述确定取脉压力Popt的步骤包括对所采集的脉搏信号进行去除基线漂移影响的处理。
15.如权利要求14所述的方法,其特征在于,其中所述去除基线漂移影响的处理包括确定所述脉搏信号中每个脉搏周期的起始点。
16.如权利要求15所述的方法,其特征在于,其中所述去除基线漂移影响的处理包括对所述起始点进行插值运算以得到所述基线的幅度值。
17.如权利要求16所述的方法,其特征在于,其中所述去除基线漂移影响的处理包括从所述脉搏信号中对应地减去所得到的所述基线的幅度值。
18.如权利要求16或17所述的方法,其特征在于,其中所述插值运算为三次样条插值运算或者分段三次Hermite插值运算。
19.如权利要求1所述的方法,其特征在于,其中所述遍历窗口的长度被设置为所述脉搏信号的采样频率的0.6~0.8倍。
20.如权利要求1所述的方法,其特征在于,其中所述基于Popt在所述压力区间中的位置来判断所述脉象浮沉度包括基于所述主峰点的幅度值随所述取脉压力变化的曲线中上升段的角度来判断所述脉象浮沉度。
21.如权利要求1所述的方法,其特征在于,其中所述取脉位置包括寸、关、尺三个部位。
22.一种用于判断脉象浮沉度的设备,其特征在于,所述设备包括以下部件:
用于在取脉位置上施加落于压力区间中的多个取脉压力的装置,所述压力区间为[0,Pmax];
用于在施加每个所述取脉压力时采集预定时间长度的反应所述脉象的脉搏信号的装置,所述预定时间长度被设置为使得所采集到的脉搏信号至少包括一个脉搏周期;
用于确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt的装置;以及
用于基于Popt在所述压力区间中的位置来判断所述脉象浮沉度的装置;
其中,Pmax为通过施压使所采集的脉搏信号的幅度为零的取脉压力,并且所述多个取脉压力之间的压力差被设置为使得所述幅度最大的脉搏周期能够被采集,
其中所述用于确定所采集到的脉搏信号中幅度最大的脉搏周期所对应的取脉压力Popt的装置包括用于确定所述脉搏信号中每个脉搏周期的主峰点的装置,以及用于将最大主峰点所在的脉搏周期确定为所述幅度最大的脉搏周期的装置,
其中用于确定所述主峰点的装置包括:
用于确定遍历窗口的装置;
用于从所采集的脉搏信号的起始点开始依次以每个点为所述遍历窗口的起点截取一段所述遍历窗口长度的脉搏信号的装置;
用于确定所截取的每段脉搏信号中幅度最大的点,并且记录该点的位置和幅度的装置;
用于比较相邻两个最大幅度值点的位置的装置,如果所述两个位置之差大于所述遍历窗口的长度,则将对应的两个最大幅度值点均确定为主峰点;如果所述两个位置之差小于所述遍历窗口的长度,则去除其中幅度较小的点并且将其中幅度较大的点确定为主峰点。
23.如权利要求22所述的设备,其特征在于,其中所述基于Popt在所述压力区间中的位置来判断所述脉象浮沉度包括如果Popt<1/2Pmax,则所述脉象为浮脉,如果Popt>2/3Pmax,则所述脉象为沉脉。
CN201110317010.4A 2011-09-30 2011-09-30 用于判断脉象浮沉度的方法和设备 Active CN103027669B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110317010.4A CN103027669B (zh) 2011-09-30 2011-09-30 用于判断脉象浮沉度的方法和设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110317010.4A CN103027669B (zh) 2011-09-30 2011-09-30 用于判断脉象浮沉度的方法和设备

Publications (2)

Publication Number Publication Date
CN103027669A CN103027669A (zh) 2013-04-10
CN103027669B true CN103027669B (zh) 2016-09-14

Family

ID=48015458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110317010.4A Active CN103027669B (zh) 2011-09-30 2011-09-30 用于判断脉象浮沉度的方法和设备

Country Status (1)

Country Link
CN (1) CN103027669B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103932687B (zh) * 2014-04-22 2017-04-12 北京印刷学院 一种脉象信号预处理方法和装置
CN107595249B (zh) * 2017-09-30 2021-03-23 深圳前海全民健康科技有限公司 基于脉搏波的怀孕女性筛查方法
TWI777099B (zh) * 2019-10-30 2022-09-11 友達光電股份有限公司 檢測裝置
CN113288075B (zh) * 2021-06-01 2022-07-19 山西和脉互联网智慧中医产业技术研究院有限公司 浮中沉脉位识别方法及系统
CN114431840A (zh) * 2022-04-08 2022-05-06 北京大学深圳研究生院 一种脉搏采集装置、脉搏采集方法和系统
CN114711733B (zh) * 2022-06-07 2023-07-28 北京大学深圳研究生院 一种脉搏信号提取方法、装置、电子设备及存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86208967U (zh) * 1986-11-12 1987-09-02 上海医疗器械研究所 中医脉象检测装置
CN1543911A (zh) * 2003-11-13 2004-11-10 李士春 中医脉诊传感系统
CN1759804A (zh) * 2005-11-02 2006-04-19 浙江大学 中医四诊合参智能辨证方法
CN1792319A (zh) * 2005-11-03 2006-06-28 浙江大学 中医脉象特征参数自动检测方法
CN1813628A (zh) * 2006-01-18 2006-08-09 哈尔滨工业大学 自动调节切脉压力的脉象检测装置及脉象特征提取方法
CN101193587A (zh) * 2006-06-21 2008-06-04 大尧Medi有限公司 一种使用压力传感器阵列的脉搏分析仪
CN201481407U (zh) * 2009-09-21 2010-05-26 上海道生医疗科技有限公司 脉象采集装置
CN101732033A (zh) * 2008-11-07 2010-06-16 中国科学院计算技术研究所 一种人体波形中的特征参数的提取方法和装置
CN101978930A (zh) * 2010-12-01 2011-02-23 岑澔 一种因人制宜的、点阵式、自适应的脉象检测装置及其方法
CN102144916A (zh) * 2011-04-21 2011-08-10 华东理工大学 一种可自动调压的多通道脉搏信号检测方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100869514B1 (ko) * 2007-03-14 2008-11-19 한국 한의학 연구원 허맥 및 실맥 진단 시스템

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86208967U (zh) * 1986-11-12 1987-09-02 上海医疗器械研究所 中医脉象检测装置
CN1543911A (zh) * 2003-11-13 2004-11-10 李士春 中医脉诊传感系统
CN1759804A (zh) * 2005-11-02 2006-04-19 浙江大学 中医四诊合参智能辨证方法
CN1792319A (zh) * 2005-11-03 2006-06-28 浙江大学 中医脉象特征参数自动检测方法
CN1813628A (zh) * 2006-01-18 2006-08-09 哈尔滨工业大学 自动调节切脉压力的脉象检测装置及脉象特征提取方法
CN101193587A (zh) * 2006-06-21 2008-06-04 大尧Medi有限公司 一种使用压力传感器阵列的脉搏分析仪
CN101732033A (zh) * 2008-11-07 2010-06-16 中国科学院计算技术研究所 一种人体波形中的特征参数的提取方法和装置
CN201481407U (zh) * 2009-09-21 2010-05-26 上海道生医疗科技有限公司 脉象采集装置
CN101978930A (zh) * 2010-12-01 2011-02-23 岑澔 一种因人制宜的、点阵式、自适应的脉象检测装置及其方法
CN102144916A (zh) * 2011-04-21 2011-08-10 华东理工大学 一种可自动调压的多通道脉搏信号检测方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于虚拟仪器的脉搏信号处理方法研究及应用;于璐;《中国医科大学硕士学位论文》;20090928;正文第4.1、4.4节 *
脉象检测中最大取脉压力和人体耐受安全值的研究;王晶晶等;《辽宁中医杂志》;20100411;第37卷(第4期);全文 *

Also Published As

Publication number Publication date
CN103027669A (zh) 2013-04-10

Similar Documents

Publication Publication Date Title
CN103027669B (zh) 用于判断脉象浮沉度的方法和设备
EP3358485A1 (en) General noninvasive blood glucose prediction method based on timing analysis
US20180350468A1 (en) Processing physiological electrical data for analyte assessments
Wu et al. Evaluating physiological dynamics via synchrosqueezing: Prediction of ventilator weaning
CN105266772A (zh) 一种生理参数的测量方法
CN105997043B (zh) 一种基于腕式可穿戴设备的脉率提取方法
CN105975780A (zh) 一种基于虚拟现实交互技术的机器问诊系统
WO2016104538A1 (ja) 呼吸状態推定装置、携帯機器、装着型デバイス、プログラム、媒体、呼吸状態推定方法及び呼吸状態推定器
CN103263271A (zh) 非接触式自动血氧饱和度测量系统及测量方法
CN103405225B (zh) 一种获取疼感评测指标的方法、装置及设备
CN108289627A (zh) 对呼吸性窦性心律不齐进行量化的方法和装置以及这种方法或这种装置的应用
CN105662598B (zh) 一种大脑皮层功能区定位装置、方法和系统
Kim et al. Novel diagnostic model for the deficient and excess pulse qualities
CN106999071A (zh) 脉波分析装置
Schrumpf et al. Derivation of the respiratory rate from directly and indirectly measured respiratory signals using autocorrelation
CN105852863B (zh) 一种呼吸率测量方法及装置
CA3113800A1 (en) Model-based sensor technology for detection of cardiovascular status
Klum et al. Minimally spaced electrode positions for multi-functional chest sensors: ECG and respiratory signal estimation
CN104305958B (zh) 一种极短时自主神经状态的光电容积波多参量分析方法
US8027717B2 (en) Method and apparatus for monitoring the sedation level of a sedated patient
Hu et al. Pulse differences and 3D pulse mapping in TPNI displacements
CN111657873B (zh) 一种基于可见光和近红外光谱技术的体质预测方法
Wolling et al. Optimal preprocessing of raw signals from reflective mode photoplethysmography in wearable devices
More et al. Development of non-invasive diagnostic tool for diseases using Photo Plethysmography
JP6482412B2 (ja) 粘弾特性取得装置、粘弾特性取得方法、粘弾特性取得プログラム、及びそのプログラムを記録する記録媒体

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant