MIT-BIH数据库PSG信号提取（含EEG）

MIT-BIH Polysomnographic

找了一圈都是针对MIT-BIH ECG数据库的，写一下另一个数据库（含脑电的PSG）

简介：

数据库由 18 个记录组成，

导联数有四个、六个和七个三种，信号采样率是 250Hz，A/D 分辨率是 12Bit，

每个记录由四个文件组成，分为头文件、数据文件和注释文件。

其中注释文件包括两个，扩展名为.ecg 的注释文件和扩展名为.st 的注释文件。

其中扩展名为.ecg 的文件主要是对心电和呼吸信号做注释，

而扩展名为.st 的文件主要对脑电导联做注释。

 

附（18记录对应导联，其他为（四））：

C4-A1：slp01a、slp01b、slp32（七）、slp37（七）、slp41（七）

O2-A1：slp02a、slp02b

C3-O1：slp03、slp04、slp14、slp16、slp45（七）、slp48（七）、slp59（七）、slp60（七）、slp61（六）、slp66（七）、slp67x（七）

 

头文件读取

用txt打开头文件 ↓ 

slp01a 4 250/0.033333333(94) 1800000 23:07:00 19/1/1989    
      %4通道、采样率250Hz，180万个采样点
slp01a.dat 212 -200/mV 12 0 -17 59911 0 ECG
slp01a.dat 212 4.77778(-477)/mmHg 12 0 -248 19332 0 BP
slp01a.dat 【1】212 【2】-6430/mV 【3】12 【4】0 【5】252 【6】49594 【7】0 【8】EEG (C4-A1)  
          【1】即每行用三个字节表示两个数m1、m2，每个数 12 bits，故又称为 212 格式
          【2】信号的增益，ADC units/mv
          【3】ADC 的分辨率为 12 位
          【4】ADC 零值为 0
          【5】采样点的值为252
          【6】校验数    【7】输入输出可以以任何尺寸的块来执行     【8】导联
slp01a.dat 212 690/l 12 0 -180 912 0 Resp (sum)
# 44 M 89 32-01-89

数据文件的读取

signald = 'D:\mit-bih\data\slp01a.dat'; 
fid2=fopen(signald,'r');
SAMPLES2READ=2500;              %2500行，采样率250Hz,即读取5s数据
A= fread(fid2, [3, SAMPLES2READ], 'uint8')'; %212格式 3字节；'uint8'用1个字节表示 
fclose(fid2);

按212格式编排的数据

.dat文件的数据格式读取为每行三个字节,即三个八位的二进制数字,其内容含义为

      0000 0000  ||             0000 0000              ||  0000 0000

sign1(L)低八位信息||左四位sign2(R)信息,右四位sign1(L)信息||sign2(R)低八位信息

将第二字节的信息处理后,后四位移至第一字节最左位即得到完整的sign1

将第二字节的信息处理后,前四位移至第一字节最左位即得到完整的sign2.

M2H= bitshift(A(:,2), -4);        % 取A得第二列，-4表示字节向右移四位，即取信号2的高4位
M1H= bitand(A(:,2), 15);          % bitand返回A和B的按位‘与’，即取信号1的高4位（15表示‘0000 1111’）
PRL=bitshift(bitand(A(:,2),8),9);              % 取出字节低四位中最高位，向左移九位（8表示‘0000 1000’）
PRR=bitshift(bitand(A(:,2),128),5);            % 取出字节高四位中最高位，向左移五位（128表示‘1000 0000’）
M( : , 1)= bitshift(M1H,8)+ A(:,1)-PRL;        % 将M1H、M2H分别左移8位，即乘以2^8，再分别加上A(:,1)，A(:,2)，
M( : , 2)= bitshift(M2H,8)+ A(:,3)-PRR;        % 由于左移时把符号位也移动了，要减去符号位的值

按212格式解析的M

 

当信号四通道时：

前两通道212三字节，后两通道212三字节，经上述解析得 ↓

ECG	BP
EEG	Resp
ECG	BP
EEG	Resp

还要将读取的ECG、BP、EEG、Resp数据分别从M提取并放置于data的1234行

for  i=1:SAMPLES2READ
     if mod(i,2)==1
         data(1,(i+1)/2) = M(i,1);     % row1：ECG
     else
         data(3,i/2) = M(i,1);         % row3：EEG
     end  
end
for  i=1:SAMPLES2READ
     if mod(i,2)==1
         data(2,(i+1)/2) = M(i,2);     % row2：BP
     else
         data(4,i/2) = M(i,2);         % row4：Resp
     end  
end

 

plot EEG、ECG、Resp

 

当信号七通道时：

M1 = [];
for i=1:SAMPLES2READ
    M1 = [M1;M(i,1)];
	M1 = [M1;M(i,2)];
end
for  i=1:2*SAMPLES2READ
    switch mod(i,7)
        case 1
            data(1,(ceil(i/7))) = M1(i);     % row1：ECG
        case 2
            data(2,(ceil(i/7))) = M1(i);     % row2：BP
        case 3
            data(3,(ceil(i/7))) = M1(i);     % row3：EEG
        case 4
            data(4,(ceil(i/7))) = M1(i);     % row4：Resp (nasal)
        case 5
            data(5,(ceil(i/7))) = M1(i);     % row5：Resp (chest)
        case 6
            data(6,(ceil(i/7))) = M1(i);     % row6：EOG
        otherwise
            data(7,(ceil(i/7))) = M1(i);     % row7：EMG
     end  
end