脑波和眼动连通性分析 python-(1)

本文介绍如何使用Python处理bdf文件,从元数据中提取信息,并将脑电信号按delta、theta、alpha、beta波段进行通道划分。
摘要由CSDN通过智能技术生成

bdf文件中提取元数据并分成18个channel,特征波分别为delta,theta,alpha,beta波

import pandas as pd
import numpy as np
import mne
from mne.io import concatenate_raws, read_raw_edf
import matplotlib.pyplot as plt
import datetime
from eeg_processing import iir,processing
import os

def mkdir(path_3):
    folder = os.path.exists(path_3)
    if not folder:  # 判断是否存在文件夹如果不存在则创建为文件夹
        os.makedirs(path_3)  # makedirs 创建文件时如果路径不存在会创建这个路径
        print("---  new folder...  ---")
        print("---  OK  ---")
    else:
        print("---  There is this folder!  ---")
def read_file(eeg_path):
    raw = mne.io.read_raw_edf(eeg_path, preload=True)
    raw.pick_types(meg=False, eeg=True, eog=False)
    eeg, times = raw[0:18]
    fs = raw.info['sfreq']
    ch_names = raw.info['ch_names'][0:18]
    return eeg, times, fs, ch_names

def csv_read():
    raw = read_raw_edf(eeg_path, preload=False)
    print(raw)
    data = pd.DataFrame(pd.read_csv(eeg_path, encoding="utf-8", engine='python'))
    Fp2= np.array(data['EEG Fp2-Cz'])  # data length : 60fps * time (video time)
    print(len(Fp2))
def start_time():
    eeg_time_filename = mid_path + '/eeg/read me.txt'
    for line in open(eeg_time_filename, encoding='utf-8').readlines():
        url = line.rstrip()  # 去掉末尾隐藏的\n,避免出现空行。
        data = url.split('\t')[0]
        startime = data.split('time:')[-1]
        start_eeg = datetime.datetime.strptime(startime, "%Y-%m-%d %H:%M:%S.%f")
        print('startime:',startime)
        print('start:',start_eeg)
        return start_eeg
def 
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值