PhysioNet生理信号处理(三)--WFDB Toolbox for Matlab

PhysioNet是一个由美国国家通用医学科学院(NIGMS)和国家生物医学成像和生物工程院(NIBIB)共同支持的面向科研人员免费提供生理信号库和其处理工具的一个资源网站或者项目。目前项目内容主要分为两个:PhysioBank和PhysioToolkit,也就是生理信号库和处理工具箱。

网站:https://physionet.org/

上一篇我们对整个PhysioToolkit的组织架构和其中一些软件包介绍了,相信小伙伴们迫不及待的要磨拳檫掌,跃跃欲试了吧,这一期就来实际上手。

一、Toolbox的安装

上一期也提到这个,可以将.m函数下载下来添加到Matlab的路径中即可。现在还有一个更简单的办法,在Matlab中直接输入一下内容:

[old_path]=which('rdsamp'); if(~isempty(old_path)) rmpath(old_path(1:end-8)); end

wfdb_url='https://physionet.org/physiotools/matlab/wfdb-app-matlab/wfdb-app-toolbox-0-10-0.zip';

[filestr,status] = urlwrite(wfdb_url,'wfdb-app-toolbox-0-10-0.zip');

unzip('wfdb-app-toolbox-0-10-0.zip');

cd mcode

addpath(pwd)

savepath

等待一会儿,WFDB的工具箱将直接下载在Matlab的当前路径中,并产生新的文件夹:

wfdb-app-toolbox-0-9-9

二、测试安装

安装好后我们就可以测试一下运行起来是否有障碍。WFDB给我们提供了两个测试方法:

方法一:运行内置的Demo:

wfdbdemo

方法二:读入一个MIT数据库的生理信号并显示:

[tm,sig]=rdsamp('mitdb/100',1);

plot(tm,sig);

如果两个运行都有figure结果而且没有报错,那么说明工具箱成功安装,鼓掌欢呼!

三、函数说明

前面提到过,WFDBToolbox for Matlab一共有20个函数,我数了一下其实有40多个。因为工具箱里面除了WFDB函数,还有其他工具箱(如ecgpuwave等)中的函数,这对于一个工具箱来说并不算多,我们可以一个个的介绍。

函数列表可以在WFDB Applications Guide中查看:https://physionet.org/physiotools/wag/wag.htm

1.读入生理信号信号及其注释(rdsamp)

[tm,signal,Fs]=rdsamp(recordName,signaList,N,N0,rawUnits,highResolution)

举例:在Matlab中输入以下内容通过rdsamp函数可以读出MITDB信号库中的一个recordName为‘100‘的信号并借助plot函数显示:

-[tm, signal]=rdsamp('mitdb/100',[],1000);

-plot(tm,signal(:,1));

其中,该信号中有两列采集信号,每列长为650000个采样点,所以获得的signal应该是一个2*650000的矩阵。但是可以在函数输入可选项:

SignalList:选择其中某列信号读入,输入为空默认为全部列,此处就是默认全部;

N:需要读入信号的长度,本身信号长为650000,此处只截取了其中1000长度;

N0:截取信号开始的采样点,此处默认为1,也就是截取了信号的1--1000采样点;

rawUnits:信号读入信号值需要的精度,具体请看帮助;

highResolution:信号分辨率,1为高,0为低,默认为0

最终结果为:

Figure1.MIT生理库记录名为‘100’信号的第一列前1000采样信号

其实这个示例程序就是rdsamp函数帮助文件中的一个示例程序,学会用Matlab的help函数,对于学习Matlab来说可以起到事半功倍的效果。一来我们可以看到我们要查看函数的所有信息和示例用法,直接copy粘贴在Matlab中就可以看到效果;而来我们可以用来检验该函数是否已经添加在PATH中。

2.读入注释信息文件(rdann)

[ann,type,subtype,chan,num,comments]=rdann(recordName,annotator,C,N,N0,type)

举一个rdann读入注释信息函数的例子:

-[tm, signal]=rdsamp('challenge/2013/set-a/a01');

-[ann]=rdann('challenge/2013/set-a/a01','fqrs');

-plot(tm,signal(:,1));hold on;grid on

-plot(tm(ann),signal(ann,1),'ro','MarkerSize',4)

本程序首先用rdsamp函数读入‘challenge’信号库中的命名为‘a01’的全部信号,再用‘rdann’函数读入与之对应的注释信息,类型为‘fqrs’,也就是QRS出现的位置,最后由plot显示所读入的信号和注释信息,如下图:

Figure2.'Challenge'生理库中命名为‘a01’信号及其QRS注释信息

图2中我们发现ECG信号显得很密集,因为我们没有选着读入信号的某些采样而是全部60000个采样点。其中红色的圆点就是注释文件中标注QRS出现的地方。其实注释信息我们可以通过一些算法求出来的,但是既然生理库中已经给出我们就可以直接用,而且他们官方的肯定经过矫正很准确的,既方便也有保证,直接拿来用何乐不为?另外如果我们想自己开发检测QRS或者R,可以将自己得到的结果和官方的注释信息对比,看看自己算法的准确率。

这一篇虽然只介绍了WFDB Toolbox for Matlab工具箱中两个简单的读入信号的函数,但是学会了使用和查看帮助文档的方法其他的函数也就很简单了。因为我们不是所有的函数都要用到,所以当遇到需要的函数在学习就好了,简单高效!

读入文件是第一步,对生理信号进行处理和有用信息提取才是重点,以后有机会我会介绍一些常用的方法以及在自己项目中用到的方法。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,142评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,298评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,068评论 0 351
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,081评论 1 291
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,099评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,071评论 1 295
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,990评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,832评论 0 273
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,274评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,488评论 2 331
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,649评论 1 347
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,378评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,979评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,625评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,796评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,643评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,545评论 2 352

推荐阅读更多精彩内容