一种对心电数据精度转换的方法

一种对心电数据精度转换的方法

技术领域

心电数据精度的转换方法属于采样控制系统领域，主要适用于心电数据的重采样，可将数据进行精度转换，得到不同采样精度的数据，用于后续处理及研究；该发明同样适用于其它类似的数据精度转换场合。 

背景技术

人体心电信号的主要频率范围为0.05～100Hz，而电压幅度为0～4mV，信号十分微弱。而且，心电信号中通常混杂有其他生物电信号等干扰，使得测量环境较差。为了能够检测到不失真的心电信号，以用于后续临床处理，因而大部分心电检测设备在硬件电路上均要使用前置放大模块等对心电原始信号进行放大。 

另外，由于人体心电信号为模拟信号，为了便于传输及存储，必须要转换为数字信号。因而心电检测设备需要将采集到的心电模拟信号通过A/D转换，得到数字信号。由于心电信号的电压幅度过小，处于mV级别，而大部分A/D转换模块中要求输入信号的幅度要在单位V的级别上，因而心电检测设备除去前面介绍的前置放大电路，还需要采用进一步放大以及偏移值处理，以达到A/D模块中对输入信号的电压级别要求。其中的增益值、偏移值等参数，必须要同检测到的心电数字信号一起提供给心电数据处理的软件，才能进行后续心电滤波及分析的处理。 

而当前市面上心电检测设备层出不穷，种类更是多种多样，在心电采样的参数中，如采样频率、增益值、偏移值、采样精度等各不相同，这导致了一种心电设备检测的数据不能使用另一种心电分析软件对其进行分析，对于心电数据的统一管理分析非常不利。目前清华大学深圳研究生院嵌入式系统实验室正在研发的“人体健康状态监测分析与信息管理云服务网络平台”，通过自主研发产品以及整合集成医院、社区、家庭现有的各种医疗监护设备，能够随时随地监测用户的心电、血压等生理参数，通过蓝牙、WIFI、有线网络等，实时发送给医用智能手机和医院服务器多参数监护中心，以供医护人员进行心电、血压等参数的特征提取与分析，为病人提供及时的诊疗服务。由于当前各种心电检测设备的参数不同，就需要在对数据进行处理之前，先将心电数据按照心电分析软件处理的数据格式，将原始心电数据中采样频率、增益值、采样精度等转换为目标数据所规定的值。本发明中，针对的是心电数据采样精度的转换，其中也涉及到使用原始及目标数据规定的增益值、偏移值等参数。 

本发明综合考虑了项目的数据处理需要及国际化标准的要求，对所取得的数据进行精度转换设置的标准为：单导联、三导联的数据采用实际项目应用中1字节的数据精度，12导联的心电数据采用MIT-BIH数据库规定的2字节精度。 

发明内容

本发明的目的在于提供一种心电数据采样精度的转换方法，可使得不同精度的心电数据之间进行精度转换，有利于对于心电数据的分析与统一管理，同样也适用于其他离散数据精度转换的类似应用中。 

本发明使用的相关理论包括：心电数据前期处理、数据的左右移位。 

本发明中涉及的心电数据的精度转换方法与上述技术的不同之处，在于将传统的数据左右移位改变数据精度，与心电数据处理的前期过程相结合，使得心电数据在精度转换的过程中，误差大幅减小。本发明的方法流程图如图1所示，可以实现不同精度的心电数据之间的转换，根据目标精度与原始数据精度的大小，分情况进行。 

在进行数据精度转换之前，需要获知原始数据及目标数据的精度、增益值与偏移值，不同心电检测设备的这些参数一般不会完全相同。另外，心电数据均被保存为二进制文件。 

(1)若目标数据精度大于或等于原始数据，则直接将心电数据前期处理过程应用于转换，即可将原始心电数据文件转换为目标心电数据文件。 

(2)若目标数据精度小于原始数据，则将心电数据前期处理过程与传统的移位方法相结合，这样可以在保证误差较小的前提小，尽量提高执行效率，完成心电数据的精度转换，生成目标要求的心电数据文件。 

上述两种情况的心电数据精度具体的转换步骤，见本说明书的第五部分：具体实现方式。 

本发明之所以考虑目标数据精度等于原始数据精度的转换情况，是因为即使两者表示心电数据的字节数相等，即精度相同，但是由于增益值、偏移值等参数的不同，也会造成两者并不能统一由其中一种心电数据格式对应的心电分析软件来进行病历分析等操作。因而本发明所涉及的方法，将这一情况与目标数据精度大于原始数据的情况合并，只需要按照心电数据处理的前期过程进行转换，便可完成。 

