《一种心电数据快速重采样方法.pdf》由会员分享，可在线阅读，更多相关《一种心电数据快速重采样方法.pdf（12页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104055512A43申请公布日20140924CN104055512A21申请号201410314071922申请日20140701A61B5/040220060171申请人深圳市岩尚科技有限公司地址518000广东省深圳市南山区蛇口南海大道科技大厦辅楼511号申请人清华大学深圳研究生院72发明人张跃王佳丽黎洪54发明名称一种心电数据快速重采样方法57摘要心电数据采样频率转换的传统方法是，先对原始信号进行整数倍插值，达到原始频率与目标采样频率的最小公倍数，再对插值后信号进行整数倍抽取，从而得到目标采样频率要求的数据。但是传统方法中存在的严重问题是，当插值倍数较大时，运行速。
2、度很慢，占用存储空间过大，达不到实际应用的要求。本发明为一种心电数据快速重采样方法，解决了这个问题。快速重采样方法的特征在于，在达到工程要求以及符合采样定律的前提下，实现心电数据的任意采样频率转换。这种方法可以很好地应用于实际工程使用，而且对于其它要满足工程需求的信号重采样也非常有意义。51INTCL权利要求书1页说明书6页附图4页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书6页附图4页10申请公布号CN104055512ACN104055512A1/1页21心电数据快速重采样方法，其特征在于使用分数倍内插的内插方法，从而减少运行时间及占用存储空间，使得在达到工程应用要。
3、求以及符合采样定律的前提下，可快速实现心电数据的采样频率任意转换。方法主要包括以下几个部分；1分数倍内插通过对原始数据进行分数倍内插，达到目标采样频率或者其大于原始采样频率的最小整数倍；2均匀分布由于采用分数倍的插值方法，造成不均匀采样的问题，所以要在插值之后，在原始数据及插入的数据基础上，对其进行重新均匀分布；3整数倍抽取与内插若原始采样频率大于目标采样频率，在完成内插之后，要对数据进行整数倍抽取，达到目标采样频率；若原始采样频率小于目标采样频率，那么需要先对数据进行整数倍的内插，使之达到小于目标采样频率的最大整数倍值此值为1时可省略，后续采用分数倍内插达到目标采样频率值。2根据权利1所述的。
4、心电数据快速重采样方法，其特征在于采用分数倍内插的方法，直接达到目标频率或其大于原始频率的最小整数倍，而不是传统方法中整数倍内插到原始频率与目标采样频率的最小公倍数。主要包括以下内容1插值比率新方法中的插值比率为分数，且小于1，其意义为每几个数值需要插入一个新值；2内插过程原始数据读取到一个数组中，使用计算得到分数倍的插值比率，在插值过程中，将原始数据数组下标与插值比率相乘，获取其整数部分的值，当整数部分增长1时，即在当前的下标处，插入一个数据，其数值等于数组在当前位置的值。3根据权利2所述的心电数据快速重采样方法，其特征在于对原始数据进行分数倍内插后，会造成某些位置处过采样，进而造成采样分布。
5、不均匀，所以要在插值后的数据基础上进行均匀分布。均匀分布的过程为两端数据设置为内插后的数据，对于中间的数据采用LAGRANGE插值公式进行重新插值，进而可以实现均匀分布，解决某些位置处过采样的问题。4根据权利3所述的心电数据快速重采样方法，其特征在于对于原始采样频率值大于目标采样频率的情况，在进行分数倍插值后，达到大于原始采样频率值的目标采样频率的最小整数倍值，然后就需要对插值后的数据进行整数倍抽取，从而实现到目标采样频率的转换；而对于原始采样频率小于目标采样频率的情况，为保证采样比率小于1，以便进行分数倍的插值，要先对数据进行整数倍内插，达到小于目标采样频率的最大整数倍此值为1时可省略，再采。
6、用分数倍内插达到目标采样频率值。5根据权利1、2所述的心电数据快速冲采样方法，其特征在于在符合采样定律的前提下，可以实现任意频率值的转换。首先判断原始频率与目标频率的大小，然后分情况进行处理，其核心是通过分数倍的内插，减小运行时间及占用存储空间。当然该方法也考虑特殊情况，如减采样中原始采样频率恰好为目标采样频率的整数倍，那么无须再进行插值，可以直接通过整数倍抽取来实现；而增采样过程中，若目标频率证号为原始频率的整数倍，那么在进行整数倍插值之后，即可达到目标频率，无须进一步分数倍插值。权利要求书CN104055512A1/6页3一种心电数据快速重采样方法技术领域0001心电数据快速重采样方法属于。
7、采样控制系统领域，主要适用于心电数据的重采样，得到不同采样频率的数据，用于进一步研究；对于其它有同样需求的离散数据重采样同样适用。背景技术0002改革开放以来，我国医院信息化建设得到了长足发展，特别是医院信息系统HOSPITALINFORMATIONSYSTEM，HIS对于提高医院科学管理水平，以及推动医院现代化建设发挥了重要作用。而实现信息标准化以完成信息的统一管理，是目前医院HIS建设中面临的首要问题。0003当前实现远程传输的数据主要有心电、体温、血压、血氧饱和度等，这些数据信息需要实现标准管理。在上述数据的格式统一过程中，最复杂的当属心电数据，心电数据检测设备多样化，不同设备的数据格式。
8、并不相同，而且心电数据量大，对于诊断病人病情至关重要，所以心电数据格式的统一势不容缓。0004心电数据格式转换中的信息包括电设备的导联数、病人信息、增益、采样频率等，而由于心电分析软件大都采用差值运算或变换的方法，来对心电原始数据进行预处理及分析，因而增益等其他参数对于心电诊断分析的影响不大；但是采样频率如果与目标分析软件不一致，那么心电波形的显示将会出现错误，甚至造成无法诊断或诊断错误。0005另外，当前研究心电方法等都是基于国际标准心电数据库进行研究的，当前全球公认的标准心电数据库有三个，分别是美国麻省理工学院提供的MITBIH数据库，美国心脏学会的AHA数据库以及欧洲ATT心电数据库。其。
9、中，MITBIH数据库的应用相对广泛，MITBIH数据库的采样频率固定为360HZ，而国内很多心电检测设备的采样频率却有很多变化，将标准数据库的采样频率进行转换，可以采用其来完成特定频率下的分析，从而可以实现与标准的对比，这对于校准非常重要。0006传统实现采样频率转换的方法即严格按照离散时间信号的内插与抽取理论，先对待转换信号进行整数倍的插值，达到原始采样频率与目标采样频率值的最小公倍数，然后再对插值后的数据进行整数倍抽取。例如若将心电数据从MITBIH数据库的固定采样频率360HZ变为250HZ，传统方法需要先将数据内插到9000HZ，即内插了25倍，然后再进行抽取36倍，这样虽然保证了重。
10、采样的数据完整，但是运行效率很低，占用存储空间过大。发明内容0007本发明的目的在于提供一种心电数据的快速重采样方法，对于心电格式转换及标准化，以及心电检测设备的校验都有非常重要的作用，同样也适用于有类似需求离散数据重采样的实际应用。0008本发明使用的相关理论包括采样定律，离散时间信号的内插与抽取整数倍，LAGRANGE插值公式。说明书CN104055512A2/6页40009我们所涉及的心电数据快速重采样与上述技术不同之处，在于将内插的倍数改为分数倍内插，然后再使用LAGRANGE插值公式进行重新均匀分布。新方法流程图如图1所示，可以实现不同采样频率间的转换，根据原始采样频率值与目标采样频。
11、率值的关系，分两种情况进行。00101若原始采样频率值小于目标采样频率值，步骤如下0011步骤1，先对数据进行整数倍内插，达到小于目标采样频率的原始值最大整数倍此值为1时可省略，设该值为NN1，那么在原始数据相邻两个数据之间插入N1个数据，数据值为两个数据中前者的值；0012步骤2，在完成整数倍内插后，就可以对其进行分数倍内插，得到采样频率为目标频率值的数据；0013步骤3，最后对分数倍内插后的结果通过LAGRANGE插值公式进行重新均匀分布。00142若原始采样频率值大于目标采样频率值，步骤如下0015步骤1，先对数据进行分数倍内插，达到大于原始采样频率值的目标值的最小整数倍，设该值为M，M。
12、必定大于1；0016步骤2，对分数倍内插的结果通过LAGRANGE插值公式进行重新均匀分布；0017步骤3，最后对结果进行整数倍抽取，抽取倍数为M倍，即从相邻M个数据中，抽取其中的第一个值作为结果，便得到采样频率为目标值的数据。0018本发明中先对数据进行分数倍内插，即将原先每相邻两个数据之间插入数据改为，每相邻多个数据插入一个数据；对于分数倍内插造成的不均匀分布，采用对数据进行重新插值计算的方式消除某些位置过采样的问题。插值并进行均匀分布的示意图如图2所示，而普通插值未均匀分布的示意图如图3所示，均匀分布的效果图，如图4所示。0019当然也要考虑到特殊情况原始采样频率为目标采样频率的整数倍，。
13、以及目标采样频率正好为原始频率值的整数倍。前者在进行插值到大于原始频率的目标频率值的最小整数倍之前，经过判定之后，原始频率正好为目标频率的整数倍，则无需进行插值，可以直接进行整数倍抽取，达到目标频率。而后者，在插值到小于等于目标采样频率的原始频率的最大整数倍之后，如果目标采样频率为原始频率的整数倍，那么经过插值后，达到的频率即为目标频率，无需继续分数倍插值及均匀分布。附图说明0020图1本发明提出的方法流程图；0021图2插值后均匀分布示意图；0022图3未执行均匀分布效果图；0023图4均匀分布效果图。0024具体实现方式0025本发明中，可以实现任意频率间的转换，但是如果要能正常显示，必须。
14、要符合采样定理，即采样频率必须大于被采样数据最大频率的2倍。0026以原始采样频率为360HZ的MITBIH数据库中的心电数据转换为其他采样频率来说明，本发明的采样频率转换分为两种情况，设置目标采样频率值分别为250HZ和800HZ，从两种情况更加清晰的说明方法的实现方式，在这里进行说明的数据都以一维数组的形式说明书CN104055512A3/6页5存放。特殊情况不再进行详细说明，如将360HZ减采样转换为180HZ，或者增采样为720HZ，前者只需要进行2倍抽取，后者只需要进行2倍内插，即可以达到目的频率值，且都可以归到一般情况进行讨论，不再赘述。00271原始采样频率大于目标采样频率360。
15、HZ到250HZ0028为了简化问题，取1S的数据大小为例，且心电数据为单导联，数据为8BIT的无符号整型数据。设原始数据数组为SRC360，分数倍内插的结果数组为MID500，最终结果的数组为RESULT250；00291对数据进行分数倍插值是之达到大于原始值的目标采样频率值的最小整数倍，此例中该值为500HZ，即先将SRC进行分数倍插值，达到采样频率为500HZ的数据MID500。0030插值比率RATE500360/3607/18；使用原始数据数组的下标值与此插值比率相乘，若所得结果的整数部分比前一项的值增长1，那么就在该位置插入一个新数据，数据值为原始数据数组在前一位置的值。例如从数组。
16、下标位置0到2，与插值比率7/18相乘的结果，其整数部分一直为0，而下标3与插值比率7/18相乘，其整数部分为1，则在新数组中下标3的位置插入一个新数据，其值为原始数据数组下标为2处的值。0031以变量POS1表示SRC数组的下标位置，COUNT表示当前已经插入新数据的个数，初始值为0；0032当POS10，1，2时，INTPOS1RATE一直为0，因而MID的值等于相应位置的SRC的值；0033MIDPOS1COUNTSRCPOS1；0034当POS13时，INTPOS1RATE1；整数部分增长1，因而需要插入新数据；0035MIDPOS1COUNTSRCPOS11；0036COUNTCOU。
17、NT1；0037MIDPOS1COUNTSRCPOS1；0038依次类推，POS14，5时，INTPOS1RATE一直为1，那么不进行插值；0039MIDPOS1COUNTSRCPOS1；0040POS16时，INTPOS1RATE2，整数部分增长1，则此处需要插入新数据；0041MIDPOS1COUNTPOS1SRCPOS11；0042COUNTCOUNT1；0043MIDPOS1COUTNPOS1SRCPOS1；0044按照这种方式，直到SEC的数组被遍历完毕。00452均匀分布对于分数倍插值后的结果进行均匀分布，消除某些位置过采样而造成不均匀分布的问题。采用的方式为当上个步骤中计算的结果。
18、需要插入新数据时，以整数部分前一次增长下标的位置为初始端，新插值位置的下标为末端，对这其中的数据进行重新分布。0046以上一步骤中的例子来说明；0047POS13时，INTPOS1RATE增长为1，此时COUNT0；0048设初始端点与末端点的变量为START，END；那么设置初始端点为0，末端位置为POS1即为3；对MID数组中03的数据进行重新均匀分布，两端位置不变，中间位置根说明书CN104055512A4/6页6据LAGRANGE插值公式进行重新计算；0049START0；ENDPOS1COUNTPOS13，COUNT0；0050MID0、MID3不变；00510052这里N3，带入上。
19、述公式可得00530054从后端到前端使用该公式重新计算，将新数据值填入MIDI中，从而实现均匀分布；以此类推，在进行插值之后就要记录均匀分布块的初始端点；0055STARTPOS1COUNT；此时POS13，COUNT1；0056POS16时，INTPOS1RATE增长为2，此处进行插值；0057ENDPOS1COUNT；0058保持两端点值不变，MIDSTARTMIDSTART，MIDENDMIDEND；0059中间值从后端到前端利用上述LAGRANGE插值公式进行重新计算，N同样为3；00600061如此类推，在每次插值之后，都对当前区段进行重新计算，防止过采样，直到SRC数组遍历完毕。。
20、00623整数倍抽取0063上述两个步骤完成之后，得到采样频率为500HZ的心电数据，表示形式为一维数组MID500，而目标采样频率为250HZ，因而需要进行整数倍抽取。0064此例当中，抽取倍数N500/2502；需进行2倍的抽取，采用的方式为对MID数组中两个相邻的数据，只取其一，方法中将数据设置为两个数据中位置在前的一个；0065RESULTIMID2I；I0I0I359；0073MID718SRC359；说明书CN104055512A5/6页70074MID719SRC359；0075由此完成整数倍的内插，获得中间结果MID720；00762分数倍内插达到目标采样频率与上一种情况中分数。
21、倍内插的方式相同，将MID720进行分数倍内插，达到目标采样频率为800HZ的结果数组RESULT800；0077插值比率RATE800720/7201/9；使用原始数据数组的下标值与此插值比率相乘，若所得结果的整数部分比前一项的值增长1，那么就在该位置插入一个新数据，数据值为原始数据数组在前一位置的值，采用的方法与上一例子中相同。0078以变量POS1表示MID数组的下标位置，COUNT表示当前已经插入新数据的个数，初始值为0；0079当POS108时，INTPOS1RATE一直为0，因而MID的值等于相应位置的MID的值；0080RESU1TPOS1COUNTMIDPOS1；0081当PO。
22、S19时，INTPOS1RATE1；整数部分增长1，因而此处需要插入新数据；0082RESU1TPOS1COUNTMIDPOS11；0083COUNTCOUNT1；0084RESULTPOS1COUNTMIDPOS1；0085依次类推，POS1917时，INTPOS1RATE一直为1，那么不进行插值；0086RESULTPOS1COUNTMIDPOS1；0087POS118时，INTPOS1RATE2，整数部分增长1，则此处需要插入新数据；0088RESULTPOS1COUNTPOS1二MIDPOS11；0089COUNTCOUNT1；0090RESULTPOS1COUNTPOS1MIDPOS。
23、1；0091按照这种方式，直到MID的数组被遍历完毕。00923均匀分布同样与之前例子中所进行的均匀分布类似，在进行分数倍内插之后，同样需要对结果进行重新计算，以解决过采样的问题。采用的方式与上一例子中相同，设置初始端点为原先插入数据的位置，末端为新插值位置的下标，两端位置的数值保持不变，中间端点根据LAGRANGE插值公式进行重新分布。0093同样以上一步骤中的例子来说明；0094POS19时，INTPOS1RATE增长为1，此时COUNT0；0095设初始端点与末端点的变量为START，END；那么设置初始端点为0，末端位置为POS1即为9；对RESULT数组中09的数据进行重新均匀分布，。
24、两端位置不变，中间位置根据LAGRANGE插值公式进行重新计算；0096START0；ENDPOS1COUNTPOS19，COUNT0；0097RESULT0、RESULT9不变；00980099这里N9，带入上述公式可得说明书CN104055512A6/6页801000101从后端到前端使用该公式重新计算，将新数据值填入RESULTI中，从而实现均匀分布；以此类推，在进行插值之后就要记录均匀分布块的初始端点；0102STARTPOS1COUNT；此时POS19，COUNT1；0103POS118时，INTPOS1RATE增长为2，此处进行插值；0104ENDPOS1COUNT；0105保持两。
25、端点值不变，RESULTSTARTRESULTSTART，RESULTENDRESULTEND；0106中间值从后端到前端利用上述LAGRANGE插值公式进行重新计算，N同样为3；01070108如此类推，在每次插值之后，都对当前区段进行重新计算，防止过采样，直到RESULT数组遍历完毕。0109经过以上3个步骤，完成了从360HZ到目标频率800HZ的转换，得到的结果为一维数组RESULT8000110通过对以上两种情况的分别说明，对心电数据快速重采样方法进行了说明，传统方法中，对于360HZ转换为250HZ，要先对原始数据数组进行插值，达到两者最小公倍数9000HZ，然后再采用抽取36倍，。
26、达到目标采样频率为250HZ的结果；对于360HZ转换为800HZ，要先进行整数倍的插值，达到7200HZ，随后再进行9倍的抽取，达到目标频率800HZ。例子中仅给出1S钟时间的简化情况，而在实际应用中，如动态心电图HOLTER记录的数据，其数据量通常为大约24小时，那么采用传统方法运行效率明显很低，而且所占用存储空间也大。0111本发明心电数据快速重采样方法中，在前面介绍中可以看到，从360HZ转换到250HZ，数据插值的倍数小于1倍，达到500HZ，随后抽取2倍，中间结果占用存储空间为传统方法的1/18；而从360HZ转换到800HZ，数据先插值2倍，达到720HZ，然后再进行小于1倍的插值，达到800HZ的要求，中间结果占用存储空间为传统方法的1/10。虽然本发明心电数据快速重采样方法还需对分数倍插值的结果进行均匀分布，但是相较传统方法，无论在运行效率还是运行期间占用空间都得到了很大的改善。说明书CN104055512A1/4页9图1说明书附图CN104055512A2/4页10图2说明书附图CN104055512A103/4页11图3说明书附图CN104055512A114/4页12图4说明书附图CN104055512A12。