医学信息技术导论2.docx
- 文档编号:2990811
- 上传时间:2022-11-16
- 格式:DOCX
- 页数:11
- 大小:152.48KB
医学信息技术导论2.docx
《医学信息技术导论2.docx》由会员分享,可在线阅读,更多相关《医学信息技术导论2.docx(11页珍藏版)》请在冰豆网上搜索。
医学信息技术导论2
课程作业
课程名称:
医学信息技术导论(matlab部分)
姓名:
梅世宇
学号:
201058002
班级:
电医1001
教师:
唐洪老师
日期:
2012年7月12日
一.题目
搜索MIT-BIH数据库的网络地址
(1)下载一例心律失常病人的心电信号
(2)下载一例正常人的心电信号
(3)利用数据库提供的matlab函数,读取心电信号
(4)查阅资料,学习心率变异的概念及应用
(5)编写matlab程序对心电信号作图,找出所有R波发生位置(数据中提供了R波的位置),计算出心动周期,画出心率随时间变化的曲线,体会正常心律和异常心律的差异。
二.原理
心电图产生
心脏不断进行有节奏的收缩和舒张活动,血液才能在循环系统中不断地流动。
心脏同生物体地其他组织一样在兴奋过程中产生电位变化。
心肌细胞的生物电现象与神经肌肉的生物电现象一样分为静息电位和兴奋时的动作电位。
人与哺乳动物的心室肌细胞的静息电位为-90mv。
当心肌细胞由静止状态进入兴奋时,即产生动作电位。
心肌的动作电位起初迅速上升,跟着是一相当长的周期,在此期间电势最初下降很慢,然后较快,最后恢复到静息电位水平。
心脏传导系统主要由窦房结、结间束、房室结、房室束、左右束支以及浦肯野氏纤维构成。
心肌动作电位开始于右心房的窦房结,然后沿结间束传至房室结,再传至房室束,沿房室束的左右束支,几乎同时传至两侧心室的全部心内膜后,经浦肯野氏纤维传至心室外膜。
由于人体体液是良导体,心脏的每个心动周期所伴随的生物电变化可以传达到身体表面的各个部位。
各部分心肌细胞产生的动作电位综合而成的电信号由引导电极从体表或胸腔测得,经放大后显示或描记下来的波形,就是心电图(Electrocardiogram,ECG),它反映了心脏兴奋的产生、传导和恢复过程中的生物电变化。
一般正常人心率为60-100次/分,心电图的频率范围在0.05Hz~100Hz,幅值为10μV~5mV。
心电波形简介
伴随着心脏收缩与舒张,心电波形依次会出现P波、PR间隔、QRS复合波、ST段、T波、U波的周期变化。
为了更准确地诊断心血管疾病,最常采用的是标准12导联系统测量。
典型的心电图由下列各波组成:
P波:
由心房的激动所产生,前一半主要由右心房所产生,后一般主要由左
心房所产生;反映左右心房的除极过程。
正常P波的宽度不超过0.11s,最高幅
值不超过2.5mm。
PR间隔:
从P波起点到QRS波群起点的时间间隔。
表示从心房去极化到心
室开始去极化的电活动过程。
正常范围是0.12-0.2s,随着年龄增长而有加长的趋势。
QRS波群:
从Q波开始至S波终了的时间间期。
代表左右心室除极电位变化。
正常人最高时限不超过0.1s。
ST段:
从QRS波群的终点到T波起点的一段。
正常人ST段是接近基线的,与基线间的距离不超过0.05mm。
T波:
代表心室除极后复极时所产生的电位影响。
T波在ST段之后发生,所占时间较长;波形比较低,在R波为主的心电图上,T波不应低于R波的1/10。
心电研究的应用
心脏的活动状态除了反映在心律上,更主要是通过ECG中P、QRS、和T形波以及它们的各种ST段的参数、RR间期、PP间期、PR间期、QT间期(幅度、宽度、斜率等)反映。
ECG信号分析,重要的是进行ECG信号的特征点参数提取和特征波形段的识别和定位,这是ECG分析诊断系统的关键,其准确性、可靠性的好坏决定着诊断与治疗心脏病患者的效果,乃至挽救病人生命的成败。
利用计算机辅助医生进行ECG分析,能够帮助医生准确高效诊断,提高诊断的可靠性。
自动检测为心电自动诊断系统提供准确的各项参数数据,心电检测主要集中在QRS波检测,P波检测,T波检测。
QRS波检测是ECG检测中的首要问题,QRS波检测不仅是心律失常的重要依据,而且只有在QRS波确定以后才能分析ECG的其它细节信息。
精确定位R波再进一步用算法判定QRS复合波群的起始点和终点。
常用的R峰检测方法有幅度法、斜率法、面积法、阈值法、小波变换法。
两种或多种的方法结合的QRS检测方法是当今的趋势。
心电自动检测与分析的研究,也已经有了一些相对比较成熟的方法。
可以考虑先试用一些相对简单的算法,得到一定的检测结果,在此基础之上,再实现比较复杂效果更好的算法。
目前,心电检测在进一步提高信号的检测率、准确率以及使得分析诊断系统小型化更加智能化方面仍有大量的研究空间。
总的来QRS波检测的途径可以分为硬件检测和软件检测。
physionet心电数据库
PhysioBank包含供生物医学研究用的生理信号数据存档
http:
//www.physionet.org/
MIT-BIHArrhythmiaDatabase心率失常数据库.
MIT-BIHSTChangeDatabaseST变化数据库
MIT-BIHSupraventricularArrhythmiaDatabase前心室心率失常数据库
MIT-BIHLong-TermDatabase长时间心电数据库
EuropeanST-TDatabase欧洲ST-T数据库
MIT-BIHNormalSinusRhythmDatabase正常窦性心律数据库
PhysioToolkit包含生理信号分析和处理的开放源代码的软件库
http:
//www.physionet.org/physiotools/wfdb-windows-quick-start.shtml
数据库中心电数据格式、标注格式、头文件格式:
心电数据库包含的就是心电数据记录的集合。
每一条记录,以记录名称标识。
大多数情况下,一条记录包含至少三个文件,分别是以记录名称为文件名,后面接不同的后缀,不同的后缀表明文件的不同类型。
如对于MIT-BIH心律失常据库(MIT-BIHArrhythmiaDatabase)中的记录100,有三个文件100.dat,100.hea,100.atr均是该条记录的文件。
头文件(.hea)是短的文本文件,长度为100字节左右,包含记录的一般的信息,如导联信息、采样频率、记录长度、数据文件格式、病人年龄、性别、用药情况病情简述等等;数据文件(.dat)记录了A/D变换后心电图的数据,采用212格式存储,即3个字节存储2个采样点的数据;注释文件(.atr)记录了专家们人工对各个心电波形的标注结果,用于校验算法有效性和准确程度。
三.程序
Abnormal.m
PATH='C:
\Users\Administrator\Desktop\matlab\重新整合02\心率失常数据';%path,wheredataaresaved
HEADERFILE='800.hea';%header-fileintextformat
ATRFILE='800.atr';%attributes-fileinbinaryformat
DATAFILE='800.dat';%data-file
SAMPLES2READ=500;%numberofsamplestoberead
%incaseofmorethanonesignal:
%2*SAMPLES2READsamplesareread
%------LOADHEADERDATA--------------------------------------------------
fprintf(1,'\\n$>WORKINGON%s...\n',HEADERFILE);
signalh=fullfile(PATH,HEADERFILE);
fid1=fopen(signalh,'r');
z=fgetl(fid1);
A=sscanf(z,'%*s%d%d%d',[1,3]);
nosig=A
(1);%numberofsignals
sfreq=A
(2);%samplerateofdata
clearA;
fork=1:
nosig
z=fgetl(fid1);
A=sscanf(z,'%*s%d%d%d%d%d',[1,5]);
dformat(k)=A
(1);%format;hereonly212isallowed
gain(k)=A
(2);%numberofintegerspermV
bitres(k)=A(3);%bitresolution
zerovalue(k)=A(4);%integervalueofECGzeropoint
firstvalue(k)=A(5);%firstintegervalueofsignal(totestforerrors)
end;
fclose(fid1);
clearA;
%------LOADBINARYDATA--------------------------------------------------
ifdformat~=[212,212],error('thisscriptdoesnotapplybinaryformatsdifferentto212.');end;
signald=fullfile(PATH,DATAFILE);%datainformat212
fid2=fopen(signald,'r');
A=fread(fid2,[3,SAMPLES2READ],'uint8')';%matrixwith3rows,each8bitslong,=2*12bit
fclose(fid2);
M2H=bitshift(A(:
2),-4);
M1H=bitand(A(:
2),15);
PRL=bitshift(bitand(A(:
2),8),9);%sign-bit
PRR=bitshift(bitand(A(:
2),128),5);%sign-bit
M(:
1)=bitshift(M1H,8)+A(:
1)-PRL;
M(:
2)=bitshift(M2H,8)+A(:
3)-PRR;
ifM(1,:
)~=firstvalue,error('inconsistencyinthefirstbitvalues');end;
switchnosig
case2
M(:
1)=(M(:
1)-zerovalue
(1))/gain
(1);
M(:
2)=(M(:
2)-zerovalue
(2))/gain
(2);
TIME=(0:
(SAMPLES2READ-1))/sfreq;
case1
M(:
1)=(M(:
1)-zerovalue
(1));
M(:
2)=(M(:
2)-zerovalue
(1));
M=M';
M
(1)=[];
sM=size(M);
sM=sM
(2)+1;
M(sM)=0;
M=M';
M=M/gain
(1);
TIME=(0:
2*(SAMPLES2READ)-1)/sfreq;
otherwise%thiscasedidnotappearuptonow!
%hereMhastobesorted!
!
!
disp('Sortingalgorithmformoret
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 医学 信息技术 导论
![提示](https://static.bdocx.com/images/bang_tan.gif)