代码分享|静息态频域各指标的计算,建议收藏

2024-06-14 17:38

本文主要是介绍代码分享|静息态频域各指标的计算,建议收藏,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

       大家好,我是茗创科技的工程师周翊,静息态数据在做频域分析的时候常用的指标有四个

实现代码如下:


% 代码由茗创科技工程师 周翊编写 并无偿分享使用 转载注明来源
% 更多需要可加微信了解
% 茗创科技专注于脑影像数据处理,涵盖(fMRI,结构像,DTI,ASL,EEG/ERP,FNIRS,眼动)等,也欢迎
% 了解茗创科技的课程及数据处理服务,可添加微信号 17373158786 进行咨询。
% Written By Yizhou
% Using the code without proper understanding the code and relevant background
% of EEG may lead to confusion, incorrect data analyses,or misinterpretations
% of results.
% The author assumes NO responsibility for inappropriate or incorrect use
% of this code.
% WX:17373158786%频域分析设置
% Fs采样率
Fs = 512;
%分段长度 采样率*分段长度
L = 1024;
% 
% 把分段长度变成值最拼接的2的n次冥值,目的是提高运算速度
NFFT = 2^nextpow2(L);
%创建频域列表
f = Fs/2 * linspace(0,1,NFFT/2+1);% 路径设置为预处理的文件夹
anal_path = 'E:\Docu\Work\1_Projects\Rest_8\anal_V2\resting';cd(anal_path);
files = dir([anal_path,filesep,'*.set']);% 新建一个文件夹用于存在静息态指标文件
if ~exist([anal_path,filesep,'fft_result'],'dir')mkdir([anal_path,filesep,'fft_result']);
end
fre_path = [anal_path,filesep,'fft_result'];%%% 频域结果储存路径
cd(fre_path);for sub_No = 1:length(files)setname = files(sub_No).name;EEG = pop_loadset([anal_path, filesep, setname]);%对于所有通道for ii=1:size(EEG.data,1)%对于所有分段for jj=1:size(EEG.data,3)%提取第ii个通道 第jj个分段 所有时间点数据y = squeeze(EEG.data(ii,:,jj));%做FFT变换 求真实幅值----------------------------------指标1:幅值  这个指标使用较少temp = fft(y,NFFT)/ L;fuzhi(jj,:) = 2*abs(temp(1:NFFT/2+1));% 做FFT变换,求功率(uV平方)%-------------------------指标2:功率 使用较多temp = fft(y,NFFT)/L;power(jj,:)=(2*abs(temp(1:NFFT/2+1))).^2; % fft results, in power% 做FFT变换,求功率谱密度(PSD)(单位 uV平方/ Hz)%----指标3:功率谱密度 使用较多temp = 2*abs(fft(y,NFFT)).^2/L/Fs;PSD(jj,:) = temp(1:NFFT/2 + 1);% 做FFT变换,求功率谱密度(PSD)(单位 dB )%----------指标4:% 对指标3做log变换,目的是压缩低频,延展高频最推荐的指标temp = 2*abs(fft(y,NFFT)).^2/L/Fs;DB_PSD(jj,:) = 10*log10(temp(1:NFFT/2+1));end%数据汇总:被试*通道*频率 mean_fuzhi(sub_No ,ii,:) = squeeze(mean(fuzhi,1)); clear fuzhi;mean_power(sub_No ,ii,:) = squeeze(mean(power,1)); clear power;mean_PSD(sub_No ,ii,:) = squeeze(mean(PSD,1)); clear PSD;mean_DB_PSD(sub_No ,ii,:) = squeeze(mean(DB_PSD,1)); clear DB_PSD;end% 保存各个指标end% 保存结果文件到结果文件的路径,四个指标分别保存
save ('fuzhi.mat', 'mean_fuzhi');
save ('power.mat', 'mean_power');
save ('PSD.mat', 'mean_PSD');
save ('DB_PSD.mat', 'mean_DB_PSD');
% %             clear mean*;% Copyright (c) 2020-2021 MCKJ. All rights reserved.

