一种心电数据快速重采样方法
【专利摘要】心电数据采样频率转换的传统方法是，先对原始信号进行整数倍插值，达到原始频率与目标采样频率的最小公倍数，再对插值后信号进行整数倍抽取，从而得到目标采样频率要求的数据。但是传统方法中存在的严重问题是，当插值倍数较大时，运行速度很慢，占用存储空间过大，达不到实际应用的要求。本发明为一种心电数据快速重采样方法，解决了这个问题。快速重采样方法的特征在于，在达到工程要求以及符合采样定律的前提下，实现心电数据的任意采样频率转换。这种方法可以很好地应用于实际工程使用，而且对于其它要满足工程需求的信号重采样也非常有意义。
【专利说明】一种心电数据快速重采样方法
【技术领域】
[0001] 心电数据快速重采样方法属于采样控制系统领域，主要适用于心电数据的重采 样，得到不同采样频率的数据，用于进一步研究；对于其它有同样需求的离散数据重采样同 样适用。
【背景技术】
[0002] 改革开放以来，我国医院信息化建设得到了长足发展，特别是医院信息系统 (Hospital Information System,HIS)对于提高医院科学管理水平，以及推动医院现代化建 设发挥了重要作用。而实现信息标准化以完成信息的统一管理，是目前医院HIS建设中面 临的首要问题。
[0003] 当前实现远程传输的数据主要有心电、体温、血压、血氧饱和度等，这些数据信息 需要实现标准管理。在上述数据的格式统一过程中，最复杂的当属心电数据，心电数据检测 设备多样化，不同设备的数据格式并不相同，而且心电数据量大，对于诊断病人病情至关重 要，所以心电数据格式的统一势不容缓。
[0004] 心电数据格式转换中的信息包括电设备的导联数、病人信息、增益、采样频率等， 而由于心电分析软件大都采用差值运算或变换的方法，来对心电原始数据进行预处理及分 析，因而增益等其他参数对于心电诊断分析的影响不大；但是采样频率如果与目标分析软 件不一致，那么心电波形的显示将会出现错误，甚至造成无法诊断或诊断错误。
[0005] 另外，当前研究心电方法等都是基于国际标准心电数据库进行研究的，当前全球 公认的标准心电数据库有三个，分别是美国麻省理工学院提供的MIT-BIH数据库，美国心 脏学会的AHA数据库以及欧洲AT-T心电数据库。其中，MIT-BIH数据库的应用相对广泛， MIT-BIH数据库的采样频率固定为360Hz，而国内很多心电检测设备的采样频率却有很多 变化，将标准数据库的采样频率进行转换，可以采用其来完成特定频率下的分析，从而可以 实现与标准的对比，这对于校准非常重要。
[0006] 传统实现采样频率转换的方法即严格按照离散时间信号的内插与抽取理论，先对 待转换信号进行整数倍的插值，达到原始采样频率与目标采样频率值的最小公倍数，然后 再对插值后的数据进行整数倍抽取。例如若将心电数据从MIT-BIH数据库的固定采样频率 360Hz变为250Hz，传统方法需要先将数据内插到9000Hz，即内插了 25倍，然后再进行抽取 36倍，这样虽然保证了重采样的数据完整，但是运行效率很低，占用存储空间过大。
【发明内容】
[0007] 本发明的目的在于提供一种心电数据的快速重采样方法，对于心电格式转换及标 准化，以及心电检测设备的校验都有非常重要的作用，同样也适用于有类似需求离散数据 重采样的实际应用。
[0008] 本发明使用的相关理论包括：采样定律，离散时间信号的内插与抽取（整数倍）， Lagrange插值公式。
[0009] 我们所涉及的心电数据快速重采样与上述技术不同之处，在于将内插的倍数改为 分数倍内插，然后再使用Lagrange插值公式进行重新均匀分布。新方法流程图如图1所 示，可以实现不同采样频率间的转换，根据原始采样频率值与目标采样频率值的关系，分两 种情况进行。
[0010] (1)若原始采样频率值小于目标采样频率值，步骤如下：
[0011] 步骤1，先对数据进行整数倍内插，达到小于目标采样频率的原始值最大整数倍 (此值为1时可省略），设该值为N(N > 1)，那么在原始数据相邻两个数据之间插入（N-1) 个数据，数据值为两个数据中前者的值；
[0012] 步骤2,在完成整数倍内插后，就可以对其进行分数倍内插，得到采样频率为目标 频率值的数据；
[0013] 步骤3,最后对分数倍内插后的结果通过Lagrange插值公式进行重新均匀分布。
[0014] (2)若原始采样频率值大于目标采样频率值，步骤如下：
[0015] 步骤1，先对数据进行分数倍内插，达到大于原始采样频率值的目标值的最小整数 倍，设该值为M，Μ必定大于1 ;
[0016] 步骤2,对分数倍内插的结果通过Lagrange插值公式进行重新均勻分布；
[0017] 步骤3,最后对结果进行整数倍抽取，抽取倍数为Μ倍，即从相邻Μ个数据中，抽取 其中的第一个值作为结果，便得到采样频率为目标值的数据。
[0018] 本发明中先对数据进行分数倍内插，即将原先每相邻两个数据之间插入数据改 为，每相邻多个数据插入一个数据；对于分数倍内插造成的不均匀分布，采用对数据进行重 新插值计算的方式消除某些位置过采样的问题。插值并进行均匀分布的示意图如图2所 示，而普通插值未均匀分布的示意图如图3所示，均匀分布的效果图，如图4所示。
[0019] 当然也要考虑到特殊情况：原始采样频率为目标采样频率的整数倍，以及目标采 样频率正好为原始频率值的整数倍。前者在进行插值到大于原始频率的目标频率值的最小 整数倍之前，经过判定之后，原始频率正好为目标频率的整数倍，则无需进行插值，可以直 接进行整数倍抽取，达到目标频率。而后者，在插值到小于等于目标采样频率的原始频率的 最大整数倍之后，如果目标采样频率为原始频率的整数倍，那么经过插值后，达到的频率即 为目标频率，无需继续分数倍插值及均匀分布。
【专利附图】
【附图说明】
[0020] 图1 :本发明提出的方法流程图；
[0021] 图2 :插值后均匀分布示意图；
[0022] 图3 :未执行均匀分布效果图；
[0023] 图4:均匀分布效果图。
[0024] 具体实现方式
[0025] 本发明中，可以实现任意频率间的转换，但是如果要能正常显示，必须要符合采样 定理，即采样频率必须大于被采样数据最大频率的2倍。
[0026] 以原始采样频率为360Hz的ΜΙΤ-ΒΙΗ数据库中的心电数据转换为其他采样频率来 说明，本发明的采样频率转换分为两种情况，设置目标采样频率值分别为250Hz和800Hz， 从两种情况更加清晰的说明方法的实现方式，在这里进行说明的数据都以一维数组的形式 存放。特殊情况不再进行详细说明，如将360Hz减采样转换为180Hz，或者增采样为720Hz， 前者只需要进行2倍抽取，后者只需要进行2倍内插，即可以达到目的频率值，且都可以归 到一般情况进行讨论，不再赘述。
[0027] 1)原始采样频率大于目标采样频率：360Hz到250Hz
[0028] 为了简化问题，取Is的数据大小为例，且心电数据为单导联，数据为8bit的无符 号整型数据。设原始数据数组为src[360]，分数倍内插的结果数组为mid[500]，最终结果 的数组为result [250];
[0029] (1)对数据进行分数倍插值：是之达到大于原始值的目标采样频率值的最小整 数倍,此例中该值为500Hz，即先将src进行分数倍插值，达到采样频率为500Hz的数据 mid[500]。
[0030] 插值比率rate = (500-360)/360 = 7/18 ;使用原始数据数组的下标值与此插值 比率相乘，若所得结果的整数部分比前一项的值增长1，那么就在该位置插入一个新数据， 数据值为原始数据数组在前一位置的值。例如从数组下标位置〇到2,与插值比率7/18相 乘的结果，其整数部分一直为〇,而下标3与插值比率7/18相乘，其整数部分为1，则在新数 组中下标3的位置插入一个新数据，其值为原始数据数组下标为2处的值。
[0031] 以变量posl表示src数组的下标位置，count表示当前已经插入新数据的个数， 初始值为〇 ;
[0032]
【权利要求】
1. 心电数据快速重采样方法，其特征在于：使用分数倍内插的内插方法，从而减少运 行时间及占用存储空间，使得在达到工程应用要求以及符合采样定律的前提下，可快速实 现心电数据的采样频率任意转换。方法主要包括以下几个部分； 1) 分数倍内插：通过对原始数据进行分数倍内插，达到目标采样频率或者其大于原始 采样频率的最小整数倍； 2) 均匀分布：由于采用分数倍的插值方法，造成不均匀采样的问题，所以要在插值之 后，在原始数据及插入的数据基础上，对其进行重新均匀分布； 3) 整数倍抽取与内插：若原始采样频率大于目标采样频率，在完成内插之后，要对数 据进行整数倍抽取,达到目标采样频率；若原始采样频率小于目标采样频率，那么需要先 对数据进行整数倍的内插，使之达到小于目标采样频率的最大整数倍值（此值为1时可省 略），后续采用分数倍内插达到目标采样频率值。
2. 根据权利1所述的心电数据快速重采样方法，其特征在于：采用分数倍内插的方法， 直接达到目标频率或其大于原始频率的最小整数倍，而不是传统方法中整数倍内插到原始 频率与目标采样频率的最小公倍数。主要包括以下内容： 1) 插值比率：新方法中的插值比率为分数，且小于1，其意义为每几个数值需要插入一 个新值； 2) 内插过程：原始数据读取到一个数组中，使用计算得到分数倍的插值比率，在插值 过程中，将原始数据数组下标与插值比率相乘，获取其整数部分的值，当整数部分增长1 时，即在当前的下标处，插入一个数据，其数值等于数组在当前位置的值。
3. 根据权利2所述的心电数据快速重采样方法，其特征在于：对原始数据进行分数倍 内插后，会造成某些位置处过采样，进而造成采样分布不均匀，所以要在插值后的数据基础 上进行均匀分布。均匀分布的过程为：两端数据设置为内插后的数据，对于中间的数据采用 Lagrange插值公式进行重新插值，进而可以实现均匀分布，解决某些位置处过采样的问题。
4. 根据权利3所述的心电数据快速重采样方法，其特征在于：对于原始采样频率值大 于目标采样频率的情况，在进行分数倍插值后，达到大于原始采样频率值的目标采样频率 的最小整数倍值，然后就需要对插值后的数据进行整数倍抽取，从而实现到目标采样频率 的转换；而对于原始采样频率小于目标采样频率的情况，为保证采样比率小于1，以便进行 分数倍的插值，要先对数据进行整数倍内插，达到小于目标采样频率的最大整数倍（此值 为1时可省略），再采用分数倍内插达到目标采样频率值。
5. 根据权利1、2所述的心电数据快速冲采样方法，其特征在于：在符合采样定律的前 提下，可以实现任意频率值的转换。首先判断原始频率与目标频率的大小，然后分情况进行 处理，其核心是通过分数倍的内插，减小运行时间及占用存储空间。当然该方法也考虑特殊 情况，如减采样中原始采样频率恰好为目标采样频率的整数倍，那么无须再进行插值，可以 直接通过整数倍抽取来实现；而增采样过程中，若目标频率证号为原始频率的整数倍，那么 在进行整数倍插值之后，即可达到目标频率，无须进一步分数倍插值。
【文档编号】A61B5/0402GK104055512SQ201410314071
【公开日】2014年9月24日 申请日期:2014年7月1日 优先权日:2014年7月1日
【发明者】张跃, 王佳丽, 黎洪 申请人:深圳市岩尚科技有限公司, 清华大学深圳研究生院