附图说明

图1：本发明提出的方法流程图； 

图2：精度转换原理图。 

具体实现方式 

本发明所涉及的心电数据精度转换方法，方法的流程图如图1所示。只要获知原始心电数据文件与目标心电数据文件各自所规定的偏移值、增益值，便可以实现任意精度之间的转换，具体的心电数据处理过程及精度转换的原理图如图2所示。 

在说明具体方式之前，首先对传统的精度转换方法进行分析。传统的精度转换方法，是直接采用左右移位的方式来实现。例如，若原始心电数据的精度为12位，目标数据精度为16位，那么直接采用向左4位，即实现数据的精度转换。但是这种方法，并未考虑数据的增益值、偏移值等参数，若原始心电数据的偏移值与目标数据的相等，并设其为0，目标数据为0x3ff0，在计算心电数据的浮点数值时，(0x3ff-0)与(0x3ff0-0)在除以各自的增益值之后，得到浮点数类型的心电信号数据，误差只发生在后4位的补齐上。但如果两者的偏移 值、增益值不等，特别是偏移值的不等，那么会造成更大的误差，如设原始数据偏移值为0，目标数据的偏移值为50，那么数据在计算时可以看到由于目标数据的增益值并没有按照左移4位的方式进行放大，造成(0x3ff-0)与(0x3ff0-50)两者的误差较大。同样在目标数据精度小于原始数据精度的转换过程中，也存在这样的问题，这也说明了为什么两种心电数据的采样精度相同，但是在偏移值等参数不同时，仍需要进行精度转换。 

由于本发明在进行实现时，讨论了两种情况，因而分两个例子来说明具体实现方式。以MIT-BIH其中一种格式的心电数据来说明，这种格式的心电数据采用1.5个字节来表示一个数据，即精度为12位。那么暂不考虑其具体的导联数，先按照本发明在权利要求书特征4中提到的两种格式对其进行格式转换。本发明在统一格式时，规定单导联和3导联的心电数据精度统一为1字节，即8位；12导联的心电数据精度统一为2字节，即16位。对于精度相同的情况，由于其采用的方式与目标数据精度大于原始精度的情况，这里不再赘述。那么该例子可以充分讨论这两种情况，具体说明实现的步骤。 

为了简化问题，描述得更加清晰，下列步骤中采用其中1个数据来说明，这个数据通过原始文件读取，并设获取到的原始、目标心电数据值分别为srcValue、destValue，原始、目标心电数据增益值分别为srcGain、destGain，而原始、目标心电数据的偏移值分别为srcOffset、destOffset；关系符合与原理对应的等式：(destValue-destOffset)/destGain＝(srcValue-srcOffset)/srcGain。 

1)目标数据精度大于原始数据精度：12位到16位 

这种情况可直接按照心电数据处理的前期过程进行处理，根据图2所示的精度转换原理图，我们为了提高运行速度，对其进行了简化运算的变换，可得步骤如下。 

(1)取得的整型原始数据值减去偏移值； 

temp＝srcValue-srcOffset； 

(2)计算目标数据与原始数据的增益比值k，k为浮点数； 

k＝destGain/srcGain； 

(3)将第一步中计算的结果与第二步中增益比值相乘； 

tempResult＝temp*k； 

(4)将上一步中的结果在计算机程序中强制转换为整型数据； 

destValue＝(int)tempResult； 

(5)上一步得到的整型数据，加上目标数据对应的偏移值destOffset，即可完成该数据的精度转换； 

destValue＝destValue+destOffset； 

如上步骤，对原始数据文件中的逐个数据进行转换，保存为目标数据精度类型的文件，即可完成。本方法严格遵循图2中的原理，只是对其进行了改动，以减少数据运算的步骤，保证等式(destValue-destOffset)/destGain＝(srcValue-srcOffset)/srcGain的成立。 

2)目标数据精度小于原始数据精度：12位到8位 

当目标数据精度小于原始数据精度时，可采用传统移位方法与心电处理前期过程相结合 的方法，在减小误差的同时，尽量提高运行效率。同样，在计算机实现时，对处理过程进行了简化，以减少运算步骤，具体如下。 

(1)获取整型原始心电数据，减去原始偏移值； 

temp＝srcValue-srcOffset； 

