发明名称 基于PT变换与线性预测结合的心电图信号无损压缩方法
摘要 本发明公开了一种基于PT变换与线性预测结合的心电图信号无损压缩方法,属于数据处理技术领域。其压缩过程是:对原始的心电图信号进行PT变换;将变换后的原始信号值与估计值相减,得到整段心电图信号的残差信号;对残差信号进行自适应可变阶RICE编码输出编码比特流,完成对心电图信号的无损压缩;对压缩的比特流进行解码;重构信号,完成解压缩。本发明作为一种无损压缩方法,保证了所压缩心电图信号的信息完整及正确性,可用于心电图信号的传输和存储。
申请公布号 CN101669819B 申请公布日期 2011.06.01
申请号 CN200910024050.2 申请日期 2009.09.25
申请人 西安电子科技大学 发明人 吴家骥;高力鑫;焦李成;石光明;张向荣;侯彪;公茂果;马文萍
分类号 A61B5/0452(2006.01)I;G06F19/00(2006.01)I 主分类号 A61B5/0452(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 1.一种基于PT变换与线性预测结合的心电图信号无损压缩方法,包括如下步骤:(1)对原始的心电图信号进行PT变换;(2)对变换后的信号按如下过程进行预测;(2a)选取一段信号根据下式进行训练,得到固定预测系数组:<maths num="0001"><![CDATA[<math><mrow><mover><mi>a</mi><mo>&RightArrow;</mo></mover><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>C</mi><mi>T</mi></msup><mi>C</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>C</mi><mi>T</mi></msup><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FSB00000432767300012.GIF" wi="28" he="47" />为预测系数组的向量,C为估计矩阵,<img file="FSB00000432767300013.GIF" wi="29" he="55" />为所选取信号段的向量;(2b)根据1)式对当前点X<sub>n</sub>之前的2k个点进行训练,得到动态预测系数组,k为预测阶数;(2c)根据设定的门限选择采用动态预测系数组或固定预测系数组,当1)式中的矩阵(C<sup>T</sup>C)<sup>-1</sup>的行列式值低于门限阈值时,则选取固定预测系数组作为预测系数,反之,采用当前的动态预测系数组作为预测系数;(2d)根据确定的预测系数,利用<u>如下</u>自回归公式,求得当前点X<sub>n</sub>的估计值:<maths num="0002"><![CDATA[<math><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>en</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>a</mi><mi>t</mi></msub><msub><mi>X</mi><mrow><mi>n</mi><mo>-</mo><mi>t</mi></mrow></msub><mo>|</mo><mo>|</mo><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>n</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>a</mi><mi>t</mi></msub><msub><mi>X</mi><mrow><mi>e</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mi>t</mi><mo>)</mo></mrow></mrow></msub><mo>|</mo><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,X<sub>en</sub>是当前点X<sub>n</sub>的估计值,k为预测阶数,a为确定的预测系数,X<sub>n-t</sub>为X<sub>n</sub>的第前t个点,X<sub>e(n-t)</sub>为X<sub>en</sub>的第前t个点;(2e)重复(2b)-(2d),得到整段信号的估计值,完成对变换后心电图信号的预测;(3)将变换后的原始信号值与估计值相减,得到整段心电图信号的残差信号;(4)对残差信号进行如下自适应可变阶RICE编码;(4a)设当前点e<sub>n</sub>,根据e<sub>n</sub>之前的n个点的数值大小的分布,求得这n个点数值大小绝对值的均值;设定以2<sup>2</sup>为门限,每2<sup>2</sup>为一阶,根据均值所在的门限,选择e<sub>n</sub>的编码阶数;(4b)根据如下定义,确定e<sub>n</sub>的前码段和阶数码段,其中k为编码阶数:前码段:e<sub>n</sub>整除2<sup>k</sup>,所得商为n<sub>a</sub>,则以n<sub>a</sub>个1和一个0作为前码段;阶数码段:e<sub>n</sub>整除2<sup>k</sup>,所得余数为n<sub>e</sub>,则以k个n<sub>e</sub>的二进制码作为阶数码段;(4c)将前码段和阶数码进行合并,完成e<sub>n</sub>的编码;(4d)重复(4a)-(4c)完成对残差信号的编码;(4e)输出编码比特流,完成对心电图信号的无损压缩。
地址 710071 陕西省西安市太白路2号