代码分享|时频分析时绘制热图进行平滑的代码

大家好，我是茗创科技的周翊，近日分析数据的时候碰到一批采样率比较低，虽然对于分析的频段已经足够了，但是画出来的时频图却不好看，如下左图，本着对客户负责的原则，就想在现有的数据基础上给客户进行平滑作图，就有了右图的效果，是不是瞬间好看了很多？

闲话不多说，直接上代码

% 代码由茗创科技工程师 周翊编写 并无偿分享使用 转载注明来源
% 更多需要可加微信了解
% 茗创科技专注于脑影像数据处理，涵盖（fMRI,结构像,DTI,ASL,EEG/ERP,FNIRS,眼动）等，也欢迎
% 了解茗创科技的课程及数据处理服务，可添加微信号 17373158786 进行咨询。
% 
% Written By MCKJ tech 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.

%% 画热图平滑的代码--左图太粗糙，对其进行平滑作图呈现在论文里，建议计算机配置比较高，图片大于几百兆，比较耗内存
% img_P_BC是做时频分析时生成的时频矩阵
[X,Y] = meshgrid(1:size(img_P_BC,2),1:size(img_P_BC,1));% 时间*频率
[Xn,Yn] = meshgrid(1:0.1:size(img_P_BC,2),1:0.1:size(img_P_BC,1));%0.1足够了，0.01生成的数据3G大，会死机
img_P_BC_map1 = interp2(X,Y,img_P_BC,Xn,Yn,'linear');
imagesc(t,f,img_P_BC_map1);axis xy;

% Copyright (c) 2020-2021 MCKJ. All rights reserved.

实现环境为matlab2013+eeglab13

更多需要扫码咨询

↓ ↓ ↓