所有g_sub*_component_ica_s1_.nii中提取几个component相加制作成mask,避免赋值需要标量结构体。matlab

%以下是一个MATLAB代码示例，可以从所有g_sub*component_ica_s1.nii文件中提取几个component并相加制作成mask，并保存为相应的mask文件。
%这里假设需要提取的component分别为1、2、4，并且所有g_sub*component_ica_s1.nii文件都存储在同一个文件夹中。

% 设置需要提取的component和NIFTI文件路径
components = [9, 10, 24];    % 需要提取的component编号
folder_path = 'E:/icac/gica40';    % NIFTI文件所在的文件夹路径

% 获取文件夹中所有g_sub*_component_ica_s1_.nii文件的文件名
file_names = dir(fullfile(folder_path, 'g_sub*_component_ica_s1_.nii'));
n_files = length(file_names);

% 读取其中一个NIFTI文件，用于后续创建共同的mask文件
nii_template = load_nii(fullfile(folder_path, file_names(1).name));

% 创建一个空的数据矩阵，用于存储所有NIFTI文件中的component数据
mask_data = zeros(size(nii_template.img));

% 循环处理每个文件
for i = 1:n_files
    % 读取NIFTI文件
    nii = load_nii(fullfile(folder_path, file_names(i).name));
    
    % 获取数据矩阵
    data = nii.img;
    
    % 提取需要相加的component
    components_data = data(:,:,:,components);
    mask = sum(components_data, 4);
    
    % 将mask数据存储到mask_data矩阵中
    mask_data = mask_data + mask;
end

% 将mask_data矩阵写入一个新的NIFTI文件中，保存为mymask.img
nii_mask = nii_template;
nii_mask.img = mask_data;
nii_mask.hdr.dime.datatype = 16;    % 修改数据类型为float32
nii_mask.hdr.dime.bitpix = 32;    % 修改位深为32位
save_nii(nii_mask, 'mymymask.img');