本文主要是介绍代码分享|时频分析时绘制热图进行平滑的代码,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
大家好,我是茗创科技的周翊,近日分析数据的时候碰到一批采样率比较低,虽然对于分析的频段已经足够了,但是画出来的时频图却不好看,如下左图。本着对客户负责的原则,就想在现有的数据基础上给客户进行平滑作图,就有了右图的效果,是不是瞬间好看了很多?
% 代码由茗创科技工程师 周翊编写 并无偿分享使用 转载注明来源
% 更多需要可加微信了解
% 茗创科技专注于脑影像数据处理,涵盖(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
更多需要请联系微信号 17373158786 进行咨询
这篇关于代码分享|时频分析时绘制热图进行平滑的代码的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!