(2)上一步的结果先进行右移，右移4位，使得数据结果为8位，与目标心电数据表示的位数相同； 

temp＝temp＞＞4； 

这一步骤相当于除以2^4，即16，从而去掉第一步结果中的后4个二进制位； 

(3)计算参数k； 

k＝destGain/srcGain*(2^4)； 

为了使等式(destValue-destOffset)/destGain＝(srcValue-srcOffset)/srcGain成立，由于第2步中使用了右移步骤，相当于原始心电数据端先除以因子(2^4)，所以在本步骤中要对其进行补偿。 

(4)将第2步的结果temp与上一步计算的参数k相乘； 

tempResult＝temp*k； 

(5)同样，将第4步中的结果强制转换为整型数据； 

destValue＝(int)tempResult； 

(6)上一步的结果加上目标心电数据所对应的偏移值destOffset，即完成该数据的精度转换； 

destValue＝destValue+destOffset； 

同样，对原始数据文件中的各个数据进行如上步骤，将最终计算的结果同样以二进制文件的形式保存，便可完成心电数据精度的转换。在第2步中，将数据对数据进行移位，可以将原先1.5个字节的数据转换为1个字节的数据，使得后续参与乘法运算的因子字节数减少，可以提高运行效率；并且，1.5字节的数据读取方式复杂，需要专门的程序处理，不如整个字节读取写入方便。 

综上所述，本发明提出的心电数据精度转换方法，可以实现心电数据间任意精度间的转换，下面讨论一般情况。设原始心电数据的精度为M位，目标心电数据的精度为N位，同样对其数值、增益值、偏移值进行假设，其值与上述具体例子中的相同，原始与目标数据的值分别为srcValue、destValue；增益值分别为srcGain、destGain；偏移值分别为srcOffset、destOffset；它们之间的关系仍符合与原理对应的等式：(destValue-destOffset)/destGain＝(srcValue-srcOffset)/srcGain。。 

根据原始心电数据与目标数据的精度大小分情况进行讨论，同样使用其中一个数据的转换步骤来进行说明，这里将步骤进行了简化。 

1)N＞＝M的情况 

(1)原始数据值减去原始偏移值； 

temp＝srcValue-srcOffset； 

(2)计算参数值k，结果为浮点型； 

k＝destGain/srcGain； 

(3)将第1步的结果与第2步计算的参数k相乘； 

tempResult＝temp*k； 

(4)将上一步的结果强制转换为整型数据； 

destResult＝(int)tempResult； 

(5)上一步结果加上目标心电数据对应的偏移值，即完成了该心电数据的精度转换； 

destResult＝destResult+destOffset； 

对表示原始心电数据格式的文件中的所有数据，进行上述步骤的运算，最终结果保存为二进制文件，即为目标心电数据格式。 

2)N＜M的情况 

(1)原始数据减去原始偏移值； 

temp＝srcValue-srcOffset； 

(2)第1步计算的结果进行移位，右移(M-N)位，其效果等价于除以2^(M-N)； 

temp＝temp＞＞(M-N)； 

(3)计算参数k，同样对第2步进行参数补偿，以保证符合原理； 

k＝destGain/srcGain*2^(M-N)； 

(4)第2步骤所得的结果temp与第3步计算的参数k相乘； 

tempResult＝temp*k； 

(5)将上一步所得的结果进行强制转换，转换为整型数据； 

destValue＝(int)destValue； 

(6)第5步的结果加上目标心电数据格式对应的偏移值，即可完成该情况下心电数据精度的转换； 

destValue＝destValue+destOffset； 

同样，这只是一个数据的转换过程，对文件中所有数据进行如上计算，便最终可完成对数据的精度转换。 

这里对本发明中涉及的心电数据的精度转换方法进行了更加一般情况的阐释，上述步骤为了表达的更加清楚，步骤较多，在实际计算机实现中，可以将其中某两个或两个以上的步骤进行合并，减少程序的长度。如M＜N的情况下，可以将第1步和2步合并，将计算得到的结果直接进行移位；也可以将前4步合并为一个步骤，这样可以大幅地减小程序长度，等等。 

根据本发明对于心电数据精度转换方法的描述，可以发现其它类似应用场合下的数据文件，只要将其实际数据处理的前期过程与传统的移位方法相结合，同样可以用这种类似方法来实现精度之间的转换。这对于不同精度下二进制文件中数据的格式统一，以及数据的压缩等，都具有非常重要的意义。