Sleep-EDF数据库实战：如何用Matlab快速处理睡眠分期标签（附完整代码）

Sleep-EDF数据库实战如何用Matlab快速处理睡眠分期标签附完整代码睡眠研究是神经科学和临床医学的重要领域而Sleep-EDF数据库作为公开可用的标准数据集为科研人员提供了宝贵的多导睡眠图(PSG)记录。但在实际应用中原始Hypnogram标签的格式往往需要经过转换才能用于分析。本文将分享一套完整的Matlab解决方案帮助您高效处理这些数据。1. 准备工作与环境配置在开始处理Sleep-EDF数据前我们需要做好以下准备工作数据获取从PhysioNet官网下载Sleep-EDF Database Expanded数据集建议选择v1.0.0或更高版本Matlab环境确保安装Signal Processing Toolbox和Statistics and Machine Learning Toolbox文件结构建议按以下方式组织项目文件夹/Sleep-EDF-Analysis ├── /data │ ├── /raw (存放原始EDF文件) │ └── /hypnograms (存放原始标签文件) ├── /processed │ ├── /labels (存放处理后的标签) │ └── /figures (存放分析图表) └── /scripts (存放Matlab代码)提示处理大规模数据时建议使用Matlab的并行计算功能可以显著提升处理速度。2. 理解Hypnogram标签格式Sleep-EDF数据库中的Hypnogram文件采用特定格式记录睡眠分期信息。原始文件通常包含三列数据列序号内容描述数据类型示例1阶段开始时间(秒)数值02阶段持续时间(秒)数值303睡眠阶段标识字符串Sleep stage W睡眠阶段标识与标准分期的对应关系如下% 睡眠阶段映射表 stage_mapping { Sleep stage W, 0; % 清醒期 Sleep stage 1, 1; % N1期 Sleep stage 2, 2; % N2期 Sleep stage 3, 3; % N3期 Sleep stage 4, 4; % N4期(已合并到N3) Sleep stage R, 5; % REM期 Sleep stage ?, -1; % 不确定阶段 };3. 标签转换的核心算法我们需要将原始标签转换为每30秒一个分期的标准格式。以下是核心处理函数function [stage_vector, time_vector] convert_hypnogram(raw_data) % 输入: raw_data - 从Hypnogram文件导入的原始数据 % 输出: stage_vector - 转换后的睡眠阶段向量 % time_vector - 对应的时间向量(秒) % 初始化参数 epoch_duration 30; % 标准分期时长(秒) % 提取原始数据中的时间和阶段信息 start_times raw_data(1:3:end); durations raw_data(2:3:end); stages raw_data(3:3:end); % 创建空向量存储结果 total_duration sum(durations); num_epochs ceil(total_duration / epoch_duration); stage_vector zeros(1, num_epochs); time_vector (0:num_epochs-1) * epoch_duration; % 填充阶段信息 current_pos 1; for i 1:length(stages) epoch_count round(durations(i) / epoch_duration); stage_vector(current_pos:current_posepoch_count-1) stages(i); current_pos current_pos epoch_count; end end注意实际应用中需要考虑边界情况处理如最后一个分期可能不足30秒。4. 批量处理与自动化流程对于包含多个受试者的大型研究手动处理每个文件效率低下。我们可以编写批量处理脚本% 批量处理Hypnogram文件脚本 data_dir data/hypnograms; output_dir processed/labels; % 获取所有Hypnogram文件 file_list dir(fullfile(data_dir, *.txt)); % 并行处理所有文件 parfor i 1:length(file_list) % 读取原始数据 raw_data importdata(fullfile(data_dir, file_list(i).name)); % 转换标签 [stage_vector, time_vector] convert_hypnogram(raw_data); % 保存结果 [~, basename] fileparts(file_list(i).name); save(fullfile(output_dir, [basename _label.mat]), ... stage_vector, time_vector); fprintf(已处理文件: %s\n, file_list(i).name); end为提高代码复用性我们可以将常用功能封装为独立函数function plot_sleep_stages(stage_vector, varargin) % 绘制睡眠阶段图 % 可选参数: % SubjectID - 受试者ID % SavePath - 保存路径 % 解析输入参数 p inputParser; addParameter(p, SubjectID, Unknown, ischar); addParameter(p, SavePath, , ischar); parse(p, varargin{:}); % 创建阶段名称映射 stage_names {W, N1, N2, N3, REM}; % 创建图形 figure(Color, white); imagesc(stage_vector); colormap(jet(5)); colorbar(Ticks, 0:4, TickLabels, stage_names); title(sprintf(睡眠阶段图 - 受试者 %s, p.Results.SubjectID)); xlabel(时间(30秒间隔)); ylabel(睡眠阶段); % 保存图形 if ~isempty(p.Results.SavePath) saveas(gcf, p.Results.SavePath); close(gcf); end end5. 常见问题与解决方案在实际应用中可能会遇到以下典型问题标签不一致问题现象不同版本的Sleep-EDF数据库可能使用略有不同的标签命名解决方案在预处理阶段添加标准化步骤function normalized_label normalize_stage_label(raw_label) % 统一处理不同版本的标签差异 raw_label strtrim(raw_label); if startsWith(raw_label, Sleep stage ) normalized_label raw_label; else % 处理其他格式的标签 normalized_label [Sleep stage raw_label]; end end时间对齐问题现象Hypnogram标签与PSG记录时间不完全同步解决方案添加时间校验步骤function is_valid check_time_alignment(hypnogram, psg_duration) % 检查Hypnogram总时长是否与PSG记录匹配 total_hypno_time sum(hypnogram(2:3:end)); time_diff abs(total_hypno_time - psg_duration); is_valid time_diff 30; % 允许30秒以内的差异 end性能优化技巧对于大型数据集可以使用内存映射文件技术% 使用memmapfile处理大文件 m memmapfile(large_hypnogram.bin, ... Format, {double, [3 inf], data}); chunk_size 1000; for i 1:chunk_size:size(m.Data.data, 2) chunk m.Data.data(:, i:min(ichunk_size-1, end)); % 处理数据块 end6. 进阶应用与扩展掌握了基础标签处理后可以进一步实现以下高级功能睡眠质量指标计算function metrics calculate_sleep_metrics(stage_vector) % 计算各类睡眠质量指标 total_sleep_time sum(stage_vector 0) * 30 / 60; % 分钟 rem_percentage sum(stage_vector 5) / length(stage_vector) * 100; wake_after_sleep_onset sum(stage_vector 0 cumsum(stage_vector 0) 0); metrics struct(... TST, total_sleep_time, ... REM, rem_percentage, ... WASO, wake_after_sleep_onset); end多受试者统计分析% 示例分析多受试者的REM比例 subject_files dir(processed/labels/*.mat); rem_percentages zeros(1, length(subject_files)); for i 1:length(subject_files) data load(fullfile(subject_files(i).folder, subject_files(i).name)); metrics calculate_sleep_metrics(data.stage_vector); rem_percentages(i) metrics.REM; end % 绘制分布图 figure; histogram(rem_percentages, BinWidth, 2); xlabel(REM睡眠比例(%)); ylabel(受试者数量); title(REM睡眠比例分布);与EEG数据同步分析function aligned_data align_with_eeg(sleep_stages, eeg_data, fs) % 将睡眠阶段与EEG数据对齐 % fs - EEG采样频率(Hz) epoch_samples 30 * fs; % 每个分期对应的样本数 num_epochs length(sleep_stages); aligned_data zeros(1, num_epochs * epoch_samples); for i 1:num_epochs start_sample (i-1)*epoch_samples 1; end_sample i*epoch_samples; aligned_data(start_sample:end_sample) sleep_stages(i); end end在实际项目中我发现使用面向对象的方式封装这些功能可以大大提高代码的可维护性。例如创建一个SleepStudy类将数据加载、处理和可视化方法整合在一起。这种方法特别适合长期研究项目可以方便地扩展新功能而不破坏现有代码结构。