保存转换后的 .bdf 文件
Saving a transformed .bdf file
我使用了 MATLAB R2011b 和 EEGLab 13.4.3b to load a Biosemi .bdf 文件(EEG 记录),其中包含几个虚假触发器。触发器的位置存储在 header.BDF.Trigger.POS
中,触发器类型存储在 header.BDF.Trigger.TYP
.
中
加载 the file 后:
header = sopen('05-AM-200-Deci.bdf','r',[1:40],'OVERFLOWDETECTION:OFF');
我对 header.BDF.Trigger.POS
和 header.BDF.Trigger.TYP
向量进行了转换并替换了它们。
注意:这些转换应该不会有什么不同,因为即使我只加载 bdf 文件然后尝试将其保存回来,我也会得到完全相同的错误。
这是 header
的结构(用 get()
提取):
header =
FileName: '05-AM-200-Deci.bdf'
FILE: [1x1 struct]
TYPE: 'BDF'
ErrNum: [1025 0]
ErrMsg: ''
keycode: [1x34 double]
NS: 41
SampleRate: 1
T0: [2015 7 6 17 41 44]
Filter: [1x1 struct]
FLAG: [1x1 struct]
EVENT: [1x1 struct]
VERSION: -1
Patient: [1x1 struct]
PID: '05-AM-200'
RID: '05-AM-200'
HeadLen: 10752
reserved1: '24BIT '
NRec: 1207
Dur: 1
AS: [1x1 struct]
Label: {41x1 cell}
Transducer: {41x1 cell}
PhysDim: {41x1 cell}
PhysMin: [1x41 double]
PhysMax: [1x41 double]
DigMin: [1x41 double]
DigMax: [1x41 double]
PreFilt: [41x80 char]
GDFTYP: [1x41 double]
SPR: 1
THRESHOLD: [41x2 double]
Cal: [1x41 double]
Off: [1x41 double]
Calib: [41x40 double]
BDF: [1x1 struct]
PhysDimCode: [41x1 double]
ELEC: [1x1 struct]
LeadIdCode: [41x1 double]
CHANTYP: ' '
InChanSelect: [40x1 double]
我尝试使用 EEGLab 函数保存转换后的文件 (header
),但无济于事。我知道如何使用 GUI 在 EEGLab 中保存 bdf 文件。但是,使用 GUI 不允许进行我需要的转换。
为了保存文件,我尝试使用 pop_writeeeg()
and swrite()
,但其中 none 有效。
我使用了以下命令:
先拍:
`pop_writeeeg(header);`
Returns:
Reference to non-existent field 'trials'.
Error in pop_writeeeg (line 38)
if EEG.trials > 1
然后我试了:
`pop_writeeeg(header, '05-new', 'TYPE', 'BDF');`
返回:
Reference to non-existent field 'chanlocs'.
Error in pop_writeeeg (line 66)
if ~isempty(EEG.chanlocs)
我的下一步是按以下方式使用 swrite()
:
`swrite(header, '05-new');`
返回:
Error SWRITE can not be applied, File 05-AM-200-Deci.bdf is not opened in WRITE mode
我很确定我在这里遗漏了一些简单的东西。
有谁知道如何将转换后的 EEG (bdf) 数据保存回 bdf 文件?
EEGLAB 工具箱中没有您需要的功能;你想要 Biosig toolbox(它是 EEGLAB 的插件)。 Biosig 函数的范例是不寻常的。它使用结构中带有文件句柄的数据结构。因此,要打开一个文件然后保存修改后的数据,您需要执行以下步骤:
1) 打开BDF文件,读入数据。此步骤检索文件元数据并创建元数据的结构(如问题中所示)。
[raw_data,meta_data] = sload('05-AM-200-Deci.bdf',[1:40],'OVERFLOWDETECTION:OFF');
2) 即使转换仅涉及元数据,您也必须使用原始数据和元数据重新写入文件。您的数据必须位于大小为 Rows x Number_of_Channels 的变量中。根据 .bdf 文件的特性,您可能必须添加一个额外的列来表示状态通道。
% do the transforms here. If it does not affect the raw data, no need to change variables.
transformed_data = raw_data; % transformed data to be saved to a file
meta_data.fZ = zeros([1 meta_data.NS]);
transformed_data(:,meta_data.NS)=0;
3) 打开文件进行写入警告:这将删除磁盘上的文件
meta_data = sopen(meta_data,'w');
4) 将您的数据保存到文件中:
swrite(meta_data,transformed_data);
5) 关闭文件:
sclose(meta_data);
我使用可用的测试文件对此进行了测试 here。 BDF 格式很旧并且有很多选项,因此此脚本中有一个额外的步骤来填充 "Status" 通道(通道 17),以使其与 Biosig 功能一起使用。
[s,hdr] = sload('Newtest17-2048.bdf'); % get the raw data from the file
hdr.fZ = zeros([1 hdr.NS]); s(:,hdr.NS)=0; % (create "missing" channel values)
hdr = sopen(hdr,'w'); % open the file for writing (ERASES FILE!)
swrite(hdr,s); % write data to the file
sclose(hdr); % close the file
问题中链接的 .bdf 文件也需要 "missing channel" 的 hack。
我使用了 MATLAB R2011b 和 EEGLab 13.4.3b to load a Biosemi .bdf 文件(EEG 记录),其中包含几个虚假触发器。触发器的位置存储在 header.BDF.Trigger.POS
中,触发器类型存储在 header.BDF.Trigger.TYP
.
加载 the file 后:
header = sopen('05-AM-200-Deci.bdf','r',[1:40],'OVERFLOWDETECTION:OFF');
我对 header.BDF.Trigger.POS
和 header.BDF.Trigger.TYP
向量进行了转换并替换了它们。
注意:这些转换应该不会有什么不同,因为即使我只加载 bdf 文件然后尝试将其保存回来,我也会得到完全相同的错误。
这是 header
的结构(用 get()
提取):
header =
FileName: '05-AM-200-Deci.bdf'
FILE: [1x1 struct]
TYPE: 'BDF'
ErrNum: [1025 0]
ErrMsg: ''
keycode: [1x34 double]
NS: 41
SampleRate: 1
T0: [2015 7 6 17 41 44]
Filter: [1x1 struct]
FLAG: [1x1 struct]
EVENT: [1x1 struct]
VERSION: -1
Patient: [1x1 struct]
PID: '05-AM-200'
RID: '05-AM-200'
HeadLen: 10752
reserved1: '24BIT '
NRec: 1207
Dur: 1
AS: [1x1 struct]
Label: {41x1 cell}
Transducer: {41x1 cell}
PhysDim: {41x1 cell}
PhysMin: [1x41 double]
PhysMax: [1x41 double]
DigMin: [1x41 double]
DigMax: [1x41 double]
PreFilt: [41x80 char]
GDFTYP: [1x41 double]
SPR: 1
THRESHOLD: [41x2 double]
Cal: [1x41 double]
Off: [1x41 double]
Calib: [41x40 double]
BDF: [1x1 struct]
PhysDimCode: [41x1 double]
ELEC: [1x1 struct]
LeadIdCode: [41x1 double]
CHANTYP: ' '
InChanSelect: [40x1 double]
我尝试使用 EEGLab 函数保存转换后的文件 (header
),但无济于事。我知道如何使用 GUI 在 EEGLab 中保存 bdf 文件。但是,使用 GUI 不允许进行我需要的转换。
为了保存文件,我尝试使用 pop_writeeeg()
and swrite()
,但其中 none 有效。
我使用了以下命令:
先拍:
`pop_writeeeg(header);`
Returns:
Reference to non-existent field 'trials'. Error in pop_writeeeg (line 38) if EEG.trials > 1
然后我试了:
`pop_writeeeg(header, '05-new', 'TYPE', 'BDF');`
返回:
Reference to non-existent field 'chanlocs'. Error in pop_writeeeg (line 66) if ~isempty(EEG.chanlocs)
我的下一步是按以下方式使用
swrite()
:`swrite(header, '05-new');`
返回:
Error SWRITE can not be applied, File 05-AM-200-Deci.bdf is not opened in WRITE mode
我很确定我在这里遗漏了一些简单的东西。
有谁知道如何将转换后的 EEG (bdf) 数据保存回 bdf 文件?
EEGLAB 工具箱中没有您需要的功能;你想要 Biosig toolbox(它是 EEGLAB 的插件)。 Biosig 函数的范例是不寻常的。它使用结构中带有文件句柄的数据结构。因此,要打开一个文件然后保存修改后的数据,您需要执行以下步骤:
1) 打开BDF文件,读入数据。此步骤检索文件元数据并创建元数据的结构(如问题中所示)。
[raw_data,meta_data] = sload('05-AM-200-Deci.bdf',[1:40],'OVERFLOWDETECTION:OFF');
2) 即使转换仅涉及元数据,您也必须使用原始数据和元数据重新写入文件。您的数据必须位于大小为 Rows x Number_of_Channels 的变量中。根据 .bdf 文件的特性,您可能必须添加一个额外的列来表示状态通道。
% do the transforms here. If it does not affect the raw data, no need to change variables.
transformed_data = raw_data; % transformed data to be saved to a file
meta_data.fZ = zeros([1 meta_data.NS]);
transformed_data(:,meta_data.NS)=0;
3) 打开文件进行写入警告:这将删除磁盘上的文件
meta_data = sopen(meta_data,'w');
4) 将您的数据保存到文件中:
swrite(meta_data,transformed_data);
5) 关闭文件:
sclose(meta_data);
我使用可用的测试文件对此进行了测试 here。 BDF 格式很旧并且有很多选项,因此此脚本中有一个额外的步骤来填充 "Status" 通道(通道 17),以使其与 Biosig 功能一起使用。
[s,hdr] = sload('Newtest17-2048.bdf'); % get the raw data from the file
hdr.fZ = zeros([1 hdr.NS]); s(:,hdr.NS)=0; % (create "missing" channel values)
hdr = sopen(hdr,'w'); % open the file for writing (ERASES FILE!)
swrite(hdr,s); % write data to the file
sclose(hdr); % close the file
问题中链接的 .bdf 文件也需要 "missing channel" 的 hack。