更多需要请联系微信号 17373158786 进行咨询

这篇关于代码分享|静息态频域各指标的计算,建议收藏的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1061075

相关文章

Python并行处理实战之如何使用ProcessPoolExecutor加速计算

《Python并行处理实战之如何使用ProcessPoolExecutor加速计算》Python提供了多种并行处理的方式,其中concurrent.futures模块的ProcessPoolExecu... 目录简介完整代码示例代码解释1. 导入必要的模块2. 定义处理函数3. 主函数4. 生成数字列表5.

Python虚拟环境与Conda使用指南分享

《Python虚拟环境与Conda使用指南分享》:本文主要介绍Python虚拟环境与Conda使用指南,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、python 虚拟环境概述1.1 什么是虚拟环境1.2 为什么需要虚拟环境二、Python 内置的虚拟环境工具

Python实例题之pygame开发打飞机游戏实例代码

《Python实例题之pygame开发打飞机游戏实例代码》对于python的学习者,能够写出一个飞机大战的程序代码,是不是感觉到非常的开心,:本文主要介绍Python实例题之pygame开发打飞机... 目录题目pygame-aircraft-game使用 Pygame 开发的打飞机游戏脚本代码解释初始化部

Java中Map.Entry()含义及方法使用代码

《Java中Map.Entry()含义及方法使用代码》:本文主要介绍Java中Map.Entry()含义及方法使用的相关资料,Map.Entry是Java中Map的静态内部接口,用于表示键值对,其... 目录前言 Map.Entry作用核心方法常见使用场景1. 遍历 Map 的所有键值对2. 直接修改 Ma

深入解析 Java Future 类及代码示例

《深入解析JavaFuture类及代码示例》JavaFuture是java.util.concurrent包中用于表示异步计算结果的核心接口,下面给大家介绍JavaFuture类及实例代码,感兴... 目录一、Future 类概述二、核心工作机制代码示例执行流程2. 状态机模型3. 核心方法解析行为总结:三

python获取cmd环境变量值的实现代码

《python获取cmd环境变量值的实现代码》:本文主要介绍在Python中获取命令行(cmd)环境变量的值,可以使用标准库中的os模块,需要的朋友可以参考下... 前言全局说明在执行py过程中,总要使用到系统环境变量一、说明1.1 环境:Windows 11 家庭版 24H2 26100.4061

pandas实现数据concat拼接的示例代码

《pandas实现数据concat拼接的示例代码》pandas.concat用于合并DataFrame或Series,本文主要介绍了pandas实现数据concat拼接的示例代码,具有一定的参考价值,... 目录语法示例:使用pandas.concat合并数据默认的concat:参数axis=0,join=

C#代码实现解析WTGPS和BD数据

《C#代码实现解析WTGPS和BD数据》在现代的导航与定位应用中,准确解析GPS和北斗(BD)等卫星定位数据至关重要,本文将使用C#语言实现解析WTGPS和BD数据,需要的可以了解下... 目录一、代码结构概览1. 核心解析方法2. 位置信息解析3. 经纬度转换方法4. 日期和时间戳解析5. 辅助方法二、L

Python使用Code2flow将代码转化为流程图的操作教程

《Python使用Code2flow将代码转化为流程图的操作教程》Code2flow是一款开源工具,能够将代码自动转换为流程图,该工具对于代码审查、调试和理解大型代码库非常有用,在这篇博客中,我们将深... 目录引言1nVflRA、为什么选择 Code2flow?2、安装 Code2flow3、基本功能演示

IIS 7.0 及更高版本中的 FTP 状态代码

《IIS7.0及更高版本中的FTP状态代码》本文介绍IIS7.0中的FTP状态代码,方便大家在使用iis中发现ftp的问题... 简介尝试使用 FTP 访问运行 Internet Information Services (IIS) 7.0 或更高版本的服务器上的内容时,IIS 将返回指示响应状态的数字代