不成功的TCA代码
2024-10-13 20:57:53
%--brain mask with the brain tissue
mask_name = 'C:\Users\Administrator\Desktop\workspace\preprocessed\masks\within_brain_mask.nii';
M = load_untouch_nii( mask_name ); % load mask NIFTI
mask = double(M.img>0); % get 3d v %--brain functional 4d data
data_4d = 'C:\Users\Administrator\Desktop\workspace\preprocessed\4d\func_3d.nii';
% data_4d = 'C:\Users\Administrator\Desktop\phycaa_workspace\phycaa_plus_2104_03_27\_PHYCAA_step1+2.nii'; V = load_untouch_nii( data_4d ); %--transform 4d array to 2d array, using brain_mask
within_brain_voxels = nifti_to_mat(V,M); nt_matrix = within_brain_voxels; n = size(nt_matrix,1);
t = size(nt_matrix,2); mean_value = mean ( mean(nt_matrix,2) ); tmp_matrix = nt_matrix - mean_value; mean_value2 = mean(tmp_matrix,2) ;
% normalise_data = abs(( nt_matrix - repmat( mean_value, 1 , t )) ./ repmat( mean_value, 1 , t )) ; normalise_data = abs(( tmp_matrix - repmat( mean_value2, 1 , t )) ./ repmat( mean_value2, 1 , t )) ; % normalise_data = abs(( nt_matrix - mean_value) )/ mean_value ; w_data = normalise_data; for i=1:n tmp_rows = normalise_data(i,:);
max_value = max(tmp_rows); tmp_rows_index = (tmp_rows == max_value);
tmp_rows1 = tmp_rows.* double(tmp_rows_index); % w_data(i,:) = double(tmp_rows_index);
w_data(i,:) = tmp_rows1; end reference_waveform = sum(w_data); plot(reference_waveform);
最新文章
- Android入门(十):界面的布局方式及其实际应用
- SQL查询——同一张表的横向与纵向同时比较
- Big Number--hdu1018(数学)
- Asp.Net Web API 2第五课——Web API路由
- linux svn安装和配置
- 【IOS学习基础】NSObject.h学习
- 微信小程序 PHP后端form表单提交实例详解
- Java数据结构和算法 - 递归
- MIUI9系统怎么启用Root超级权限的经验
- 详解javascript事件绑定使用方法
- iOS--------获取当前连接的WiFi以及IP地址
- bootstrap 三个文件的引入
- WebSocket原理与实践
- C#USB设备枚举Kernel32的PInvoke
- js sort
- js-权威指南学习笔记14
- django restframework 简单总结
- RTC终于tm的通了
- RESTORE DATABASE命令还原SQLServer 2005 数据库
- OFBiz实战(1):整合Groovy+FreeMaker